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This  thesis  describes  a  theoretical  and  experimental  study  of  the 
effects  of  medium  inhomogeneities  on  the  performance  of  a  parametric 
receiver. 

A  review  is  made  of  the  basic  principles  of  parametric  receiver 
operation  and  of  wave  propagation  in  an  inhomogeneous  medium. 

A  theoretical  analysis  is  presented  for  the  case  of  a  signal  wave  source 
located  on  the  axis  of  the  parametric  receiver  in  a  weakly  inhomogeneous 
medium.  Expressions  are  developed  that  predict  the  level  of  amplitude  fluc¬ 
tuations  in  the  second-order  pressure  at  the  hydrophone  of  the  parametric 
receiver.  Both  collimated  and  spherically  spreading  pump  waves  are 
considered. 

The  experimental  study  reported  in  this  thesis  was  conducted  in  a  model 
tank  in  which  an  array  of  immersion  heaters  and  a  perforated  screen  produced 
a  thermal  microstructure.  Measurements  were  made  of  the  coefficient  of 
amplitude  variation  for  the  signal,  pump,  and  interaction  frequency  waves 
associated  with  both  a  nearfield  and  a  farfield  parametric  receiver. 

Results  are  compared  to  theory  and  discussed, 
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CHAPTER  1.  INTRODUCTION 

The  subject  considered  in  this  thesis  is  the  performance  of 
a  parametric  acoustic  receiving  array  in  an  inhomogeneous  medium.  Develop¬ 
ment  of  this  subject  requires  the  synthesis  of  two  branches  of  acoustics: 
(1)  nonlinear  acoustics,  from  which  the  concept  of  the  parametric  receiver 
has  been  developed,  and  (2)  wave  propagation  in  an  inhomogeneous  medium. 
Both  of  these  topics  have  received  considerable  study  in  recent  years,  and 
it  is  reasonable  to  introduce  the  topics  separately  before  considering 
their  interrelation.  In  this  chapter  an  introduction  is  given  to  nonlinear 
acoustics  and  to  the  effects  of  medium  inhomogeneities.  Following  this  is 
an  introduction  to  the  present  investigation,  including  a  description  of 
the  aims  and  the  organization  of  this  investigation. 

1.1  Nonlinear  Acoustics  and  Parametric  Arrays 

It  has  been  known  theoretically  since  Stokes'  analysis1  in  1848  that 
an  acoustic  wave  distorts  as  it  propagates  through  a  fluid.  This  distor¬ 
tion  occurs  because  the  propagation  speed  of  a  sound  wave  is  a  function  of 

2 

the  particle  velocity;  that  is, 

=  c0  +  u(x,t)  +  y  |  u(x,t)  ,  (1.1) 


where 


^  is  the  propagation  speed, 

cq  is  a  constant  equal  to  the  sound  speed  for  waves  of  infinitesimal 
amplitude, 

—  is  a  constant  that  characterizes  the  nonlinearity  of  the  pressure- 
A  density  relation,  [defined  in  Eq.  (1.2)],  and 


u  is  the  particle  velocity  in  the  x  direction. 

It  can  be  seen  from  Eq.  (1.1)  that  there  are  two  effects  which  contribute 


1 


2 

to  the  variability  of  the  propagation  speed.  One  effect  is  'convection' 
which  arises  because  the  particles  of  the  supporting  medium  are  moving  in 
the  same  direction  as  the  propagating  wave.  convection  is  taken  into  ac¬ 
count  in  Eq.  (1.1)  by  the  term  u(x,t),  whicli  indicates  that  the  propaga¬ 
tion  speed  is  linearly  dependent  on  the  particle  velocity.  A  second 

cause  of  variations  in  propagation  speed  is  die  to  nonlinearity  in  the 

3 

pressure -density  relation.  The  relation  between  acoustic  pressure  and 
excess  density  in  a  liquid  mav  he  wrii:t«n  is’ 


2 


where 

P  =  Pj  “  Pn  is  acoustic  pressure, 

P  =  p„,  -  P  is  excess  density, 

1  o  J 


S  is  entropy,  and 

the  subscripts  zero  and  T  denote  equilibrium  and  total  values, 
respectively. 

The  quantity  B/A  is  a  measure  of  the  nonlinearity  of  a  fluid.  In  Eq.  (1.1) 

the  effect  of  this  nonlinearity  on  the  propagation  speed  is  accounted  for 

1  B 

by  the  term  —  u(x,t).  For  water,  B/A  has  a  value  of  the  order  of  5.2 

3 

while  for  air,  B/A  approximately  equals  0.4.  It  can  be  seen  by  substitu¬ 
ting  these  values  of  B/A  into  Eq.  (1.1)  that  convection  is  the  dominant  ef¬ 
fect  in  distorting  airborne  sound  waves,  while  fluid  nonlinearity  domi¬ 
nates  in  the  distortion  of  sound  waves  in  water. 

As  a  result  of  the  nonlinear  character  of  acousLie  wave  propagation, 
when  two  waves  propagate  simultaneously  in  a  fluid,  interaction-frequency 

3 

components  are  generated.  Difference  frequency  tones  have  been  observed 


3 


6_q 

by  musicians  since  1745,  and  it  was  the  observation  of  a  difference 

9  10 

frequency  sound  beam  that  led  Westervelt  to  formulate  the  concept  of  the 
parametric  array  in  1960.  This  concept  can  be  summarized  as  follows.  If 
pressure  waves*  of  frequencies  f^  and  insonify  a  common  region  of  fluid, 
there  arise  in  that  region  secondary  sound  sources  of  frequency  |f^  +  f ^ | • 
These  secondary  sources  act  as  an  array  which  can  be  used  to  transmit  or 
receive  directive  beams  of  sound. 

In  the  transmitting  application,  the  'primary  waves'  (i.e.,  the  waves 
of  frequency  f^  and  f^)  are  launched  from  a  common  transducer,  and  the  ef¬ 
fective  array  length  is  controlled  by  attenuation  in  the  medium.  Such  a 
parametric  transmitter  can  be  used  to  generate  a  directive  beam  with  negli¬ 
gible  sidelobe  levels11  at  frequencies  | f ^  ±  f^l. 

In  the  receiving  application,  one  primary  wave  is  usually  a  low 
frequency  acoustic  signal  arriving  at  the  parametric  receiver  from  some  dis¬ 
tant  source.  This  signal  interacts  with  a  locally  generated  'pump'  wave  to 
form  sum  and  difference  frequency  sources  in  the  region  common  to  both 
waves.  The  radiation  from  these  sources  is  detected  by  a  hydrophone  placed 
in  the  pump  beam.  The  array  length  in  this  instance  is  determined  by  the 
separation,  L,  between  the  pump  transducer  and  hydrophone.  The  parametric 
receiver  acts  as  a  conventional  end-fire  linear  array  of  length  L,  operating 
at  the  frequency  of  the  signal  to  be  detected.  The  advantages  of  direc¬ 
tivity  offered  by  an  end-fire  array  (i.e.,  reduction  of  interference  due  to 
ambient  noise  and  multipath  signals)  are  achieved  by  the  parametric  receiver, 
but  only  two  transducers  are  required  to  'construct'  the  array. 

Further  discussion  of  the  basic  principles  of  the  parametric  acoustic 
receiving  array,  as  well  as  an  outline  of  the  history  of  its  development, 
is  given  in  Chapter  2. 

*The  waves  are  assumed  to  be  such  that  their  propagation  vectors  are  at 
angles  of  less  than  90°  to  one  another. 
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1.2  Effects  of  Medium  Inhomogeneities 

The  variability  of  sound  transmission  in  the  sea  is  a  topic  that  has 

received  a  great  deal  of  study  in  the  last  30  years.  There  is  a  variety  of 

mechanisms  that  can  cause  fluctuations  in  the  amplitude  and  phase  of  an 
12 

acoustic  signal.  For  example,  if  an  acoustic  source  or  receiver  is 
mounted  on  a  ship,  then  irregular  motion  of  the  platform  can  give  rise  to 
fluctuations.  If  surface  reflections  contribute  to  the  received  signal, 
then  changes  in  the  surface  due  to  wave  motion  will  cause  phase  and  ampli¬ 
tude  variations.  Even  when  the  source  and  receiver  are  fixed  and  there 
are  no  surface  reflections  involved,  signal  fluctuations  will  still  occur. 

Fluctuations  having  frequencies  of  the  order  of  cycles  per  hour  to  cycles 

13 

per  day  will  result  from  internal  waves  and  internal  tides,  while  fluc¬ 
tuations  with  frequencies  in  the  range  of  cycles  per  second  to  cycles  per 
minute  will  result  from  scattering  by  random  inhomogeneities  in  the  sea.^ 
Although  these  inhomogeneities  may  be  due  to  turbulence,  thermal  micro- 
structure,  air  bubbles,  or  biological  matter,  they  may  all  be  treated  as 
variations  in  the  sound  velocity  (or  refractive  index)  of  the  medium.  The 
sea  can  then  be  crudely  modelled  as  a  medium  containing  'patches'  of  vari¬ 
able  refractive  index  which  scatter  a  propagating  acoustic  wave.  Usually 
the  principal  cause  of  scattering  is  the  thermal  microstructure,  as  we  will 
see  in  Chapter  3. 

Numerous  theoretical  studies  have  appeared  that  predict  the  amount  of 

amplitude  and  phase  fluctuations  in  an  acoustic  wave  when  the  statistical 

nature  of  the  thermal  microstructure  is  known  (see,  for  example,  refer- 
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ences  ).  Some  results  of  these  theoretical  studies  which  are  pertinent 
to  the  present  investigation  will  be  discussed  in  Chapter  3. 

As  discussed  in  the  previous  section,  there  are  three  acoustic  waves 
associated  with  the  parametric  receiver:  the  incoming  'signal'  wave,  the 
pump  wave  with  which  the  signal  interacts,  and  the  interaction  frequency 


wave  which  arises  from  the  nonlinear  interaction.  It  is  reasonable  to 
expect  that  each  of  these  waves  would  be  scattered  in  an  inhomogeneous  me¬ 
dium,  and  consequently,  there  would  be  three  components  of  fluctuation  in 
the  output*  of  the  parametric  receiver.  These  fluctuations  can  have  a 
significant  effect  on  the  performance  of  a  parametric  receiver  in  detecting 
and  resolving  low-level  acoustic  signals.  In  anticipating  practical  appli¬ 
cations  in  the  ocean,  it  is  therefore  desirable  to  know  the  effects  that 
medium  inhomogeneities  have  on  the  performance  of  a  parametric  receiver. 

It  is  the  aim  of  the  investigation  reported  in  this  thesis  to  determine 
these  effects. 

1.3  An  Overview  of  the  Thesis 

These  introductory  considerations  of  the  parametric  receiver  and  of 
the  effects  of  medium  'inhomogeneities  are  extended  in  Chapters  2  and  3,  re¬ 
spectively.  Both  of  these  chapters  contain  a  literature  survey  and  a  review 
of  basic  principles  for  their  respective  topics. 

In  Chapter  4,  the  analytical  methods  described  in  Chapters  2  and  3  are 
used  to  develop  a  theory  which  predicts  amplitude  fluctuations  in  the  inter¬ 
action  frequency  pressure  at  the  hydrophone  of  the  parametric  receiver.  It 
is  assumed  in  this  analysis  that  the  amplitude  and  phase  fluctuations  are 
small  compared  to  unity,  and  that  the  signal  source  is  located  on  the  main 
beam  of  the  parametric  receiver.  With  these  assumptions,  integral  expres¬ 
sions  are  derived  for  the  coefficient  of  amplitude  variation  for  the  cases 
of  collimated  planar  and  spherically  spreading  pump  waves.  Approximate  so¬ 
lutions  to  these  integral  expressions  are  obtained  by  making  a  number  of 
assumptions  regarding  the  spatial  correlation  functions  contained  in  these 
expressions. 

*The  parametric  receiver  'output*  may  be  taken  to  be  the  sum  or  difference 
frequency  pressure  at  the  hydrophone. 
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An  experimental  investigation  is  reported  in  Chapter  5.  The 
experiments  were  conducted  in  a  model  tank  in  which  an  array  of  immersion 
heaters  and  a  screen  mesh  were  used  to  generate  a  thermal  microstructure. 
Measurements  are  discussed  in  Chapter  5  that  were  made  to  determine  the 
coefficient  of  amplitude  variation  for  the  signal,  pump,  and  interaction  fre¬ 
quency  waves  of  a  model  parametric  receiver.  These  results  are  then  compared 
to  theoretical  predictions  based  on  the  analysis  developed  in  Chapter  4. 

In  the  final  chapter,  Chapter  6,  the  results  of  this  study  are 
summarized  and  discussed. 


CHAPTER  2. 


THE  PARAMETRIC  ACOUSTIC  RECEIVING  ARRAY 


2.1  Introduction 

The  purpose  of  this  chapter  is  to  explain  the  basic  principles  of  the 
parametric  acoustic  receiving  array,  and  to  outline  the  history  of  its  de¬ 
velopment.  The  theoretical  development  of  the  parametric  acoustic  array  was 
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based  upon  Lighthill's  study  *  of  sound  generated  aerodynamically . 

Westervelt  applied  Lighthill's  results  to  the  problem  of  scattering  of 
22  23  10 

sound  by  sound  '  and  later  extended  this  work  to  the  parametric  array. 
In  formulating  the  theory  of  the  parametric  array,  Westervelt  developed  an 
inhomogeneous  wave  equation  that  describes  the  second-order  sound  field.  A 
brief  sketch  of  the  development  of  this  wave  equation,  as  well  as  a  discus¬ 
sion  of  the  physical  significance  of  Westervelt' s  results,  will  be  given  in 
the  next  section  of  this  chapter.  Solutions  of  the  second-order  wave  equa¬ 
tion  for  various  configurations  of  the  parametric  receiving  array  will  then 
be  discussed.  Finally,  practical  matters  such  as  signal  processing,  shad¬ 
owing,  finite  amplitude  effects,  and  the  effects  of  noise  on  the  performance 
of  the  parametric  receiver  are  considered. 

2.2  Basic  Theoretical  Concepts 

We  will  begin  discussion  of  the  parametric  acoustic  array  by  turning 
briefly  to  some  of  the  first  principles  of  acoustics.  Acoustic  wave  propa¬ 
gation  in  a  lossless  medium  is  governed  by  the  following  equations  of  fluid 
mechanics: ^ 

Continuity  Equation 


3PT 

t r+v '  Wu)  - 0 


(2.1) 
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Momentum  Equation 


(PTU)  +  PTu(V-u)  +  (u*V)pTu  +  Vp  =  0 
. sentropic  Equation  of  State 


P  =  p(p)s 


(2.2) 


where 

is  the  total  density, 

V  is  the  vector  differential  operator  del, 
u  is  the  particle  velocity, 

p  is  the  acoustic  pressure,  and 

the  subscript  S  denotes  constant  entropy. 

The  isentropic  equation  of  state  may  be  expanded  in  a  Taylor  series  about 
the  ambient  density  of  the  fluid,  giving 


.  2,  .  1  Co  /B) 

PT  ‘  »o  +  co  <PT  -  V  +  2  T^b.} 


(PT  -  !>„) 


(2.3) 


where  the  terms  higher  than  second-order  have  been  neglected.  It  may  be 
noted  that  Eqs.  (2.3)  and  (1.2)  are  identical,  except  for  a  change  of 
variables. 

The  field  variables  p^,  u,  and  p  appearing  in  these  fundamental 
equations  can  be  expanded  in  series  such  that 


PT  =  PQ  +  pi  +  P, 


+  ...  +  p 


U  =  ll  +  U-+U-  +  ...+U 

o  1  2  n» 


lr  =  Po  +  pi  +  P2  +  •••  +  Pn 


and 


(2.4) 
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The  physical  meaning  of,  for  example,  the  various  pressure  terms  in  Eqs. 

(2.4)  is  as  follows.  The  zeroth-order ,  or  static  component,  p  ,  is  the 
pressure  of  the  fluid  in  the  absence  of  a  sound  wave.  The  first-order  term, 
p^,  is  the  pressure  associated  with  a  wave  propagating  linearly  through  the 
fluid.  The  higher  order  components  are  due  to  the  nonlinear  nature  of  acous¬ 
tic  wave  propagation,  as  discussed  in  Section  1.1.  Under  the  limitations  of 
perturbation  theory  for  which  the  acoustic  signals  are  sufficiently  small 
that 

P0  »  PX  »  P2  »  •••  »  Pn  • 
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terms  higher  than  second  order  are  neglected.  Neglecting  these  terms 
is  effectively  the  same  as  assuming  that  no  nonlinear  interactions  occur 
beyond  those  which  give  rise  to  the  second-order  field;  i.e.,  the  second- 
order  radiation  is  assumed  to  propagate  linearly. 

First  and  second-order  wave  equations  can  be  derived  by  substituting 
the  series  expressions  in  Eq.  (2.4)  into  the  fundamental  equations,  Eqs. 
(2.1)-(2.3),  and  then  retaining  only  terms  of  first  or  second  order,  re¬ 
spectively.  The  parametric  acoustic  array  utilizes  the  second-order  field 

variables,  so  the  second-order  wave  equation  is  of  particular  interest,  and 

,  .  ,  10,26 
is  given  by 


□  P2  =  - 


.2  2 

B  3  P1 


p  c  A  3t2 
o  o 


(2.5) 


where 


o  2  2  ^ 

|  |  =  V  -  — j  — -  is  the  D ' Alembertian  operator,  and 

co  at2 


6  =  1  +  2l 
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This  inhomogeneous  wave  equation  has  a  source  term  on  the  right  hand  side 
that  is  a  function  of  the  first-order  pressure  field  and  some  constants. 
Apparently  to  emphasize  the  physical  meaning  of  this  wave  equation, 
Westervelt  rewrote  Eq.  (2.5)  in  terms  of  the  source  strength  density,  q,  as 
follows : 


where 


□  P2  -  -  »o-j? 


3P. 


2  4  3t 


P  c 
o  o 


(2.6) 


Equation  (2.6)  has  the  same  form  as  the  linear  (first-order)  wave  equation 
in  which  there  are  real  mass  sources  present,  the  sources  having  source 
strength  density,  q.  By  analogy  with  such  a  linear  wave  equation,  Eq.  (2.6) 
describes  the  propagation  of  second-order  acoustic  waves  that  are  generated 
by  the  nonlinear  interaction  of  first-order  acoustic  waves.  The  sources 
that  generate  the  secondary  pressure  field  are  not  real,  but  are  'virtual' 
sources  produced  by  the  nonlinear  properties  of  acoustic  wave  propagation. 
These  sources  are  distributed  throughout  the  region  of  space  in  which  the 

nonlinear  interaction  of  the  primary  field  occurs.  The  general  solution  for 

27 

the  pressure  p2  in  Eq.  (2.6)  is 


P2(r’C)  =  ^  J?^q(ro'T)dv 

V 

where 

-j 

r  is  the  observation  point, 
r  is  the  source  point. 
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T  =  t  -  —  IS  the  retarded  time, 
c 

o 

V  is  the  interaction  volume,  and  the  geometry  of  Fig.  2.1 
applies. 

This  result  can  be  summarized  qualitatively  as  follows.  Whenever  some 
first-order  pressure  p^  is  present  at  a  point  in  the  medium,  a  second-order 
virtual  source  will  arise  at  that  point  having  a  source  strength  qdv.  The 
second-order  pressure  at  an  observation  point  r  will  be  the  sum  of  all 
sources  in  the  interaction  volume;  this  pressure  is  given  by  Eq.  (2.7).  In 
the  next  few  sections,  we  will  discuss  specific  solutions  for  p^  for  given 
first-order  pressure  fields. 


2.3  The  Parametric  Array 

Two  practical  applications  for  the  second-order  sound  field  described 
by  Eq.  (2.7)  were  proposed  by  Westervelt. 10  He  noted  that  if  two  plane 
harmonic  waves  of  frequencies  and  w2  are  projected  by  a  common  trans¬ 
ducer,  an  array  of  virtual  sources  will  be  established  in  the  beam  of  the 
transducer,  radiating  at  sum  and  difference  frequencies  ±  The  array 

would  be  phased  in  an  end-fire  manner,  and  its  radiation  pattern  would  have 

negligible  sidelobe  levels.  This  concept  of  a  parametric  transmitting  array 

28 

was  confirmed  experimentally  by  Beilin  and  Beyer,  and  led  to  further 

29  30  31 

studies  of  this  application  by  Berktay,  Berktay  and  Smith,  Hobaek, 

32  33  11  25  34 

Zverev  and  Kalachev,  Muir  and  Blue,  Muir,  and  Smith.  ’ 

In  addition  to  its  use  as  a  transmitter,  Westervelt  also  suggested 

using  the  parametric  array  as  an  acoustic  receiver.  In  this  application 

some  distant  signal  would  interact  nonlinearly  with  a  locally  generated 

signal.  The  second-order  pressure  arising  from  this  interaction  would 

then  be  detected  by  a  hydrophone  placed  in  the  vicinity  of  the  locally  gen- 
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erated  wave.  Several  studies  were  also  prompted  by  this  possibility. 


The  findings  of  some  of  these  studies  will  be  discussed  in  the  following 
sections . 
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2.3.1  The  Nearfiold  Receiving  Array 

As  an  illustration  of  the  parametric  receiving  array,  consider  the 
arrangement  shown  in  Fig.  2.2.  A  pump  transducer  of  dimensions  2b  x  2b  is 
situated  at  the  origin  and  radiates  a  wave  of  frequency  along  the  x-axis. 
A  signal  wave  with  frequency  u>  propagates  in  a  direction  given  by  the  unit 
vector 


K  =  i  cosO  +  j  sinO  ,  (2.8) 

where  i  and  j  are  unit  vectors  in  the  x  and  y  directions,  respectively. 

The  second-order  pressure  produced  by  interaction  of  the  pump  and 

signal  waves  is  detected  by  a  hydrophone  placed  at  x  =  L.  It  is  assumed 

that  the  detector  is  in  the  nearfield  of  the  pump  transducer,  so  that  the 
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pump  radiation  may  be  approximated  by  a  plane  wave.  The  signal  source  is 

assumed  to  be  located  so  distantly  that  the  signal  wave  is  planar  in  the 

vicinity  of  the  pump.  With  these  considerations,  the  pump  and  signal  pres¬ 
sures  are  described  by 

j  (w  t  -  k  x)  -ax 

p  (x,t)  =  P  e  p  P  e  P  (2.9) 

P  P 

j  (a)  t  -  k  x)  -a  R 

pg(x,t)  =  Pg  e  S  S  e  (2.10) 

P  is  the  pressure  amplitude  of  the  pump  wave  in  the  nearfield 
P 

of  the  pump  transducer, 

Pg  is  the  pressure  amplitude  of  the  signal  wave  at  the  origin, 

a  and  a  are  attenuation  coefficients  for  the  pump  and  signal 
P  s 

waves,  respectively,  and 
R  =  x  cos0  +  y  sin© 

In  writing  Eqs.  (2.9)  and  (2.10),  it  has  been  assumed  that  the  attenuation 
coefficients  are  constants  because  the  waves  are  each  of  single  frequency. 


Pump: 


Signal: 


where 
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Also,  the  phase  of  the  signal  wave  has  arbitrarily  been  set  to  zero  at  the 
origin. 

The  total  first-order  pressure  at  any  point  in  the  pump  column  (i.e., 
the  region  indicated  by  dashed  lines  in  Fig.  2.2)  is 


P1  = 


+  P, 


(2.11) 


The  second-order  pressure  detected  by  the  hydrophone  at  x  *  L  can  be 
calculated  by  substituting  Eq.  (2.11)  into  Eqs.  (2.7)  and  (2.8).  If  only 
the  sura  and  difference  frequency  terms  are  retained,  the  result  is 


P+OL)  = 


W±26PPPS 


L  b  b 

Iff  ^ 


-j (k  x  ±  k  R)  -(a  x  +  a  R) 
p  s  p  s 

e 


4irp  c 

o  o  o  -b  -b 


r(ai  *  jy/q—  z)2  +  v2  +  z2 

/(L  -  Z  ) 2  +  y2  +  Z2 


dxdydz 


(2.12) 


where  w  =  w  ±  w  , 

i  p  s 

a+  is  the  attenuation  coefficient  at  frequency  w+  (a  constant), 
k+  is  the  acoustic  wavenumber  at  frequency  w+,  and  the  time  de¬ 
pendence  has  been  suppressed. 
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This  integral  is  not  easy  to  evaluate  directly,  but  Berktay  has 

devised  a  method  for  simplifying  the  integration  which  will  now  be  adopted. 

The  interaction  region  in  front  of  the  pump  transducer  can  be  regarded  as  a 

series  of  elemental  wafers  of  thickness  dx  and  cross-sectional  dimensions 

2b  x  2b  (see  Fig.  2.3).  Each  of  these  wafers  is  insonified  by  the  first- 

order  sound  field  and  will  radiate  at  the  frequency  u>+.  These  wafers  will 

produce  nearfield,  planar  radiation  at  a  range  less  than  Rq  from  the  wafer, 

and  farfield,  spherically  spreading  radiation  at  ranges  greater  than  R  . 

45 

Using  Freedman's  description  of  radiation  from  a  square  piston,  Berktay 
38 

and  Al-Temimi  choose  R  to  be 
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y 


L  -R , 


AREA  (2b)* 


x  =  L 


PUMP 


dx 


FIGURE  2.3 

ELEMENTAL  WAFERS  IN  THE  INTERACTION  REGION 
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R  =  1.3  b/X ,  ,  (2.13) 

O  3 

where  X  is  the  acoustic  wavelength  at  frequency  w  • 

s  s 

The  detected  second-order  pressure  may  be  considered  to  have  two 

components:  (1)  nearfield  radiation  from  wafers  within  a  distance  Rq  of  the 

detector  and  (2)  farfield  radiation  from  wafers  at  a  distance  greater  than 

R  from  the  detector.  For  this  illustration  assume  that  L  <  R  so  that 
o  o 

all  wafers  are  in  the  nearfield  of  the  detector.  With  the  further  assump¬ 
tion  that  absorption  along  the  v-axis  is  negligible  compared  to  that  in  the 

x  direction,  so  that  agR  can  be  approximated  by  a^x  cos6,  Berktay  and 
38 

Al-Temini  evaluate  Eq.  (2.12)  with  the  result 

-jio  P  Pf5Lr  r-i-tim 

PN(L,0)  =  - - ~  p  •  [exp  -(a±  +  jk±)L,J  •  exp[+ JmJ  *  ^ — 

2P°C°  (2.14) 


where 

the  subscript  N  indicates  that  the  result  is  valid  only  when  the 
detector  is  in  the  nearfield  of  all  elemental  wafers,  and 
M  =  j  kgL(l  ~  cos9)'  (2.15) 

The  pressure  detected  by  the  hydrophone  can  be  expressed  in  terms  of 
the  on-axis  pressure  amplitude,  P^(L,0),  and  the  directivity,  D(0),  as 
follows : 


where 


pNa,0)  =  pNa,o)  •  d (0 ) 

to  P  P  0L 

pNa,°)  =  ~  P  S3-  exp(-a±L) 


D(0) 


sin  M 
M 


(2.16) 

(2.17) 


and 


(2.18) 
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Experimental  confirmation  of  the  validity  of  Eqs.  (2.17)  and  (2.18)  is 

37  38  47 

reported  by  Berktay  and  Al-Temini.  ’  ’ 

The  directivity  function  is  identical  to  that  of  a  continuous  end-fire 
array  of  length  L,  realized  from  elements  operating  at  frequency  w  .  It  is 
this  directional  property  of  the  parametric  receiver  that  makes  it  an  at¬ 
tractive  device  for  the  reception  of  underwater  signals.  It  offers  the  same 
capabilities  as  the  conventional  end-tire  array  as  regards  directional  de¬ 
tection  of  signals;  yet  the  parametric  device  requires  only  two  transducers 
in  its  construction. 

These  conclusions  have  been  based  on  results  obtained  assuming  that 
the  hydrophone  is  placed  in  the  nearfield  of  the  pump  transducer.  In  appli¬ 
cations  requiring  high  directivity,  and  consequently  large  values  of  L,  it 
is  likely  that  the  hydrophone  will  be  placed  well  in  the  farfield  of  the 
pump.  The  next  section  extends  the  discussion  to  this  situation. 

2.3.2  The  Farfield  Receiving  Array 
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In  1972,  Barnard  ex  a_l.  considered  a  case  in  which  the  hydrophone  is 
placed  in  the  farfield  of  the  pump  transducer.  Their  analysis  assumed 
interaction  between  a  plane  signal  wave  and  a  spherical  pump  wave  produced 
by  a  baffled  circular  piston.  Prediction  of  the  interaction  frequency  pres¬ 
sure  was  given  by  a  numerical  solution  to  Eq.  (2.12).  The  results  of  the 
integration  were  tested  against  experiment  for  the  sum  frequency,  with  good 
agreement  between  theory  and  experiment. 

A  closed-form  solution  for  parametric  reception  with  spherically 
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spreading  pump  waves  is  given  by  Berktay  and  Shooter.  The  analysis  is 
similar  to  that  for  the  nearfield  receiving  array,  but  rather  than  having 
planar  elemental  wafers,  the  interaction  frequency  sources  are  assumed  to 
be  cophasal  along  the  spherical  wavefronts.  The  received  pressure  at  the 
interaction  frequency  w+  can  be  written 


l‘l 


-w  i’’P 


UJ  A  A  ^  \J  r 

PF(L,0)  =  — - — E-—  exp  +(a±  +  jk±)L~jM 
2n  c ;  L  J 


D(0) 


2p  c 
o  o 


(2.19) 


where 

p  (L,0)  is  the  pressure  of  the  interaction  frequency  signal  at  x  =  L, 

F 

for  the  case  of  spherical  pump  waves, 

P'  is  the  on-axis  farfield  pump  pressure  amplitude  referred  to  1  meter, 

P 

D(0)  =  sin  M/M  is  the  directivity  function,  and 

M  =  k  L(1  -  cos0). 

2  s 

In  comparing  results  for  the  collimated  pump  [Eq.  (2.14)]  and  the 
spherically  spreading  pump  [Eq.  (2.19)],  it  can  be  seen  that  the  directivity 
patterns  for  the  two  cases  are  identical.  The  pressure  amplitude  for  the 
case  of  the  spherically  spreading  pump  is,  with  the  exception  of  absorption, 
independent  of  array  length,  whereas  for  the  collimated  pump,  the  pressure 
amplitude  increases  with  L. 

The  directivity  function  of  the  pump  transducer  does  not  enter  into  the 

result  for  the  farfield  parametric  receiver,  but  it  has  been  assumed  in  the 

analysis  that  the  pump  beamwidth  is  less  than  the  beamwidth  of  the  parametric 
43 

receiver.  The  directivity  function  of  the  parametric  receiver  for  both 
the  collimated  and  spherically  spreading  pump  is  independent  of  pump  fre¬ 
quency,  and  is  the  same  whether  the  sum  or  difference  frequency  signal  is 
used . 

A  3 

Berktay  and  Shooter  compared  the  closed-form  solution  [Eq.  (2.19)] 

39 

to  results  obtained  by  the  numerical  integration  of  Barnard  et  al.  and 

found  that  agreement  between  the  predictions  was  good.  Computed  results  of 

Eq.  (2.19)  were  also  found  to  be  in  good  agreement  with  experimental  mea- 
42  43 

surements.  ’  These  results  show  that  the  amplitude  of  the  interaction 


frequency  signal  component  can  be  made  approximately  equal  to  the  original 
signal  amplitude  by  proper  choice  of  pump  source  level  and  array  length  so 
that  no  loss  in  signal  level  is  experienced  in  parametric  reception. 
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Truchard  ’  solved  Che  problem  of  a  parametric  receiving  array  with 
an  arbitrarily  shaped  planar  piston  pump  as  follows.  He  first  analyzed  the 
interaction  of  a  high  frequency  spherical  wave  from  a  point  source  and  a 
low  frequency  plane  wave.  The  solution  for  pump  transducers  of  two  or  three 
dimensions  was  then  obtained  by  summing  the  point  source  solution  over  the 
active  face  or  volume  of  the  pump  transducer.  In  this  manner  Truchard 
analyzed  the  following  situations:  (1)  a  point  pump  source  with  a  line 
hydrophone,  (2)  a  rectangular  pump  with  a  point  hydrophone,  and  (3)  a  cir¬ 
cular  piston  pump  with  a  point  hydrophone.  Truchard  also  studied  these  sit^_ 
uations  experimentally  and  obtained  very  good  agreement  between  his  theoret¬ 
ical  and  experimental  results. 

2.4  Practical  Considerations 
2,4.1  Signal  Processing 

In  practical  applications  there  are  factors  which  make  detection  of  the 

interaction  frequency  signal  difficult.  Since  u>  <<  m  ,  the  interaction 

s  p 

frequency  oi  ±  u>  is  very  close  to  the  pump  frequency  iu  .  Also  the  pump 
p  s  p 

amplitude  P  is  generally  much  greater  than  either  of  the  amplitudes,  P  and 

r  S 

P+.  This  situation  is  illustrated  by  the  sketch  of  Fig.  2.4.  The  quantity 

/  2 

ex  shown  in  the  figure  is  the  'pump  excess,'  defined  by  Berktay  and  Muir 

to  be  the  ratio  of  the  pump  pressure  at  the  hydrophone  to  the  amplitude  of 

A  2 

the  interaction  frequency  signal.  For  the  case  of  spherical  pump  waves 


e 

x 


=  2p  c  / (u 

0  0  p 
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s  s 


(2.20) 


The  signal  processing  problem  is  to  reject  the  high  amplitude  pump  signal 
and  retain  one  or  both  of  the  interaction  frequency  signals.  Since  a  pump 
excess  on  the  order  of  160  dB  may  be  encountered  in  practice, ^  careful  de¬ 
sign  of  the  signal  processing  system  may  be  required  in  applications  of  the 
parametric  receiver. 


relative  pressure  amplitude 


FIGURE  2.4 

SKETCH  OF  THE  FILTERING  PROBLEM 
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There  are  three  basic  types  of  signal  processing  systems  that  have 
been  used  with  the  parametric  receiving  array:  the  band  pass  receiver,  the 
carrier  rejection  or  band  elimination  receiver,  and  the  phase-locked  loop 
receiver.  The  band  pass  receiver  is  relatively  simple  to  design  and  con¬ 
struct  and  is  adequate  in  applications  where  the  pump  excess  is  low  and 
the  interaction  frequency  is  well  separated  from  the  pump  frequency.  This 
type  of  receiver  was  used  in  the  present  study,  and  is  described  in  Chapter 
5. 

For  applications  involving  a  high  pump  excess  and  low  signal  frequencies, 
the  band  elimination  receiver  is  suitable.  This  type  of  receiver  may  be 
described  using  the  block  diagram  of  Fig.  2.5.  The  design  uses  a  cascade  of 
two  crystal  notch  filters  and  two  preamplifiers  to  reduce  the  level  of  the 
carrier,  and  to  amplify  the  interaction  frequency  (i.e.,  sideband)  compo¬ 
nents.  The  signal  is  then  split  into  its  in-phase  and  quadrature  components, 

and  these  components  are  processed  to  give  the  upper  and  lower  sideband  sig- 
48 

nals.  Rohde  et  al.  report  this  method  of  signal  processing  to  be  useful  in 
a  parametric  receiver  with  carrier  to  sideband  ratios  approaching  180  dB. 

A  third  type  of  signal  processing  system  is  the  phase-locked  loop  (PLL) 

receiver,  shown  in  Fig.  2.6.  This  system  operates  as  a  closed  loop  servo 

49 

system  with  phase  as  the  controlled  variable.  The  input  of  the  PLL  re¬ 
ceiver  is  connected  to  the  hydrophone  of  the  parametric  receiving  array. 

With  the  voltage  controlled  oscillator  operating  nominally  at  the  pump  fre¬ 
quency,  and  with  the  loop  bandwidth  less  than  the  lowest  modulation  frequency 
of  interest,  the  error  voltage  r(t)  contains  the  demodulated  signal.  This 

type  of  receiver  is  useful  in  minimizing  the  effects  of  low  frequency  noise, 
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such  as  due  to  transducer  vibration.  '  This  point  will  be  further  dis¬ 
cussed  in  Section  2.4.5. 

2.4.2  Shadowing 

Any  time  that  an  obstacle  significantly  large  with  respect  to  wavelength 
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is  placed  in  the  path  of  a  propagating  wave,  shadows  may  form  in  the  region 
behind  the  obstacle.  In  the  case  of  the  parametric  receiver,  the  pump 
transducer  and  its  mounting  obstruct  the  incoming  signal  wave  and  form  a 
shadow  in  the  array  interaction  region.  This  shadow  has  the  effect  of  trun¬ 
cating  the  array  length,  with  consequent  effects  on  the  sensitivity  and  the 

directivity  function  of  the  array. 
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Al-Temimi  ’  has  done  a  theoretical  and  experimental  study  of  the 
effects  of  shadows  on  the  performance  of  parametric  receivers.  One  impor¬ 
tant  result  of  his  work  is  a  description  of  the  condition  for  which  acoustic 
shadows  cause  no  appreciable  deterioration  of  array  performance  and  hence 
may  be  neglected.  This  condition  is  that  the  shadow  length  d  be  kept  below 
5%  of  the  array  length  L;  i.e., 


d  <  0.05  L,  (2.21) 

where 

d  =  1.25  b2/X  ,  (2.22) 

s 

and  2b  is  the  length  of  each  side  of  a  square  transducer. 

2.4.3  Finite  Amplitude  Effects 

It  was  noted  in  Chapter  1  that  a  finite  amplitude  acoustic  wave  distorts 
as  it  propagates,  this  distortion  being  due  to  the  dependence  of  the  phase 
velocity  on  the  particle  velocity;  i.e., 

^  =  cq  +  8  u(x,t)  ,  (2.23) 


cq  is  the  sound  speed  with  respect  to  the  fluid  particles. 


where 
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is  the  phase  (or  propagation)  velocity  of  a  given  point  on  a 
sound  waveform,  and 
u  is  particle  velocity. 

The  distortion  of  an  initially  sinusoidal  acoustic  wave  gives  rise  to 

harmonic  components.  Energy  is  transferred  from  the  fundamental  component 

of  the  wave  to  the  harmonics  as  the  wave  is  progressively  distorted,  so  the 

fundamental  is  attenuated  by  this  effect.  If  the  wave  experiences  enough 

distortion  for  a  discontinuity  or  'shock'  to  develop,  the  situation  becomes 

more  complex.  When  the  acoustic  wave  is  'weak,'  i.e.,  u  /c  <  0.1,  then  a 

o  o 

method  known  as  'weak-shock  theory'  can  be  used  to  describe  the  wave  propa- 
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gation.  ’  We  will  not  discuss  this  theory  in  any  detail,  but  merely  note 

that  when  a  shock  forms,  dissipation  may  be  assumed  to  occur  at  the  shock 

52 

front  to  the  degree  that  the  waveform  never  becomes  multivalued. 

In  general,  there  are  three  mechanisms  that  attenuate  an  acoustic  wave: 
geometrical  spreading,  absorption  by  the  medium,  and  the  'excess  attenuation' 
due  to  energy  transfer  from  fundamental  to  harmonic  components  as  the  wave 
is  distorted.  The  spreading  and  absorption  mechanisms,  by  reducing  the  amp¬ 
litude  of  an  acoustic  wave,  tend  to  reduce  the  nonlinear  distortion  process 
considerably.  For  example,  a  spherically  spreading  wave  will  have  a  value 
of  u  that  reduces  as  -p,  so  that  the  importance  of  the  second  term  in  Eq. 
(2.23)  diminishes  with  range.  Absorption  will  similarly  reduce  the  effects 
of  nonlinearity.  It  can  be  seen,  therefore,  that  excess  attenuation  will 
only  be  significant  in  the  case  of  waves  with  initially  large  amplitudes. 

The  excess  attenuation  for  a  plane  wave  propagating  through  a  lossy 
medium  is  calculated  by  Blackstock"^  in  terms  of  the  Gol'dberg  number,  T, 
and  the  normalized  range  parameter,  o.  These  quantities  are: 

f  ek 


r 


a 


(2. 24) 
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where 


u 

o 


o 


uq  being  the  particle  velocity  at  the  transducer  face, 

and 


a 


x 

gek 


(2.25) 


The  value  of  the  excess  attenuation  EXDB  in  the  fundamental  component  can 
be  read  from  curves  given  in  reference  54.  An  approximate  solution,  said 
to  be  accurate  to  within  1  dB  of  Blackstock's  results,  is  given  by^’^ 


/p\2  2 

EXDB  -  10  log 

1  4  [1  -  exp(-2a/F)  ] 

Some  effects  of  finite  amplitude  attenuation  on  the  performance  of  a 
parametric  receiver  have  been  studied  by  Berktay  and  Al-Temimi. They 
found  that  for  the  case  of  the  nearfield  array,  excess  attenuation  of  the 
pump  wave  can  bring  about  saturation  in  the  sensitivity  of  the  parametric 
receiver. 
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It  was  found  by  Al-Temimi  that  finite  amplitude  attenuation  had  little 
effect  on  the  parametric  receiver  directivity.  The  combination  of  shadowing 
and  finite  amplitude  attenuation,  both  dependent  on  frequency,  can  reduce 
the  bandwidth  of  the  receiver.  The  effect  is  particularly  significant  when 
a  great  deal  of  extra  attenuation  occurs  in  the  shadow  region  of  the  pump, 
so  that  the  pump  wave  is  well  attenuated  before  entering  the  interaction 
volume. ^ 


2.4.4  Sources  of  Noise  that  Affect  Parametric  Reception 

Any  noise  in  the  acoustic  environment  at  frequencies  to  and  to  ±  to 
J  s  p  s 

will  be  detected  by  the  parametric  receiver.  The  signal  frequency  noise 
component  (at  to^)  will  be  up-converted  along  with  the  desired  signal,  and 
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the  interaction  frequency  noise  component  (at  w  ±  w  )  will  be  detected 

p  s 

directly  by  the  receiving  transducer.  The  parametric  receiver,  by  virtue  of 

its  directivity  at  ta  ,  discriminates  against  the  low  frequency  noise.  Some 

rejection  of  the  interaction  frequency  noise  is  provided  by  the  receiving 

transducer  directivity.  An  analysis  of  the  effects  of  ambient  noise  on 
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parametric  reception  may  be  found  in  Berktay's  papers.  ’  ’  Effects  of 

ambient  noise  were  an  important  consideration  in  the  parametric  receiver  de¬ 
sign  described  by  Goldsberry ,  ^  ^  who  discusses  in  Ref.  61  an  extensive  ex¬ 
perimental  program  which  studied  the  performance  of  the  parametric  receiver 
in  a  fresh  water  lake. 

In  addition  to  noise  in  the  acoustic  environment,  electronic  noise  in 

the  pump  signal  source  will  contribute  to  the  overall  noise  level  of  the 

parametric  receiver.  This  noise  is  due  to  the  fact  that  the  pump  transducer 

does  not  radiate  at  a  single  frequency,  but  over  a  frequency  band,  and  if 

co  <<  on  ,  then  there  is  very  likely  to  be  energy  radiated  by  the  pump  at 
s  p 

co  ±(o,  which  will  be  detected  by  the  hydrophone.  This  radiation  depends 
P  s 

primarily  on  the  sideband  noise  produced  by  the  pump  oscillator,  and  must 

be  minimized  if  low-level,  low  frequency  signals  are  to  be  detected  bv  the 

parametric  receiver.  A  crystal  oscillator  developed  for  the  purpose  of 

minimizing  the  noise  of  the  pump  source  is  reported  to  have  a  sideband  noise 

level  (measured  in  a  1  Hz  band)  of  -160  dB  referenced  to  the  carrier  level, 

61  62 

measured  at  100  Hz  from  the  carrier  frequency.  ’ 

Another  kind  of  noise  may  arise  due  to  transducer  vibration.  If  the 
pump  transducer  and  hydrophone  are  mounted  such  that  they  move  relative  to 
one  another  because  of  vibration,  then  the  array  length  will  be  varied  by 
this  motion.  As  a  result  the  pump  and  interaction  frequency  signals  will 
be  phase  modulated.  This  modulation  produces  undesirable  sidebands  in  the 
hydrophone  output  which  may  appear  as  spurious  low  frequency  signals  that 
are  indistinguishable  from  acoustic  signals. ^  The  response  of  a  parametric 
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receiver  to  transducer  vibration  has  been  studied  analytically  and 
experimentally  by  Reeves  et^  ^.63,64  and  by  Lamb.^  Lamb  has  demonstrated 
the  usefulness  of  the  phase-locked  loop  receiver,  discussed  in  Section  2.4.1, 
in  minimizing  the  noise  due  to  transducer  vibration. 

2.4.5  Effect  of  Turbulence  and  Medium  Inhomogeneities 

Another  factor  which  may  affect  receiver  performance  is  the  signal 
fluctuation  due  to  turbulence  and  medium  inhomogeneities.  Both  turbulence 
and  inhomogeneities  in  the  medium  produce  random  scattering  of  an  acoustic 
signal.  The  various  scattered  components  interfere  at  an  observation  point 
to  produce  amplitude  and  phase  fluctuations  in  the  observed  signal.  For 
the  parametric  receiver  there  are  three  signal  components  which  undergo 
fluctuations:  the  incoming  low  frequency  signal,  the  pump  signal,  and  the 

interaction  frequency  signal.  In  a  severely  inhomogeneous  or  turbulent 
medium,  the  composite  effect  of  these  fluctuations  could  prove  prohibitive 
to  the  detection  of  low-level  signals.  This  is  a  problem  that  has  received 
little  attention  to  date,  and  it  is  the  principal  topic  with  which  this 
thesis  will  be  concerned. 

2.5  Summary 

In  this  chapter,  the  basic  concepts  of  the  parametric  receiver  have 
been  discussed  and  an  outline  of  the  historical  highlights  of  its  develop¬ 
ment  presented.  Of  necessity,  topics  have  been  treated  briefly.  Some 

42 

topics,  such  as  the  use  of  arrays  of  parametric  receivers,  the  phase  mod- 

41 

ulation  model  of  the  parametric  receiver,  and  the  effects  of  having  the 
signal  source  located  in  the  nearfield  of  the  parametric  receiver,65  have 
been  omitted,  as  they  do  not  directly  bear  upon  the  present  study. 

It  has  been  shown  in  this  chapter  that  the  parametric  array  makes  use 
of  the  second-order  sound  field  generated  by  nonlinear  interaction  of  a  low 
frequency  signal  wave  and  a  locally  generated  pump  wave.  A  second-order 
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wave  equation  [Kq.  (2.6)]  was  presented,  the  solution  to  which  is  dependent 
upon  the  geometry  of  the  interaction  region.  Solutions  to  the  wave  equation 
for  the  cases  of  the  nearfield  and  farfield  parametric  receiver  were  dis¬ 
cussed.  It  was  seen  from  these  solutions  that  the  parametric  receiver  has 
the  same  directivity  characteristics  of  an  end-fire  array  of  equal  length, 
realized  by  continuous  elements.  Thus  the  advantages  of  directivity  offered 
by  the  end-fire  array  are  achieved  by  the  parametric  receiver,  which  re¬ 
quires  only  two  transducers  for  its  construction.  Finally,  practical  con¬ 
siderations  regarding  signal  processing,  shadowing,  finite  amplitude  effects, 
and  the  effects  of  noise  were  discussed  briefly.  It  was  noted  that  one  re¬ 
maining  topic  for  study  regarding  practical  application  of  the  parametric 
receiver  is  the  effects  of  medium  inhomogeneities;  this  is  the  topic  for 


this  thesis. 


CHAPTER  3.  ACOUSTIC  WAVE  PROPAGATION  IN  AN  INHOMOGENEOUS  MEDIUM 

3.1  Introduction 

The  discussion  of  the  parametric  receiver  in  the  previous  chapter 
implicitly  assumed  that  the  acoustic  medium  is  homogeneous.  This  assumption 
is  usually  unrealistic  in  practice.  The  ocean,  for  example,  may  contain  one 
or  more  of  a  variety  of  inhomogeneities:  fish,  bubbles,  algae,  thermal 
microstructure,  and  turbulent  eddies.  These  inhcmogeneities  scatter  the 
energy  of  an  acoustic  wave  and,  as  they  move  about  in  the  medium,  produce 
random  fluctuations  in  the  amplitude  and  phase  of  the  transmitted  wave.  In 
many  instances  the  effects  of  fish,  bubbles,  and  algae  can  be  ignored  be¬ 
cause  these  kinds  of  scatterers  occur  in  small  quantities.  More  widespread 

are  thermal  inhomogeneities  and  turbulent  eddies  which  can  cause  significant 

14 

fluctuations  in  a  transmitted  signal.  Both  of  these  inhomogeneities  can 

be  modelled  as  regions  or  'patches'  of  variable  refractive  index. 

Thermal  patches  are  usually  the  principal  cause  of  scattering  in  the 
14 

ocean.  The  dimensions  of  these  patches  are  generally  large  compared  to  the 

acoustic  wavelength,  in  which  case,  as  a  first  approximation,  they  may  be 

treated  as  spherical  sound  sources  of  radius  a,  where  a  is  the  mean  radius 
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of  the  thermal  patches.  ’  Scattered  sound  from  the  patches  will  propa¬ 
gate  in  the  same  direction  as  the  incident  wave,  and  with  a  farfield  beam 

angle  of  1/ka  radians,  where  k  is  the  acoustic  wavenumber  (see  Fig.  3.1). 

2 

The  region  in  front  of  the  patch  out  to  a  distance  ka  is  taken  to  be 
the  'ray  region'  of  the  scattered  radiation. ^  In  this  region  the  patch  will 
behave  like  a  lens,  focusing  or  defocusing  the  scattered  rays  according  to 
whether  the  sound  velocity  in  the  patch  is  smaller  or  larger  than  the  aver¬ 
age  value,  respectively.^  In  1946,  a  'ray  theory'  was  developed  by 

Bergman^  which  predicted  that  a  wave  propagating  over  a  distance  less  than 
2 

ka  would  undergo  amplitude  fluctuations  proportional  to  the  three-halves 
power  of  the  range. 
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The  region  beyond  a  distance  ka  from  the  patch  is  the  ,'wave  region' 

2 

of  the  scattered  radiation.  At  ranges  greater  than  ka  from  a  source,  inter¬ 
ference  will  occur  between  waves  scattered  by  various  patches  in  the  medium. 
For  these  ranges  the  theory  of  Mintzer^  predicts  mean  square  amplitude 
fluctuations  in  the  acoustic  signal  which  increase  linearly  with  range. 
Mintzer's  prediction  of  range  dependence  is  in  good  agreement  with  the  ex¬ 
perimental  data  of  Sheehy,^  who  measured  fluctuations  in  a  24  kHz  pulsed 
signal  transmitted  over  ranges  of  approximately  30  to  3000  meters.  Further 

verification  of  this  theoretical  range  dependence  was  given  by  the  model 
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tank  experiments  of  Stone  and  Mintzer.  Chernov  reproduced  the  Bergmann 
and  Mintzer  results  and  in  addition  calculated  the  transverse  and  longitudi¬ 
nal  spatial  correlation  of  the  fluctuations.  Chernov's  analysis  of  fluctua¬ 
tions  and  their  correlation  was  extended  to  the  case  of  spherical  waves  by 
Karavainikov. 68 

All  of  the  theoretical  studies  mentioned  above  assumed  that  the 
refractive  index  field  in  the  acoustic  medium  could  be  described  by  an  ex¬ 
ponential  or  Gaussian  spatial  correlation  function.  Measurements  of  thermal 
microstructure  in  the  Pacific  Ocean  were  reported  by  Liebermann^  m  1951, 

— X 1 3i 

and  his  results  agreed  reasonably  well  with  the  exponential  function,  e  , 

where  a  is  the  mean  radius  of  the  thermal  patches.  (Liebermann  measured  a 
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mean  patch  size  of  0.6  meters.)  When  Whitmarsh  et  al.  made  simultaneous 
measurements  of  acoustic  fluctuations  and  thermal  raicrostructure  in  1957, 
however,  they  found  that  the  Kolmogorov  theory  of  turbulence  provided  a 
more  accurate  description  of  the  thermal  structure  of  the  sea  than  did  the 
exponential  or  Gaussian  functions.  Whitmarsh  et^  aJL.  suggest  that  the  agree¬ 
ment  of  the  Kolmogorov  theory  with  their  experimental  results  was  a  conse¬ 
quence  of  the  'freezing'  of  the  thermal  patches  after  their  turbulent  motion 
has  been  damped  out.  They  reason  as  follows: 


Since  the  static  heat  conductivity  is  almost  1000  times 
smaller  than  the  turbulent  conductivity,  the  patches  continue 
to  exist,  even  after  the  turbulent  motion  has  disappeared. 

This  stability  seems  to  lead  to  excellent  statistical  condi¬ 
tions  and  to  a  temperature  distribution  that  agrees  with  the 
experimental  results .... ^ 

This  'Kolmogorov'  model  of  the  medium  has  been  used  by  Tatarski^  and 
Medwin  ^  to  predict  sound  phase  and  amplitude  fluctuations  in  a  turbulent 
medium.  The  Kolmogorov  model  will  be  used  in  the  present  study,  so  a  brief 
description  of  turbulence  theory,  and  some  of  the  results  of  Tatarski  and 
Medwin,  will  be  presented  later  in  this  chapter. 

The  research  cited  above  all  deals  with  situations  in  which  the  total 


fluctuations  imposed  on  the  propagating  wave  are  sufficiently  small  that 
only  single  scattering  is  involved.  By  the  end  of  the  1960s,  work  had  be¬ 
gun  on  the  more  difficult  problem  of  multiple  scattering,  where  fluctuations 

in  the  wave  parameters  can  become  large.  Multiple  scattering  is  discussed 
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at  some  length  in  the  books  of  Tatarski,  Uscinski,  and  Ishimaru.  In 

this  thesis  we  will  be  concerned  only  with  single  scattering.  More  detailed 

surveys  of  theoretical  work  for  both  single  and  multiple  scattering,  in- 
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eluding  extensive  bibliographies,  may  be  found  in  Ishimaru  and  Barabanenkev 


3.2  Methods  of  Describing  the  Inhomogeneous  Medium 

A  theoretical  analysis  of  the  propagation  of  linear  acoustic  waves  in 
an  inhomogeneous  medium  requires  statistical  methods  for  describing  the 
medium.  The  medium  may  be  modelled  as  a  randomly  varying  field  of  refrac¬ 
tive  index,  or  equivalently,  of  sound  speed.  The  sound  speed  is  then  given 
by 


c  =  c (x,y,z,t) 


and  the  refractive  index  is 


f 
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n(x,y,z,t) 


c 

o 


c(x,y,z,t) 


1  +  P  (x,y,z,t) 


(3.1) 


where 

=  <c>  is  the  mean  sound  speed, 

<  >  denotes  ensemble  averaging,  and 

p(x,y,z,t)  is  the  deviation  from  unity  in  the  refractive  index. 

Variations  in  the1 ref ractive  index  are  usually  very  small;  for  example,  Urick 

and  Searfoss^  measured  a  value  of  <y^>  =  8  x  10  ^  at  a  depth  of  about  6 
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meters  in  a  mixed  layer  off  Key  West,  Florida.  Liebermann  observed  mean 

2  -9 

square  refractive  index  variations  of  <p  >  =  5  *  10  at  a  depth  of  50  meters 
off  the  coast  of  California.  Because  the  quantities  c  and  n  are  random 
variables,  a  description  of  their  spatial  dependence  requires  statistical 
parameters  such  as  the  correlation  function  and  its  Fourier  transform,  the 
spatial  wavenumber  spectrum.  Both  the  correlation  function  and  the  wave- 
number  spectrum  are  used  frequently  in  what  follows,  so  we  will  define  each 
of  these  parameters. 

If  p(x,y,z,t)  may  be  assumed  to  be  an  ergodic  process,  then  the 
fluctuations  may  be  characterized  by  the  correlation  function^ 


l 

I 


N^2  “  V (x^ ,  y^,  ^i >  b)  P (X2 ,  y2»  z2 ‘  ^ )  »  (3.2) 

where  the  overbar  denotes  time  averaging.  Because  of  the  ergodic  hypothesis, 
ensemble  averaging  and  time  averaging  are  equivalent.  If  it  is  further  re¬ 
quired  that  the  process  be  spatially  homogeneous,  then  the  correlation  func¬ 
tion  depends  only  upon  the  separations  Ax,  Ay,  and  Az,  where  Ax  ■  x9  -  x^. 

Ay  =  y2  “  yv  and  Az  =  z2  -  Zy  When  points  (x^  y^,  z^  and  (x2»  y2»  z?) 


are  the  same,  so  that  Ax  =  Ay  =  Az  =  0,  then  the  correlation  function  has  a 
~2 

maximum  value  of  y  .  As  the  separation  between  the  two  points  is  increased, 
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the  correlation  function  decreases  in  value  until,  at  infinite  separation, 
h'12  becomes  zero. 

If  the  variations  in  refractive  index  are  isotropic  as  well  as 
homogeneous,  then  the  correlation  function  will  depend  only  upon  the  mag¬ 
nitude  of  the  separation  between  points.  This  condition  may  be  written  as 

N12  ' 

where 

p  =  /(Ax)^  +  (Ay)2  +  (Az)2 

— 2 

It  is  useful  to  normalize  the  correlation  function  by  dividing  by  p  . 
This  results  in  the  correlation  coefficient  R,  given  by 

R  a  N12^y2 

Two  correlation  coefficients  frequently  used  to  describe  the  refractive 
index  variations  in  the  ocean  are  the  exponential  function, 

R(p)  =  exp(-p/a)  ,  (3.3) 


and  the  Gaussian  function, 


R(p)  =  exp(-p2/a2) 


(3.4) 


where  a  is  a  constant  corresponding  to  the  mean  patch  radius.  It  should 
be  noted  that  the  derivative  of  the  exponential  function,  which  may  be 
found  from  Eq.  (3.3)  to  bo 


dR 

dp 


-  exp(-p/a) 
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has  a  nonzero  value  at  p  =  0.  This  implies  a  discontinuity  in  the 
correlation  function  at  p  =  0,  which  is  possible  only  if  the  variation  in 
refractive  index,  p(p),  is  discontinuous,  a  condition  which  is  unrealistic 
in  the  ocean. ^  The  Gaussian  function,  Eq.  (3.4),  does  not  exhibit  this 
discontinuity,  and  in  this  respect  is  a  more  realistic  representation  of  the 
refractive  index  fluctuations  than  is  the  exponential  function. 

An  alternate  way  to  describe  spatial  variations  in  refractive  index  is 
to  use  a  spatial  wavenumber  spectrum.  Physically,  this  spectrum  character¬ 
izes  the  distribution  of  patch  sizes  that  are  present  in  an  inhomogeneous 
medium.  As  mentioned  in  Section  3.1,  the  distribution  of  thermal  patches 
in  the  ocean  is  a  result  of  turbulent  mixing.  Similarly,  the  thermal  micro¬ 
structure  in  a  model  tank  such  as  used  in  the  present  study  is  generated  by 
turbulent  mixing.  Before  returning  to  our  discussion  of  acoustic  wave 
propagation  in  an  inhomogeneous  medium,  we  will  consider  some  basic  concepts 
of  turbulence  theory. 

3.3  Turbulence 

The  theory  of  turbulence^ ’ ^ ^  used,  in  general,  to  describe 
viscous  fluid  flow  for  large  Reynolds  numbers.  The  Reynolds  number.  Re,  is 
defined  as 


Re 


where 

L  is  the  characteristic  scale  of  flow, 
v  is  the  characteristic  flow  velocity,  and 
v  is  the  kinematic  viscosity  of  the  fluid. 

For  small  values  of  Re,  the  fluid  flow  is  orderly  or  laminar.  When  the 
Reynolds  number  exceeds  a  critical  value,  Kecr>  the  flow  becomes  unstable 
and  breaks  up  into  turbulent  eddies.  Each  eddy  of  size  i  will  have 
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associated  with  it  a  local  Reynolds  number.  Re, .  If  Re.  also  exceeds  Re  , 

12.  cr 

the  eddy  will  further  break  down  into  smaller  eddies.  This  process  will 
continue  until  the  eddies  are  small  enough  that  viscous  dissipation  balances 
out  the  energy  being  supplied  to  the  eddies  from  the  external  source.  When 
this  occurs,  the  turbulence  reaches  steady  state,  and  a  range  of  eddy  sizes 
exists  such  that 


1  <  2  <  L 

o  o 


where 

L  is  the  outer  scale  of  turbulence,  and 
o 

1  is  the  inner  scale  of  turbulence, 
o 

The  outer  scale,  Lq,  is  determined  by  the  boundary  conditions  of  the  fluid 

flow.  In  an  acoustics  application,  Lq  may  generally  be  taken  to  be  the 

dimension  from  the  acoustic  source  or  receiver  to  the  nearest  boundary  of 
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the  medium.  The  inner  scale,  1  ,  is  shown  by  Tatarski  to  be 

o 

10  =  ,  (3.5) 


where 

e  is  the  energy  dissipated  as  heat  per  unit  mass  per  unit  time. 

This  qualitative  description  may  be  summarized  in  terms  of  energy  as  follows 
an  external  energy  source  supplies  energy  to  eddies  of  size  Lq  which  trans¬ 
fer  this  energy  down  a  chain  of  successively  smaller  eddies,  the  energy 
finally  being  dissipated  as  heat  due  to  the  viscous  losses  associated  with 
the  smallest  eddies. 

In  describing  turbulence  it  is  often  convenient  to  deal  with  the 
wavenumber,  x,  associated  with  an  eddy  rather  than  use  its  characteristic 
dimension,  l.  The  wavenumber  is  inversely  proportional  to  the  eddy  size 
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and  may  be  written  as 
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k 


2tt 

X 


(3.6) 


It  can  be  seen  that  large  wavenumbers  correspond  to  small  eddy  sizes,  and 
vice  versa ■  Because  the  kinetic  energy  of  the  flow  is  distributed  through¬ 
out  a  spectrum  of  eddy  sizes,  it  is  possible  to  define  a  power  spectral 

density  for  the  flow.  The  average  kinetic  energy  of  the  flow  is  then  the 
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sum  of  the  energy  associated  with  all  the  eddies;  i.e.. 


T=/V‘ 


)  dtc 


where 

T  is  the  average  kinetic  energy  per  unit  mass  per  unit  time, 

■J^Ck)  is  the  one  dimensional  power  spectral  density,  and 
the  flow  is  assumed  to  be  homogeneous  and  isotropic. 

The  quantity  <J>v(k)  is  related  to  the  velocity  fluctuations  in  the  flow, 
but  because  temperature  may  usually  be  treated  as  a  passive  additive  to  the 
turbulence,  the  spectra  of  the  temperature  and  refractive  index  variations 
are  usually  assumed  to  have  the  same  shape  as  $  (k).^ 

A  spectrum  proposed  by  Medwin^  for  describing  refractive  index 
variations  is  shown  in  Fig,  3.2.  In  the  figure,  the  range  of  wavenumbers 
corresponding  to  large,  anisotropic  eddies  is  labelled  the  'source'  sub¬ 
range,  as  these  eddies  supply  energy  to  the  entire  spectrum.  The  energy  of 
the  source  eddies  gradually  becomes  less  anistropic  (directional)  in  the 
transition  subrange,  until,  in  the  inertial  subrange,  the  energy  is 
isotropic  and  homogeneous.  In  the  inertial  subrange,  the  spectral  density 
is  given  by  the  simple  relation: ^ 


<J>  (k)  =  b< 
v 


•5/3 


(3.7) 


TYPICAL 
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where 

b  is  a  function  of  the  viscosity,  v,  and  the  energy  dissipation  rate, 


The  -5/3  power  law  expressed  by  Eq.  (3.7)  was  originally  proposed  by 
83 

Kolmogorov,  and  has  been  experimentally  verified  in  a  number  of 
70  71  78 

studies.  ’  ’  The  turbulent  energy  of  the  flow  is  finally  lost  to  heat 

in  the  'dissipation'  subrange,  where  viscous  forces  become  dominant. 

The  various  subranges  of  the  spectrum  in  Fig.  3.2  are  bounded  by  the 

wavenumbers  k  ,  k  ,  and  k  .  Medwin  defines  these  boundary  wavenumbers  as 
m  t  o 


follows 
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k  =  where  D  =  depth,  (3.8) 

ra  zu 

kq  =  (e/v^)^^  ,  and  (3.9) 

k  =  0. 5 (k  k  )1/2  .  (3.10) 

t  o  m 

The  wavenumber  k  .  defined  by  Eq.  (3.8),  depends  only  upon  the  depth 
m ' 

of  the  acoustic  experiment.  In  the  upper  ocean,  the  sea  surface  will  be 
the  nearest  physical  boundary  to  the  acoustic  experiment,  so  it  seems  rea¬ 
sonable  that  this  boundary  will  determine  the  maximum  significant  eddy  size. 
The  effect  on  the  acoustic  experiment  of  eddy  sizes  corresponding  to  wave- 
numbers  less  than  may  be  assumed  to  be  negligible.2*^ 

The  lower  limit  of  the  inertial  subrange,  is  an  empirical  value 

that  is  shown  by  Medwin  to  obey  the  relation  expressed  in  Eq.  (3.10).  The 
upper  limit  of  the  inertial  subrange  is  the  Kolmogorov  wavenumber,  kq,  given 
by  Eq.  (3.9).  Medwin  assumes  that  the  spectrum  can  be  truncated  at  kq,  and 
that  the  effects  of  the  dissipation  subrange  upon  an  ccoustic  experiment 
can  be  ignored.2^*  However,  the  results  of  Chotiros  and  Smith^^  suggest  that 


for  the  model  tank  used  in  the  present  investigation  it  is  more  reasonable 
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to  assume  that  the  truncation  wavenumber  of  Pao  is  valid.  Pao's  theory 
assumes  that  the  spectrum  of  the  thermal  raicrostructure  in  the  dissipation 
subrange  departs  from  that  of  turbulence  and  is  controlled  by  the  diffusiv- 
ity.  Pao  suggests  that  the  -5/3  power  law  continues  to  be  valid  in  the 
dissipation  subrange,  with  a  truncation  wavenumber  given  by 


K 


P 


<e/DV/4 


(3.11) 


where 

D  is  the  diffusivity. 

When  k  is  used  to  truncate  the  spectrum  rather  than  <  ,  the  wavenumber  k 
P  o’  t 

will  become: ^ 


k  =  0.5(k  k  )1^2 
t  p  m 


(3.12) 


In  this  section  and  the  previous  section,  two  alternative  methods  of 
describing  an  inhomogeneous  medium  have  been  presented:  the  method  of  the 
correlation  function  and  that  of  the  spectral  density  function.  The  two 
methods  are  not  independent,  however;  they  are  related  by  the  Fourier 
transform  theorem.  In  three  dimensions,  for  a  homogeneous  medium  this  re¬ 
lation  may  be  written  as: 


Ay,  Az) 


(3.13) 


X  exp{  j  (k^Ax  +  K2Ay  +  k^Az)  Jdk^di^dtc^ 


W  k2’  k3^  t*1e  t*iree  dimensional  spectrum  of  the  refractive 
index  variations,  and 


where 
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k^,  k 2 »  and  k3  are  wavenumbers  in  the  x,  y,  and  z  dimensions, 
respectively. 

For  an  isotropic  medium  Eq.  (3.13)  becomes ^ 

00 

N12(p)  =  f  ®(K)dK  ,  (3.14) 

J  Kp  P 

O 

where  the  one-dimensional  spectral  density  is  related  to  the  isotropic  three- 
dimensional  spectral  density  by 

2 

4>u(k)  =  4ttk  S  (k) 

3.4  The  Wave  Equation  for  an  Inhomogeneous  Medium 

In  the  preceding  sections  we  have  discussed  some  of  the  causes  of 

acoustic  signal  fluctuations  and  have  presented  two  methods  that  are  useful 

in  describing  the  inhomogeneous  medium  in  which  these  fluctuations  occur. 

Now  we  will  consider  in  more  detail  the  way  in  which  medium  inhomogeneities 

produce  variations  in  the  amplitude  and  phase  of  a  propagating  wave.  -The 

acoustic  wave  equation  for  an  inhomogeneous  medium  has  been  derived  by  a 
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number  of  people,  including  Rayleigh.  The  discussion  in  this  section  will 

87 

be  based  on  Chernov's  development. 

It  is  assumed  that  the  medium  is  lossless  and  is  in  a  state  of 
equilibrium  so  that,  ignoring  the  force  of  gravity,  the  ambient  pressure  po 
is  constant  throughout  the  medium.  The  only  inhomogeneities  present  in  the 
medium  are  those  in  temperature  and  density;  the  effects  of  bubbles,  bio¬ 
logical  matter,  etc.,  are  neglected.  Temperature  and  density  variations 
are  not  independent  quantities.  In  order  for  the  pressure  to  remain  con¬ 
stant,  the  changes  in  temperature  at  any  point  in  the  medium  will  be  accom¬ 
panied  by  a  related  change  in  density.  In  the  ocean,  it  would  be  expected 
that  the  temperature  inhomogeneities  would  be  moved  about  in  the  medium  by 
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current  flow  and  convective  motion.  It  is  assumed  that  these  types  ol 
motion  are  slow  compared  to  the  velocity  of  sound  so  that  the  inhomoge¬ 
neities  can  be  regarded  as  fixed  during  the  time  required  for  an  acoustic 
wavefront  to  propagate  through  a  patch  of  diameter  2a.  This  condition 
is  what  Chernov  calls  the  'quasi-static'  condition. 

The  wave  equation  for  the  situation  described  above  is 

2 

—2  "  y2P  +  vloS  PT  *  VP  =  0  *  (3.15) 

c  3t 


where 

c  =  c(x,y,z)  =  cq  +  Ac  is  the  sound  speed, 

c  =  <c>  is  the  mean  sound  speed, 
o 

Ac  =  Ac(x,y,z)  is  the  random  variation  in  sound  speed, 

P  is  the  acoustic  pressure, 

PT  =  PT(x,y,z)  =  pQ  +  Ap  is  the  total  density, 

PQ  =  <PT>  is  the  mean  density,  and 

Ap  =  Ap(x,y,z)  is  the  random  variation  in  density. 

This  equation  differs  from  the  homogeneous  wave  equation  in  two 

respects:  (1)  the  sound  velocity  is  not  constant,  but  varies  spatially  in 

a  random  manner;  and  (2)  there  is  an  additional  terra,  Vlog  p  .  Vp,  which 

arises  as  a  result  of  the  random  variations  in  density.  It  is  assumed  that 

the  variations  in  sound  speed  and  density  are  small;  i.e.,  Ac<<c  and  Ap<<p  . 

o  o 

When  some  'primary'  acoustic  wave  is  transmitted  through  the  inhomo¬ 
geneous  medium  described  above,  each  element  of  the  medium  becomes  a  source 
of  'secondary'  scattered  waves.  For  example,  suppose  a  harmonic  plane  wave 
propagates  in  the  x  direction.  This  primary  or  incident  wave  can  be  de¬ 
scribed  by 


PA  =  PA  exp[ j (wt  -  kx) ] 


(3.16) 
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where 

a)  is  the  angular  frequency,  and 

k  =  —  is  the  wave  number, 
c 

o 

If  the  scattering  elements  are  associated  with  the  coordinates  (£,n,C)  and 
the  observation  point  is  designated  by  (x,y,z),  then  the  distance  from  a 
scattering  element  to  the  observation  point  will  be 

r  =  /(x  -  O2  +  (y  -  n)2  +  (z  -  c)2  •  (3.17) 


The  total  scattered  pressure  at  the  observation  point  will  be  given  by  the 
solution  to  Eq.  (3.15),  which  Chernov  expresses  as: 


sc 


■  “  /  P 


Ac 

c 


_  ifc 

p_ 


3(AP1 
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1- 


exp[-jk(r  +  £)]  dv 


(3.18) 


where 

^"denotes  integration  over  the  volume  of  scatterers. 

V 

[In  Eq.  (3.18)  and  hereafter  the  factor  is  dropped.]  The  first  term  in 

curly  brackets  [Eq.  (3.18)]  is  associated  with  scattering  of  sound  by  spatial 
variations  in  the  sound  velocity,  while  variations  in  the  density  give  rise  to 
the  second  term.  In  the  ocean,  the  density  fluctuations  may  be  assumed  to  be 
negligible  in  comparison  to  the  sound  speed  fluctuations,  so  the  second  term 
in  curly  brackets  [Eq.  (3.18)]  can  be  omitted  in  approximating  the  scattered 
pressure. 

One  further  modification  of  Eq.  (3.18)  can  be  made  as  follows.  From 
Eq.  (3.1)  we  can  write  the  refractive  index  variations  as 


c 

o 


+  Ac) 


-  1 


(3.19) 


where  c  =  c  +  Ac.  Some  algebraic  manipulation  of  Eq.  (3.19)  gives 
o 


P 


Ac  .  Ac 


(c  +  Ac)  c 
o  o 


(3.20) 


where,  as  mentioned  above,  Ac  <<  c  .  Substitution  of  this  result  for  —  in 

o  c 

o 

Eq.  (3.18)  leads  to  the  following  expression  for  the  scattered  pressure: 

psc  "  ~2~  f  £gP,IzJj^r  *  n,  C)  dv  .  (3.21) 

The  physical  significance  of  this  result  can  be  summarized  as  follows:  A 
scattering  element  with  volume  dv  and  refractive  index  variation  p  is  lo¬ 
cated  at  point  (£,  n,  ?)•  Due  to  the  influence  of  an  incident  wave  p^,  the 
scattering  element  acts  as  a  secondary  sound  source  which  radiates  a  pres¬ 
sure  wave  dPsc-  The  total  scattered  pressure,  Psc>  received  at  the  obser¬ 
vation  point  (x,y,z)  is  the  sum  of  the  contributions  from  all  scattering 
elements  in  the  volume  V.  The  result  expressed  by  Eq.  (3.21)  can  be  used 
to  calculate  the  amplitude  and  phase  fluctuations  of  the  pressure  wave  at 
the  observer,  as  will  now  be  discussed. 


3.5  Fluctuations  in  the  Amplitude  and  Phase  of  the  Observed  Pressure 
The  total  pressure  at  the  observation  point  will  be  the  sum  of  the 
primary  pressure  and  the  scattered  pressure, 


P  e 


-i* 


Pi  + 


where 

P  and  <j>  are  the  amplitude  and  phase  of  the  total  pressure,  p. 


As  the  inhomogeneities  move  about  in  the  medium,  the  amplitude  and  phase  of 
the  observed  pressure  will  vary  randomly  about  their  mean  values.  It  is 


47 


convenient  to  define  the  fluctuation  in  observed  pressure  amplitude  and 
phase  about  their  mean  values  as 


B  m 


,  and 


(3.22) 


S  ■  (p  —  <  <p  > 


(3.23) 


where 

B  is  the  logarithmic  fractional  variation  of  the  instantaneous  pressure 
amplitude,  P,  from  its  mean  value,  <P>,  and 
S  is  the  variation  of  the  instantaneous  phase,  <{>,  from  its  mean  value, 
<$>. 

The  &n  operator  in  Eq.  (3.22)  is  a  consequence  of  using  Rytov's  method,  also 

known  as  the  method  of  smooth  perturbations,  in  calculating  the  fluctua- 
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tions.  In  this  thesis  we  will  be  concerned  with  amplitude  fluctuations 
sufficiently  small  to  make  the  approximation 


The  average  values  of  B  and  S  are  zero,  so  it  is  necessary  to  use  the  mean 

2  2 

square  fluctuations,  <B  >  and  <S  >,  to  characterize  the  level  of  fluctua¬ 
tions  caused  by  the  inhomogeneities.  There  are  a  variety  of  analytical  ex¬ 
pressions  in  the  literature  that  relate  the  mean  square  fluctuations  to  the 

statistical  parameters  used  to  describe  the  inhomogeneous  medium.  In  the 

2  2 

next  two  sections  we  consider  expressions  for  <B  >  and  <S  >  when  the  medium 
is  described  by  a  Gaussian  correlation  function,  and  when  the  medium  is  de¬ 
scribed  by  a  spectral  density  function. 
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3.5.1  Fluctuations  in  a  Medium  Described  by  a  Correlation  Function 

In  the  case  of  plane  waves  propagating  through  a  statistically  isotropic 

89  2  2 

Inhomogeneous  medium,  Chernov  obtained  solutions  for  <B  >  and  <S  >  in 

terms  of  a  wave  parameter,  D,  given  by 


where 

L  is  the  propagation  distance  of  the  acoustic  wave  through  the  medium. 
Chernov's  analysis  assumes  that  the  inhomogeneities  are  large  compared  to 
an  acoustic  wavelength  (ka  >>  1),  and  that  the  propagation  distance  is  large 
compared  to  the  scale  of  the  inhomogeneities  (L  >>  a). 

For  large  values  of  the  wave  parameter,  which  corresponds  to  the 
wave  region  of  the  patches,  Chernov  obtains: 

00 

<B2>  =  <s2>  =  <P2>  k2L J R(p')  dp  ,  D>>1  .  (3.24) 

o 

If  the  correlation  coefficient  is  Gaussian,  given  by  Eq.  (3.4),  then  Eq. 
(3.24)  becomes 


<B2>  =  <SZ>  =  <p2>  k2aL  ,  D»1  .  (3.25) 

Mintzer  has  shown^’^  that  these  same  results,  Eqs.  (3.24)  and  (3.25),  are 
valid  for  a  spherical  wave  in  the  wave  region  when 

k2<p2>aL  «  1 

For  small  values  of  the  wave  parameter,  corresponding  to  the 
ray  region  of  the  patches,  Chernov  shows  the  mean  square  amplitude 
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fluctuations  for  a  plane  wave  to  be 


2  12  3 

<B  >  =  ^  <y  >  LJ 


/ 


V2V2R(p)dp 


D<<1 


(3.26) 


which,  for  a  Gaussian  correlation  coefficient,  becomes 


<B2>  =  <y2>  ^3  ,  D<<1  .  (3.27) 

a 

89 

The  phase  fluctuations  for  D<<1  have  been  shown  to  be  twice  their 
values  when  D>>1  (i.e.,  double  the  results  in  Eqs.  (3.24)  and  (3.25)). 

Finally,  when  the  wave  parameter  is  of  the  order  of  unity,  Chernov  obtains 


<B  >  -  2  <P2>  k2a  L^l  -  i  arctan  ,  (3.28) 

<s  >  =  —  <y2>  k2a  L^l  +  i  arctan  ,  (3.29) 

where  a  Gaussian  correlation  coefficient  has  been  assumed.  For  small  and 
large  values  of  D,  these  last  two  results  become  equivalent  to  the  ray.  and 
wave  region  solutions,  respectively. 
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Using  the  same  conditions  and  assumptions  as  Chernov,  Karavainikov 
obtained  a  similar  set  of  results  for  the  amplitude  and  phase  fluctuations 
in  a  propagating  spherical  wave.  For  the  wave  region  of  the  inhomogeneities, 
D  »  1,  Karavainikov' s  result  is  identical  to  Eq.  (3.25). 


3.5.2  Fluctuations  in  a  Turbulent  Medium 

In  the  case  of  a  medium  whose  random  refractive  index  field  is 

19  71 

determined  by  turbulence,  Tatarski  ’  has  found  expressions  for  amplitude 
fluctuations  in  terms  of  a  turbulence  parameter,  C^.  For  propagation  dis¬ 
tances  L  such  that  L  >>  /\L  >>  1  ,  where  1  and  L  are  the  inner  and  outer 

o  ooo 

scales  of  turbulence,  respectively,  and  \  is  the  acoustic  wavelength, 


Tatarski  shows 


<B2>  =  0.13  C  2  k7/6  L11/6 
n 


(3.30) 


for  spherical  waves,  and 


<B2>  -  0.31  C  2  k7/6  L1176 
n 


(3.31) 


for  plane  waves.  The  quantity  is  related  to  the  spectral  function  S^(k) 


by 
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S  (*0  0.033  2  K  K  <  K  <  IC 

y  n  t  p 

An  assumption  made  in  deriving  Eqs.  (3.30)  and  (3.31)  is  that  only  the 
inertial  range  of  the  spectrum  contributes  to  the  acoustic  fluctuations. 
If  the  spectrum  is  of  the  form  shown  in  Fig.  3.2,  so  that  there  is  an  ap¬ 
proximately  flat  transition  region,  then  there  is  an  additional  contribu¬ 
tion  from  the  transition  region  to  the  amplitude  fluctuations  of  an 
84 

amount 


_  2  .  it  3.  4  4. 

<B  >  -  ——  $  L  (k  -  k  ) 
t  480  m  t  m 


(3.32) 


for  spherical  waves,  and 


> 


.3,  * 

48  'm  L  (Kt 


•  *  * 

=  —  S> 


(3.33) 


for  plane  waves.  As  shown  in  Fig.  3.2,  $  is  the  maximum  value  of  the 

spectral  density  function  $  .  The  total  mean  square  amplitude  fluctuations, 

2  2 

when  there  is  a  contribution  due  to  the  transition  region,  is  <B  >  +  <B  >. 
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Often  the  amplitude  fluctuations  are  expressed  in  terms  of  the 
coefficient  of  amplitude  variation,  CAV,  because  this  parameter  is  easy  to 
calculate  from  a  record  of  pulsed  data.  The  coefficient  of  amplitude  varia¬ 
tion,  which  is  essentially  the  standard  deviation  of  the  pressure  amplitude, 
is  defined  as^2 


CAV 


2 


<P2> 


-  <P> 


:P> 


(3.34) 


where  P  is  the  amplitude  of  the  observed  pressure.  For  small  amplitude 
2  2 

fluctuations,  <B  >  and  CAV  are  approximately  equal. 

3.6  Correlation  of  Fluctuations 

In  addition  to  the  mean-square  fluctuations,  it  is  useful  to  know  the 
spatial  correlation  of  the  fluctuations.  As  will  be  seen  in  Chapter  4, 
these  correlations  are  important  quantities  in  developing  a  theory  of  para¬ 
metric  reception  in  a  random  medium.  In  this  section,  we  present  some  of 
the  correlation  functions  that  will  be  usefu]  in  this  study.  There  are  two 
types  of  spatial  correlations  that  are  of  interest:  longitudinal  correla¬ 
tion,  which  is  the  correlation  of  fluctuations  between  two  point  receivers, 
separated  a  distance  AL,  and  both  located  on  the  axis  of  the  propagating 
wave;  and  transverse  correlation,  which  is  the  correlation  of  fluctuations 
between  two  point  receivers  separated  a  distance  Ap  in  a  direction  trans¬ 
verse  to  the  acoustic  axis. 

For  a  plane  wave  in  a  Gaussian  medium,  the  longitudinal  correlation 
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coefficients  for  the  amplitude  and  phase  fluctuations  are  given  by 


,  D  >>  1 


for  the  wave  region  of  the  inhomogeneities,  and  by 
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(3.34) 


rl  1 


. . 

©TV  g'1" . . 

in  the  ray  region.  In  these  expressions  refers  to  the  propagation  dis¬ 
tance  to  t, e  receiver  nearest  the  signal  source.  Similar  expressions  for 

the  transverse  correlation  coefficients  for  a  plane  wave  in  a  Gaussian 

. .  92 

medium  are: 


R*  =  r!,  =  exp[- (Ap)2/a2]  ,  !)  >>  1  '  (3.37) 

I)  b 

Kj  ^  exp[-(Ap)2/a2]  •  jl  -  2[  (A,.)2/a2  J  +  ~[  (:\n)2/a2]2J  ,  D  «  1 

*  »  (3.38) 

R^  =  exp[-(Ap)2/a2] ,  0  <<  1  .  (3.39) 

There  are  also  published  results  describing  the  correlation  of 
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fluctuations  for  spherical  waves  and  for  turbulent  media.  ’  These  re¬ 
sults  in  general  involve  complicated  expressions  from  which  it  is  difficult 
to  gain  a  physical  understanding,  and  we  will  not  reproduce  them  here. 

3.7  Summary 

In  this  chapter,  a  review  of  the  basic  theory  of  acoustic  wave 
propagation  in  an  inhomogeneous  medium  has  been  presented.  The  material 
has  been  selected  on  the  basis  of  its  relevance  to  the  theoretical  analysis 
of  parametrically  received  signal  fluctuations  to  be  discussed  in  the  fol¬ 
lowing  chapter. 

It  was  seen  in  Section  3.1  that  the  principal  cause  of  scattering  in 


the  ocean  is  the  thermal  microstructure.  Thermal  patches  can  be  treated  as 
sound  sources  of  radius  a.  The  scattered  radiation  from  these  patches  has 
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a  nearfield  extending  to  a  distance  ka  ;  from  that  distance  onward  the 
radiation  is  spherically  spreading  with  a  beamwidth  of  1/ka  radians.  At 
any  point  in  the  medium  the  observed  pressure  is  the  sum  of  the  unscattered 
pressure  p^  and  the  scattered  pressure  Psc*  Due  to  random  changes  in  the 
characteristics  of  the  scattered  pressure  field,  the  observed  pressure  is  a 
randomly  fluctuating  quantity. 

In  Section  3.2  two  methods  were  discussed  that  can  be  used  to  describe 
an  inhomogeneous  medium.  One  method  uses  a  correlation  function  for  the 
variations  of  refractive  index,  and  the  other  method  uses  a  spatial  wave- 
number  spectrum.  Because  the  spectral  composition  of  the  inhomogeneities 
is  determined  by  turbulence,  Section  3.3  gave  a  brief  summary  of  some  con¬ 
cepts  from  the  theory  of  turbulence. 

The  wave  equation  for  an  inhomogeneous  medium  and  its  integral  solution 
were  discussed  in  Section  3. A.  A  summary  of  expressions  for  the  amplitude 
and  phase  fluctuations  of  the  observed  pressure,  as  well  as  for  the  corre¬ 
lation  of  these  fluctuations,  was  presented  in  the  remainder  of  the  chapter. 

The  background  material  discussed  in  Chapters  2  and  3  will  now  be.  used 
to  develop  a  theory  of  parametric  reception  in  an  inhomogeneous  medium. 


CHAPTER  4. 


THEORETICAL  ANALYSIS 


4.1  Introduction 

In  Chapter  2,  expressions  were  given  for  the  second-order  pressure 
associated  with  both  nearfield  and  farfield  parametric  receivers.  The 
results  assumed  the  parametric  receiver  to  be  operating  in  a  medium 
containing  no  inhomogeneities.  In  this  chapter,  some  effects  of  medium 
inhomogeneities  will  be  taken  into  account,  and  an  expression  will  be 
derived  that  describes  amplitude  fluctuations  in  the  second-order 
pressure  wave  of  a  parametric  receiver. 

Before  beginning  a  formal  analysis,  it  may  be  helpful  to  consider 
some  general  aspects  of  the  problem.  There  are  three  separate  signals 
which  will  be  affected  by  the  inhomogeneities.  One  is  the  low-frequency 
signal  that  is  to  be  detected  by  the  receiver.  In  practice  this  signal 
may  be  expected  to  propagate  a  significant  distance  through  the  ocean, 
and  as  a  result  of  scattering,  it  will  arrive  at  the  receiver  as  a 
randomly  varying  signal.  The  pump  wave  will  also  be  scattered  as  it 
travels  through  the  interaction  region,  so  nonlinear  interaction  will 
be  occurring  between  two  fluctuating  first-order  signals.  The  resulting 
second-order  signal  can  be  expected  to  vary  in  a  manner  that  is  related 
to  the  first-order  fluctuations.  Also,  the  second-order  signal  will 
be  scattered  as  it  propagates  from  each  source  point  to  the  observer, 
and  further  fluctuations  will  result. 

Analyses  have  been  reported  by  Smith‘d’ ^  and  by  Chotiros  and 
Smith^’^  describing  the  performance  of  a  parametric  transmitting  array 
in  a  random  medium  but,  to  date,  there  is  no  similar  study  for  receiving 
arrays.  In  this  chapter,  the  methods  developed  by  Smith  and  Chotiros 
are  used  to  study  the  effects  of  medium  inhomogeneities  on  parametric 
reception.  It  will  be  assumed  in  this  analysis  that  the  signal  source 
is  located  in  the  main  beam  of  the  parametric  receiver.  Expressions 
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predicting  the  amplitude  fluctuations  for  the  second-order  pressure  will 
be  developed  for  both  collimated  and  spherically  spreading  pump  waves. 


4.2  Amplitude  Fluctuations  for  a  Nearfield  Receiving  Array 
We  will  begin  our  analysis  with  the  simple  case  of  the  nearfield 
parametric  receiver  shown  in  Fig.  4.1.  A  signal  wave  originates  at 
a  source  located  on  the  acoustic  axis  of  the  parametric  receiver  at 
z  =  -Z  ;  the  source  is  sufficiently  far  removed  from  the  pump  transducer 
that  its  radiation  may  be  treated  as  being  planar  in  the  vicinity  of  the 
parametric  receiver.  In  rectangular  coordinates,  the  signal  wave  may 


be  represented  by 

ps  (z,t)  =  Ps[l  +  Bs(x,y ,z,t)]e  “sZ 

x  exp{  j  [oj  t  -  k  z  -4  -  S  (x,y  ,z  ,t)  ]  }  ,  (4.1) 

s  sos 

where 


P 

B 


s 

s 


=  Pg(t)  is  the  mean  pressure  amplitude  at  the  pump  transducer, 
P  (t)  -  P 

=  _s _ s  is  the  amplitude  fluctuation, 

P 


s 

S  =  Aft)  -  k  z  is  the  phase  fluctuation, 
s  s 

a  is  the  attenuation  coefficient  at  the  signal  frequency, 

3 

cb  is  the  phase  of  the  wave  at  z  =  0,  which  may  be 
o 

set  to  zero,  and 

the  overbar  indicates  time  averaging. 

Now  assume  that  the  hydrophone  is  situated  in  the  nearfield  of 


the  pump  transducer;  i.e.,  that 

L  <  4b2/Ap  » 

where 

z  =  L  is  the  location  of  the  hydrophone, 

2b  x  2b  are  the  dimensions  of  the  pump  transducer,  and 


X  is  the  wavelength  of  the  pump  wave. 

P 

For  this  case  the  pump  radiation  in  the  interaction  volume  of  the 
parametric  receiver  may  be  assumed  to  be  approximately  planar  and  may 


be  written  as 


GEOMETRY  FOR  ANALYZING  NEARFIELD  RECEIVER 
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time  derivatives  to  be  neglected  when  calculating  the  source  density. 
Such  an  approximation  is  reasonable  in  situations  where  the  signal 
frequencies  involved  are  much  higher  than  any  frequency  component 
in  the  fluctuation  spectra. 

The  interaction  frequency  pressure  produced  at  the  observation 
point,  r,  due  to  radiation  of  the  virtual  source,  qdv,  will  be  given  by 

dp+  (r , t)  =  -  Po  (1  +  B+)e  dS±  e  a±r 
4nr' 


x  q(ro>T)J  dv> 


(4.4) 


where 


B+  and  S+  are  the  amplitude  and  phase  fluctuations  produced  in 

the  interaction  frequency  wave  as  it  propagates  from 

source  point  to  observation  point,  and 
r ' 

T  =  t  -  is  the  retarded  time. 

c 

o 

An  attenuation  term,  e  a±r  ,  has  been  included  ad  hoc  in  Eq.  (4,4) 
to  account  for  attenuation  of  the  interaction  frequency  wave.  Substi¬ 
tution  of  Eq.  (4.3)  for  the  source  density  in  Eq.  (4.4),  and  integration 
over  the  interaction  volume,  V,  gives  the  following  expression  for  the 
second-order  pressure  at  the  observer: 


4  V‘ 

P±(r,t)  =  -  2~ 


P  c 
0  o 


f. 


X  (1  +  Bs)  (1  +  Bp) (1  +  B±)  e  j(Ss  ±  sp  +  S±) 


-ik+z  -ik+r'  1 

x  e  J  -  e  J  - r  dv. 

4irr 


(4.5) 


This  expression  for  p+  can  be  simplified  in  the  following  way. 
First,  assume  that  the  amplitude  and  phase  fluctuations  are  small 


B  «  1,  S  «  1, 
s  s 


compared  to  unity;  i.e. 
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B  «  1,  S  «  1, 

P  P 

B+  «  1,  S+  «  1. 

Then  the  amplitude  fluctuation  terms  can  be  expanded  and,  retaining 

only  terms  of  first  order  in  fluctuations, 

(1  +  B  )(1  +  B  )(1  +  B.)  =  1  +  B,  (4.6) 

s  p  ± 


where 


B  =  B  +  B  +  B , . 
s  p  ± 


Similarly,  the  phase  fluctuation  term  can  be  approximated  as 

e-j(Ss  ±  Sp  +  S+)  •  J  _  js 


(4.7) 


where 


s  =  s  ±  s  +  s  . 

s  p  ± 

Next,  assume  that  the  signal  frequency  is  much  lower  than  the  pump 
frequency,  and  that  the  pump  and  interaction  frequencies  are  approximately 
equal; 


w  «  to  ,  and 
s  p’ 

• 

U)  = 

p  ± 

It  may  then  be  expected  that  the  absorption  coefficients  corresponding" 
to  these  frequencies  will  obey  the  relations 

(4.8) 

(4.9) 


a  «  a  ,  and 
s  p 


a  =  0L  . 


Furthermore,  if  the  transverse  dimensions  of  the  interaction  region 
remain  small  compared  to  the  longitudinal  dimension,  then 


r '  =  (L  -  z) , 

so  that 

e"«±r’  •  e-a±(L  -  z)<  (4.10) 

Using  these  approximations  [Eqs.  (4.8)  -  (4.10)],  the  attenuation 
terms  in  Eq.  (4.5)  can  be  written  as 
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a-(as  +  ap)z  e-a±r'  ^  e-asz  (rxp  -  a+)z  e 


-ci  i  L 


-a+L 
=  e  -  . 


(A. 11) 


Now  Eqs.  (4.6),  (4.7),  and  (4.11)  can  be  used  to  rewrite  the 
expression  for  the  second-order  pressure  [Eq.  (4.5)]  as 


P±(r,t)  =  A  f(l  +  B)(l  -  jS)  ‘  '  tt  dv, 


„-jk±(z  +  r')  J_ 

r' 


(4.12) 


where,  for  convenience. 


SPgPpkt  e"0t+L 


^oco 


In  addition  to  Eq.  (4.12),  we  can  also  write  p+  explicitly  in  terms 

of  the  fluctuations  of  the  second-order  pressure.  Thus,  if  BpR  and 

S  are  the  amplitude  and  phase  fluctuations  of  the  interaction  frequency 
PR 

pressure  at  the  observation  point,  then  p+  can  be  written  as 

,"jsPR  (4 . 13) 


p+(r,t)  =  pH(l  +  BpR)  e 


where 


/*-jk+(z  +  r')  i 

ph  “‘Jo  '  -  dv 

is  the  second-order  pressure  that  would  be  produced  at  the  observation 
point  if  the  medium  were  homogeneous.  For  phase  fluctuations  small 
compared  to  unity,  the  exponential  term  in  Eq.  (4.13)  can  be  approxi¬ 
mated  by 

e~ j  SPR  i  i  _  1S  , 
e  J-  JaPR» 

so  that  the  interaction  frequency  pressure  becomes 

4P„  o  +  bpr-  JSPR),  (4-l«) 

where  only  terms  to  first-order  in  fluctuations  are  retained. 

It  is  now  possible  to  solve  for  the  amplitude  and  phase  fluctuations, 

B  and  S„„,  as  follows.  First  we  equate  Eqs.  (4.12)  and  (4.14) 

PR  PR 


PH(1  +  BpR  -  jSpR)  =  +  B  -  jS)  e"jk±(Z  +  r  5  r-  dv  ,  (4.15) 

V 
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where  cmiy  terras  to  first  order  in  fluctuations  have  been  retained  from 
Eq.  (4.12).  Next  the  homogeneous  component  of  the  pressure,  p  ,  is 


subtracted  from  both  sides  of  Eq.  (4.15),  leaving 
B 

where 


PR 


JSPR  -U(B  -  JS> 
R- 


-jk± (z  +  r’)  1  , 

e  7T  dv. 


(4.16) 


dv . 


The  amplitude  and  phase  fluctuations  of  the  interaction  frequency 
wave  at  the  hydrophone  may  now  be  found  by  equating  the  real  and  imaginary 
parts  of  Eq.  (4.16);  i.e.. 


B 


PR 


=  Re{  ±  /  (B  -  jS)  e~jk±(z  +  r,)  ±  dv  }  , 


and 


SpR  -  -Im  { 


\fa  - 1 

r 


(4.17) 


jS)  e-^<2  +  r,)  i.dv)  . 


In  this  analysis  we  will  be  concerned  only  with  the  amplitude 
fluctuations,  given  by  Eq.  (4.17).  Since  B  is  a  random  variable,  it  will 

r  K 

2 

be  useful  to  find  its  mean-square  value,  <Bpp  by  assuming  that  Bpp 
is  ergodic,  it  will  be  possible  to  take  the  (spatial)  ensemble  average 
rather  than  a  time  average.  Thus,  we  form  the  average  of  the  product, 

BPR1BPR2, 


where 

Re  { 

1 

fa  i 

-  jS1) 

-jk+(z. 

'PRl 

H 

e 

•/ 

fv 

‘PR2 

Re  { 

1 

H 

faz 

-  js2) 

-jk+(z 

e 

f  ^  v« 
r2  2 


and  the  geometry  of  Fig.  4.2  applies.  The  mean-square  amplitude 


fluctuation,  <B  >,  is  given  by* 
r  K 


*In  forming  the  product  BpR^BpR2  the  following  relation 
between  complex  numbers  a  and  b  is  used: 

Re  {  a  }  •  Re  {  b  }  =  -^  Re  {  a(b  +  b*)}, 

where  the  asterisk  denotes  the  complex  conjugate. 


(,  i 


<BPR  >  "BPR1BPR2' 


7  -  jSL>(H2  -  js,) 


x  e 


i k « (^i  +  7-p  +  ri  +  r2}  l  .  ^ 

7^  dvidv2 


f<(Ql  "  JSlUBi  +  JS2)>  ° 


-jk±(zL  "  z2  +  r[  "  rp 


x  Trirr  dv  dv ,} 

L  ^ L  9  -*•  “ 


(4.18) 


This  result  is  not  easy  to  use  as  it  will  require  numerical  integra¬ 
tion  for  its  evaluation.  However,  Eq .  (4.18)  can  be  simplified  considerably 
if  it  is  assumed  that  the  amplitude  and  phase  fluctuations  have  complete 
transverse  correlation  in  the  interaction  volume.  With  this  assumption, 

S^,  and  become  independent  of  the  variables  x  and  y,  and 


the  expressions  for  <B  >  become 

r  K 


<B  2> 
PR 


~jk±(zl  +  Z2). 


-^Re^yj^B^  -  jB1S2  -  jB2Sx  -  S1S2)>  e  "  x  “I,I„dz,dz, 


12  12 


+  i\s2  -  jB2s1  +  S1S2)>  e 


jk±(zi  ~  z2)l,l*,dz1dz„}  , 


0 

b  b 


12  12 

(4.19) 


where  I  =  [  (L  -  z) 2  +  x“  +  y2]‘}[(L  z)“  +  x  +  y  ]  dxdy 

-b  -b 

An  approximate  expression  for  I,  valid  in  the  nearfield  of  the  pump 
transducer,  is  shown  in  Appendix  1  to  be 


1  =  -  e 
k 


-jk(  (L  -  z) 


(4.20) 


A  similar  result  is  obtained  in  Appendix  1  for  H;  namely, 

-jkl. 


H 


IttL 


(4.21) 
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Substitution  of  Eqs.  (4.20)  md  (4.211  into  h. 


(4.19)  leads  to  Llm  result 


2  1  1 
<BDn  >  =  4  Re{ 


PR 


f  Abib->  "  iBis9  -  JV*.  -  s,s  )>  dz 

0  **  “  “ 

-jffa B1B2  +  J»1S2  -  1B2S1  +  S1 

’  I  Re  {F?/ _/2  '  »i»2  2  -) 


9  1  9 

L  0 


l2"o 


S2)>  dzJdz2  } 


12  <  B2S1  >)  dZldz2  }  . 


Because  <B^B2>  and  <B2;’p>  are  real  numbers,  this  expression  reduces  to 


<BPR  *  l2^AB1B2>  dzldz2 “ 


(4.22) 


This  result  for  the  mean-square  amplitude  fluctuations  of  the 

second-order  pressure  wave  can  be  evaluated  by  developing  an  expression 

for  the  correlation  term,  <R^B2  ’  as  a  function  of  z^  and  z2>  and 

performing  the  integration.  We  will  return  to  the  matter  of  evaluating 

and  interpreting  Eq.  (4.22)  below,  but  first  a  similar  expression  for 
2 

<BPR  ?  wili  be  develoPed  assuming  that  the  pump  wave  is  spherically 
spreading  rather  than  collimated  and  planar. 

4,3  Amplitude  Fluctuations  for  a  Farfield  Receiving  Array 
An  analysis  similar  to  that  in  the  previous  section  can  be 
developed  for  a  parametric  receiver  with  spherically  spreading  pump 
waves.  The  geometry  used  is  shown  in  Fig.  4.3.  As  for  the  nearfield 
analysis,  the  signal  source  is  located  on  the  axis  of  the  parametric 
receiver  and  is  sufficiently  far  from  the  pump  that  its  radiation  may 
be  assumed  to  be  planar  in  the  vicinity  of  the  parametric  receiver. 

The  signal  wave  may  be  represented  as 


Ps  =  V1  +  Vz  +  7-s)l  G  aSZ  GXP{  i[wst  -  kgz  -  Ss(z  +  Zg)]  }  , 

and  the  pump  wave  is 

Pp  =  j  (Pp/r)D(y)  [1  +  (r ) ]  e  ‘Pr  exp  |j[o)^t  -  k^r  -  (r)  ]  j  , 


•<  <  U  co 


b(S 


where 

P1  is  the  mean  pump  wave  pressure  amplitude  at  y  -  0,  i  -  1  m,  and 
P 

D(y)  is  the  directivity  function  for  the  pump  transducer. 

In  writing  the  signal  wave  fluctuation  terms,  and  Sg,  as  independent 

of  y,  complete  transverse  correlation  of  Bg  and  in  the  interaction 

volume  has  been  assumed.  Similarly,  it  is  assumed  that  the  pump  wave 

fluctuation  terms,  B  and  S  ,  are  completely  correlated  along  the 

P  P 

spherical  wavefronts  within  the  pump  beam  so  that  B^  and  are 
independent  of  y.  These  assumptions  are  discussed  in  more  detail  in 
Appendix  2. 

The  source  density  function  at  point  (r,y)  can  be  found  in  the  same 
way  as  in  the  nearfield  case  [F.q.  (A. 3)];  the  result  is 


-gP  P'  u.  -a  z  -a  r 

q  =  — |-V  c  8  c  *’  D(Y) 

per 
o  o 

-j(S  ±  Ss)  -j(k  r  ±  ksz)  ju)+t 

x  (1  +  B  ) (1  +  B  )e  P  e 

v  s  p 


e  '  .(4.23) 
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This  expression  can  be  used  to  find  the  second-order  pressure  at  the 
observer  in  the  same  way  as  in  the  previous  section.  In  the  case  of 
spherically  spreading  pump  waves,  however,  a  solution  can  be  obtained 
more  simply  by  adopting  a  procedure  developed  by  Berktay  and  Shooter. 

They  assume  that  the  sphericity  of  the  pump  wave  within  the  beam  is  small 
compared  to  the  wavelength  at  the  signal  frequency,  so  that,  for  the 
signal  source  located  on  the  Z  axis,  the  source  functions  q+  can  be 
assumed  to  be  cophasal  on  the  spherical  wavefronts.  The  frequency  of  the 
second-order  radiation  is  nearly  equal  to  the  pump  frequency,  so  the 
second-order  waves  will  radiate  spherically,  with  the  same  beam  pattern  as 
for  the  pump  wave.  These  assumptions  are  used  to  calculate  the  second-order 
pressure  at  the  observer  as  follows. 


The  elemental  particle  velocity  at  r  due  to  a  spherical  shell  of 


sources  of  thickness  5r  is 

S”r  =  o  q(r)  5r, 

where  the  source  density  function  is  given  by  Eq .  (4.23).  The  contri¬ 
bution  of  these  sources  to  the  particle  velocity  at  (L,0)  will  be 

-js 

=  (r/I.)  £ur(l  +  Bu)  e  "  e,p|  -  (V  +  jkJ(L  -  r)], 

where 

B  =  — ' — ■ -  — ^  is  the  amplitude  fluctuation, 

u  U(t) 


S  =  (t)  -  k  r  is  the  phase  fluctuation,  and 

U(t)  and  C t )  are  the  amplitude  and  phase  of  the  particle 

velocity,  respectively. 

The  total  particle  velocity  at  (L,0)  is  the  sum  of  contributions  from 
all  sources  in  the  interaction  region,  i.e., 


U±(L,0)  =  (2L)  exp  [  -  (a(  +  jk 

-js 


JL]f 

Jo 


(1  +  B  ) 


q(r)  r  exp  [  (Ct+  +  jk4)r]  dr. 


(4.24) 


The  second-order  pressure  at  the  point  (1,0)  can  now  be  found  by 
using  the  far-field  relation,  p  =  Pqcu,  in  connection  with  Eq.  (4.24). 
The  result  is 


-BPsP'us+ 

p  (L,0)  =  - — exp[-(a+  +  jk+)L] 

2p  c  L  -  - 

o  o 


■/ 

r\ 


(1  +  B  )  (1  +  B  )  (1  +  Bje 
s  P 


-J(s„ 


s  +  SJ 


x  exp[-(a  +  a  -  a+)r]dr,  (4.25) 

s  P  — 

where  B+  and  S+  are  amplitude  and  phase  fluctuations,  respectively, 
for  the  interaction  frequency  pressure  wave.* 


*It  may  be  shown  that  the  amplitude  fluctuations  in  the  pressure 
and  particle  velocity  are  equal:  i.e.,  B  =  B+ .  Similarly, 

SA.  " 


S  = 
u 
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For  small  fluctuations  this  result  becomes 

P+(L,0)  -  A  /(I  +  B)  (1  -  j  S)dr , 


(4.26) 


where 


"0 

-6F  P'w.  -(a  +  jk  )  L 

.  s  p  ±  -  - 

A  = - L—  e 

s  3. 

2.  c  L 
o  o 

B  =  B  +  B  +  B,,  and 
s  p  i 

S  =  S  ±  S  +  Sx. 
s  p  ± 

In  obtaining  Eq.(4.26),  the  attenuation  terms  have  been  approximated 
by  Eq. (4.11). 

As  in  the  case  of  the  nearfield  parametric  receiver,  the  second- 
order  pressure  can  be  written  in  terms  of  a  homogeneous  and  fluctuating 
component , 

(4.27) 

-  n  r  i\ 

where 


P±  =  PH(1  +  V  C 

rL 

/  dr  =  A  L 

*Jr 


PH  =  \ 

Eqs.(4.26)  and  (4.27)  may  be  equated,  and  if  the  fluctuations  are 
small  compared  to  unity,  we  obtain 


PH(1  +  SPR  jSPR^ 


*/0 


+  B  -  jS)dr. 


(4.28) 


The  homogeneous  component  is  subtracted  from  both  sides  of  Eq.(4.28), 
and  the  real  and  imaginary  components  of  the  result  are  equated  to 
give 


B 


PR 


*  L  fMV' 
J0 

r 

/  Sdr . 

Jo 


and 


PR  hJQ 

2 

The  mean-squared  amplitude  fluctuations,  <8^  >,  for  the  farfield 

r  K 


receiver  are  therefore 


<  B  > 
PR 


7 if  f<  B1B2  >  dtldr2' 


(4.29) 


L  ■'o 

This  result  differs  from  the  solution  for  the  nearfield  receiver 


(Eq.(4.22)]  only  in  that  <  B^Bn  >  here  is  the  correlation  function  for 
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spherically  spreading  waves  rather  than  plane  waves.  In  form  the 
solutions  are  identical. 

In  the  next  section  we  will  examine  the  components  of  the  corre¬ 
lation  term,  <  B^B2  >’  Ecls-(4.22)  and  (A. 29)  so  that  these  results 

may  be  evaluated. 

4.4  Approximating  the  Spatial  Correlation  Functions 
The  evaluation  of  Eqs.(4.22)  and  (4.29)  requires  the  calculation 
of  nine  correlation  terms;  viz . , 

<  B1B2  >  =  <[Bg(z1)  +  B  (zx)  +  B±(L  -  z1)][Bg(z2)  +  Bp(z2)  +  B±(L  -  z2)l> 
=  <Bs(z,)Bg(z2)>  +  <Bs(z1)Bp(z2)>  +  <Bs(Zl)B±(L  -  z2)> 

+  <Bp^'l)Bs(z2)>  +  +  <Bp^l)B±(L  '  z2>> 

+  <B±(L  -  z1)Bs(z2)>  +  <B±(L  -  Z;L)Bp(z2)> 

+  <B±(I.  -  z1)B±(L  -  z2)>.  (4.30) 

These  correlation  terms  fall  into  three  categories:  (1)  cross 
correlation  between  the  low  frequency  signal  fluctuations  and 
fluctuations  in  the  pump  or  interaction  frequency  waves;  (2)  cross 
correlation  between  pump  fluctuations  and  interaction  frequency 
fluctuations;  and  (3)  autocorrelation  of  the  pump,  low  frequency,  and 
interaction  frequency  amplitude  fluctuations.  These  categories  are 
discussed  separately  below. 

(1)  The  signal  wave  and  pump  wave  traverse  different  parts  of 

the  medium,  which  implies  that  the  volume  of  scatterers  associated 

with  these  waves  will  be  significantly  different.  Also,  it  has 

been  assumed  that  to  »  oj  ,  which  implies  that  there  will  be  little 
p  s’ 

frequency  correlation  between  signal  and  pump  wave.  For  these  reasons 
it  may  be  assumed  that  the  signal  and  pump  amplitude  fluctuations  will 
be  uncorrelated  in  the  interaction  region;  i. e ■ , 
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<B  (z1 ) B  (z.)>  =  <B  (z. ) B  (z„)>  =  0.  (4.31) 

slp2  p  1  s  2 

A  similar  argument  can  be  applied  to  the  signal  and  interaction 
frequency  fluctuations,  so  that 

<Bs(zl)B±(L  "  Z2)>  *  <B±(L  '  z1)Bs(z2):>  ~  °‘  (4.32) 

(2)  The  term  <Bp(z^)B+(L  -  z2)>  represents  cross  correlation 
between  pump  pressure  amplitude  fluctuations  at  the  point  (z^)  and 
interaction  frequency  amplitude  fluctuations  at  (L  -  z „) •  The 
frequencies  of  the  two  waves  are  approximately  equal  (co^  =  o)+)  so  we 
may  assume  that  complete  frequency  correlation  exists  between  B^ 
and  B  . 

The  propagation  paths  associated  with  the  pump  and  interaction 
frequency  waves  are  shown  in  Fig.  4.4.  The  term  B^(z^) 
is  due  to  scattering  of  the  pump  wave  as  it  propagates  from  the 
origin  to  z  =  Zy  The  fluctuations  in  the  interaction  frequency  wave, 
B+(L  -  z2),  are  due  to  scattering  of  the  second-order  radiation  as 
it  propagates  from  a  source  at  point  z  2  to  the  hydrophone  at  z  =  L. 

Both  and  z2  may  vary  between  0  and  L,  so  there  will  be  situations 
in  which  the  two  propagation  paths  overlap  (z^  >  z^)  and  situations 
in  which  they  are  separate  (z^  <  z2). 

An  exact  evaluation  of  Eqs.  (4.22)  and  (4.29)  would  require  an 
expression  for  the  correlation  term,  <B^B+>,  for  all  values  of  z^  and 
z2>  and  for  both  nearfield  and  farfield  receiving  arrays.  This  is 
essentially  a  problem  of  calculating  the  correlation  of  amplitude 
fluctuations  at  two  receivers  when  there  are  two  sources  generating 
separate  waves.  The  solution  to  this  problem  is  not  available  in  the 
literature,  nor  is  it  readily  obtained,  so  we  will  resort  to  the 
following  simplified  approximation. 
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PUMP  HYDROPHONE 


i 5  0  *  =  i-z  2  * s  L 


FIGURE  4.4 

PROPAGATION  PATHS  OF  PUMP  AND  SECOND- ORDER  WAVE 
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We  will  consider  two  conditions,  one  in  which  the  array  length 


L  is  much  less  than  the  longitudinal  correlation  distance  1  ,  and 

P 

the  other  in  which  the  array  length  is  much  greater  than  the  pump  cor¬ 
relation  distance.  First  assume  that  L  <<  1  .  In  this  case  the 

P 

separation  between  the  'receivers'  at  z  and  L  is  less  than  the  pump 
correlation  distance;  i.e.. 


L  -  z.  «  1  ,  (4.33) 

1  P 

where  1  is  the  distance  at  which  the  correlation  coefficient 
P 

for  the  pump  wave  amplitude  fluctuations  equals  1/e.  When  Eq.(4.33) 
applies,  the  amplitude  fluctuations  of  the  pump  wave  at  z-^  and  the 
interaction  frequency  wave  at  z  £  L  will  be  highly  correlated.  We 
approximate  their  correlation  coefficient  by  unity;  i.e.. 


<B  (z  )B+(L  -  z~)> 

i  =  - *— ! - * - - - v.*  1  ,  L  «  1 

p’~  f <Bp(z1>>  <VL  ~  *2)>l'a  P 


From  this  expression  the  cross  correlation  between  pump  and  interaction 
frequency  fluctuations  can  be  written  as 

<Bp(z1)B+(L  -  z 2 ) >  =  t<Bp(z1)>  <B^(L  -  z2)>]2  .  (4.34) 

This  equation  is  in  terms  of  mean-squared  amplitude  fluctuations,  and 
can  be  calculated  using  the  results  discussed  in  Chapter  3.  A 
similar  approximation  can  be  made  for  the  remaining  cross  correlation 
term  in  Eq.(4.30),  namely 

<Bp(z2)B±(L  -  zj>  a[<B2(z2)B*(L  -  zx)>  ]\,±, 

=  [<Bp(z2)>  <B+(L  -  z1)>]2.  (4.35) 


Again  we  note  that  Eqs.(4.34)  and  (4.35)  apply  only  for  short  array 

lengths  (L  «  lp).  In  general,  if  the  array  length  is  not  sufficiently 

short  that  L  «  1  ,  then  R  ,  is  less  than  unity,  and  we  have  the 
P  P>± 

less  restrictive  approximations 
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and 


"Bp(z2)Bi (L  "  zl)>  "  l<Bp(z9):>  <Bp(L  ~  z1)>]  2 

<B  (z1)B4  0.  -  z2)>  =  (<B^(z1)>  <B^(L  -  z2)>]  2  Rn  + 

Now  consider  the  situation  when  the  array  length  is  much  greater 


R  . . 

P> 

(4.36) 

R  . 
P.± 

(4.37) 

much 

greater 

than  the  pump  longitudinal  correlation  distance  (L  >>  Ip)*  Most 
separations,  L  -  z^,  will  be  greater  than  the  correlation  distance 
lp,  and  therefore  the  fluctuations  will  (on  the  average)  have  very 
little  correlation.  For  L  >>  1 ^  we  therefore  approximate  the  corre¬ 


lation  coefficient  bv  zero  and  obtain 


<Bn(z1)Bi(L  -  z2)>  =  0,  (4.38) 

and 

<B  (z,)B+(L  -  /  )>  =  0  ,  L  »  1  .  (4-39) 

P  -  *  1  p 

These  results  are  used  in  Eq.(4.30)  for  long  array  lengths. 

(3)  The  third  category  of  correlation  terms  in  F.q.(4.30)  are 

autocorrelation  functions  for  the  signal,  pump,  and  interaction 

frequency  fluctuations. 

The  pump  autocorrelation  function  will  be  given  by 

<B  (z  )B_(z,)>  =  (<B^(z1)>  <B^(z,)>]'i  R  ,  (4.40) 

pipz  pi  Pi  p 

where  is  the  longitudinal  correlation  coefficient  for  pump  wave 

fluctuations  at  z^  and  For  short  lengths  (L  <<  1^),  Rp  can  be 

approximated  as  unity.  For  longer  array  lengths  the  results  of 
18  94 

Chernov  or  Eliseevnin  can  be  used  to  estimate  Rp. 

Similarly,  the  autocorrelation  function  for  the  signal  wave  will 
be 

<B  (Z  +  z,)B  (Z  +  z„)>  =  [<B2(Z  +  z.)>  <B2(Zc  +  z  )>]  '2  R  ,  (4.41) 

S  S  1  S  S  2  ssl  S  S  2  S 

where  Rg  is  the  longitudinal  correlation  coefficient  for  signal 
wave  fluctuations  at  z^  and  z Here  again  the  correlation  coefficient 
can  be  approximated  as  unity  for  short  array  lengths  (L  <<  lg;  1 
is  the  signal  wave  correlation  distance).  For  L  >  lg,  theoretical 
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results  available  in  the  1 iterature^ ’ ^  can  be  used  to  predict  R  . 

s 

The  autocorrelation  function  for  the  interaction  frequency  wave  is 
somewhat  different  from  that  for  the  pump  or  signal  wave.  Rather  than 
originating  at  a  common  source  point  and  being  received  at  different  obser¬ 
vation  points,  as  are  the  signal  and  pump  waves,  the  interaction  frequency 
waves  originate  at  different  source  points  and  are  received  at  a  common 
observation  point.  This  is  sketched  in  Fig.  4.5.  The  waves  originate  at 

z ^  and  Z£,  and  are  both  received  at  z  =  L.  While  there  is  no  explicit 

91  97 

analysis  of  this  situation  in  the  literature,  Chotiros  and  Smith  ’  have 
demonstrated  that  the  principle  of  reciprocity  applies  as  follows.  If  a 
wave  of  frequency  w+  is  projected  from  the  transducer  at  z  =  L,  then  the 
correlation  of  amplitude  fluctuations  received  at  points  z^  and  z^  will  be 
<Br(L  -  z^B^CL  -  z2)>.  By  the  reciprocity  principle,  this  correlation 
will  be  identical  to  that  for  waves  originating  at  z^  and  z^,  and  received 
at  z  =  L;  i.e. , 


<B+(L  -  zx)  B+(L  -  z 2 ) >  =  <Br(L  -  z ^  Br(L  -  z2)>, 
where  the  subscript  r  indicates  that  the  positions  of  sources  and 
receivers  have  been  interchanged.  This  result  is  useful  because  it  allows 
us  to  write 

<B+(L  -  Z])  B+(L  -  z2)>  =  [<B2+(L  -  z^)  ><B2+(L  -  z2)>J  Rr  ,  (4.42) 

where  R^_  is  the  longitudinal  correlation  coefficient  for  interaction 

frequency  waves  originating  at  z  =  L  and  received  at  points  z ^  and  z2. 

When  the  array  length  is  short  (L«lp),  the  correlation  coefficient  is 

=  1.  For  longer  array  lengths,  the  results  of  Chernov^or  Eliseevnin^ 
may  be  used  to  calculate  R^.  Because  complete  frequency  correlation  has 
been  assumed  for  the  pump  and  interaction  frequency  waves,  the  coefficient 
Rr  will  be  equal  to  the  longitudinal  correlation  coefficient  for  a  wave  at 
frequency  Wp  received  at  ranges  L  -  z^  and  L  -  z2;  i.e., 


i 


FIGURE  4.5 

GEOMETRY  FOR  AUTOCORRELATION  OF 
SECOND -ORDER  WAVE  FLUCTUATIONS 


ARL:UT 
AS-80-1005 
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3.  31-80 
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The  mean-squared  fluctuation  terms  in  Eq.  (4.42)  are  for  waves  of  frequency 

iD+  traversing  paths  of  lengths  L  -  and  L  -  Z2  and  may  be  calculated 

~  19 

using  Tatarski's  formulae. 


4.5  Evaluation  of  Results 

In  the  previous  section  we  examined  the  nine  correlation  terms  that 
are  required  in  the  evaluation  of  Eqs.  (4.22)  and  (4.29).  Now  we  will 
develop  expressions  for  the  amplitude  fluctuations  of  the  second-order  wave 
for  specific  types  of  nearfield  and  farfield  receiving  arrays.  Before 
proceeding,  we  recall  that  for  small  fluctuations  the  rms  amplitude  fluc¬ 
tuations  are  approximately  equal  to  the  coefficient  of  amplitude  variation. 
Therefore,  it  is  possible  to  rewrite  Eqs.  (4.22)  and  (4.29)  in  terms  of  the 
coefficient  of  amplitude  variation  as  follows: 


CAV 


PR 


J. 

4  ^//<W  dz!dz j 


for  the  nearfield  receiving  array,  and 

L 

,2 


CAV 


PR 


J  //<B]V  drldr2 


(4.43) 


(4.44) 


i  i 


for  the  farfield  receiving  array,  where  CAVRR  is  the  coefficient  of  amplitude 
variation  for  a  parametric  receiver. 


4.5.1  Near field  Receiving  Array 

The  coefficient  of  amplitude  variation  for  the  nearfield  receiver  can 
be  evaluated  by  substituting  Eq.  (4.30)  tor  <B^B2>  in  Eq.  (4.43),  and  then 
making  the  approximations  given  in  Eqs.  (4.31),  (4.32),  (4.36),  (4.37),  and 
(4.40)  -  (4.42).  The  result  is 


CAV 


Wjj[<»p2(zl,>  <B±2(L  ‘  Rp  *  [<Bp2<*2)>  <Bi2ft  -  !1)>]\ 

+  [<»p2<zl»  <8p2(*2)>]'’  RP  +  K'1.  +  *1”  <B92fe=  *  Z2)>]H 
+  f<B+2(L  -  z,)>  <B,2(L  -  za)>]Js  R_  i  dz,  dz„ 


(4.45) 
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where,  due  to  the  high  frequency  correlation  of  the  pump  and  interaction 

frequency  waves,  we  have  assumed  that  R  =  R  . 

P>+  P 

An  example  of  the  evaluation  of  Eq.  (4.45)  will  now  be  given  for  the 
restricted  condition  that  the  array  length  is  much  less  than  either  cor¬ 
relation  length  1  or  1  (L«l  and  L«1  ).  For  this  case  R  and  R  are 
P  t>  s  p  p  s 

approximated  by  unity.  If  the  pump  and  interaction  frequency  waves  are 
planar,  and  if 

L  »  /TT  »  1 
o  p  o’ 


then  the  amplitude  fluctuations  will  be  of  the  form 

<B”>  =  0.31  C  2  k  7 76  z1176 
P  n  p  * 


(4.46) 


<B^>  =  0.31  C  2  k7/6(L  -  z)11/6 
'  n  + 


(4.47) 


We  further  stipulate  that,  for  purposes  of  illustration,  the  signal 

source  propagates  a  distance  Zg  +  z  that  obeys  the  relation 

L  »  (Z  +  z)  »  1  , 

o  s  s  o 


where  Lq  and  1q  are  outer  and  inner  scales  of  turbulence,  respectively. 


Then  the  signal  wave  amplitude  fluctuation  will  be  of  the  form 

<Bs>  ■  °-13  4  ks7/6«s  +  ‘>11/6  +  410  VZs  +  Z)3(S  -  4). 


(4.48) 


It  may  be  noted  that  Eq.  (4.48)  is  for  a  spherical  wave.  This  is  because 
the  low  frequency  signal  wave  will  propagate  as  a  spherical  wave,  although  it 
is  assumed  to  be  sufficiently  far  from  the  parametric  receiver  that  the 
wavefront  curvature  and  spreading  loss  are  negligible  in  the  interaction 
region.  Making  use  of  Eqs.  (4 . 46) - (4. 48)  and  assuming  that  =  k+,  we 
write  Eq.  (4.45)  as 


■ '  ://(•■ 


>  .  //<>r  1  1/1.’  ,1 1/17  .  _  11/12  x  1.1/12 


_  //Of 

31  >'  k  U 

n  I'  ! 


.  H/12  11/12  . 

+  z1  z  2  +  (1-  -  z  ) 


(!•  -  z.,) .  +  ?.. 

11/12  .11/12-1 
(1.  -  z.,)  J 


(I.  zx) 


■?  7/(1  ,  ,11/6  .71  .  ,3,4  4,  1 


+  [0.13  C  “  k  (/■  +  z.) 

L  ns  si 

2,  7/6,. 


+ 


[0.13  C  k  '  (Z  +  z  )  1  17  +  .  ...  . 

*-  ns's  2.  480  ms  2 


480  *m  (Zs  +  Zl}  (<t  "  Km}J 

^  K  (7,.  +  z,)3(^  -  /)]  j 


dz^dz9 


This  expression  appears  somewhat  unwieldy,  but  the  integration  is 
straightforward  and  gives 


CAVp  =  0.3375  c  2  k  7/6  L11/b  +  K  -U<a  +  bL)372  -  a372'2 
PR  “  p  L2  9b2 


J'- 


(4.49) 


where 


n  _  2  ,  7/6  11/6  |  3  it  .  4  4.3  , 

a  -  0.13  Cn  ks  zs  +  4^  ♦„(«,.  -  •  „>  Z, .  and 

.  —(0.13  C  2  k  776>;'  576  +  -|l  *  (K  «  -  „■  4;z  2 
b  =  b  n  s  s  480  m  t  m  s  . 


This  result  will  be  useful  in  Chapter  5,  where  we  will  compare  Eq.  (4.49) 
to  data  from  an  experiment  that  approximates  the  conditions  assumed  in 
deriving  this  equation. 


4.5.2  Farfield  Receiving  Array 

An  expression  similar  to  Eq.  (4.45)  can  be  derived  for  the  farfield 
receiving  array  by  substituting  Eq.  (4.30)  for  <B^B2>  in  Eq.  (4.44), 
and  then  making  the  approximations  given  in  Eqs.  (4.31),  (4.32),  (4.38), 
(4.39),  and  (4.40)  -  (4.42).  The  coefficient  of  amplitude  variation 


obtained  is 
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+  Zj)>  <B«2te» +  ^>j' 

+  f<Bp2^l)>  <Bp2<=2>]  Rp  +  K 


<‘l  <v>  <  a-v>|\ 


[ej’lL-Zj)  Bj;  (7.2)>  |'‘Rp>1  +  |<BB  (L  -Zj)>  <BB(L  -  z2)>|'iRp 


(4.50) 


where  dz  =  dr  on  the  axis  of  the  array.  We  will  examine  the  restricted 

case  where  the  pump  and  interaction  frequency  fluctuations  are  completely 

correlated  in  the  interaction  region  so  that  R  =  R  =  1  ;  we 

P  P  »  — 

assume  that  the  range  of  the  signal  source  is  much  greater  than  the  array 

length  so  the  Z  >>z.  and  Z  »z?.  This  allows  us  to  make  the 
S  I  S  t 


approximation 


<Bs2(Zs  +  2)>  4  <Bs2(zs)^  2«Zs. 


In  computing  the  pump  and  interaction  frequency  fluctuations,  it  is 
assumed  that 


L  >>  /  A  L  >>  1  , 
op  o 


so  that  the  amplitude  fluctuations  are  of  the  form 

_  2 .  »  _  ...  _  2  .  7/6  11/6  , 

<B  (z)>  =  0.13  C  k  z  ,  and 
P  n  p 

„  2,t  ,  n  „  2  ,  7/6  ,T  .11/6 

<B.  (L  -  z)>  =  0.13  C  k  (L  -  z) 

+  n  + 

The  signal  source  is  assumed  to  be  at  a  sufficiently  long  range  that 


Z  >>  L  , 
s  s  o 


so  that  the  signal  wave  amplitude  fluctuations  are  given  by 
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<B  (Z  )>  =  <)i  >  k  1  Z  a  , 

s  S>  *■  s  s 

where  a  is  Mcdwin's7^  effective  Gaussian  refractive 
distance. 

Finally,  we  approximate  the  correlation  coefficient,  R 


result 
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index  correlation 


by  Chernov's 


where  AL  is  the  separation  between  receivers. 

With  these  approximations  the  coefficient  of  amplitude  variation 
becomes 


V  TT 
2 


O 

<u"> 


dz 


2 


+  0.1415  C  2  k  7/6 
n  p 


By  making  the  change  of  coordinates 


“4  -  §  4  *3  <4 


-z^,  and  integrating,  we  obtain 

0.1415  C2  k  7/6  L11/6 
n  p 


(4.51) 


As  a  numerical  example,  consider  a  parametric  receiver  with  pump  frequency 

f  =  100  kHz  used  to  detect  a  distant  signal  source  in  the  ocean.  For  the 
P 

2  -9 

medium  parameters,  we  will  use  the  values  <y  >  =  5  *  10  ,  as  measured  by 

69  -5 

Lieberman  ,  and  a  =  0.25  m  and  =  9.30  x  10  as  computed  by  Medwin's 

methods7^.  A  plot  of  CAV  as  a  function  of  array  length  is  shown  in  Fig. 

PR 

4.6.  It  can  be  seen  from  the  figure  that  there  is  little  difference  in 

the  value  of  CAV  for  the  1  kHz  signal  when  the  range,  Z  ,  is 

extended  from  10  km  to  100  km.  Also,  at  the  longer  array  lengths  (L>50m), 

there  is  little  difference  in  CAV  when  the  signal  frequency  is  increased 

P  R 

from  1  kHz  to  5  kHz,  with  Z  =10  km.  These  results  indicate  that  CAV  „  is 

s  PK 

independent  of  f  and  Zg  at  the  longer  array  lengtiis,  at  least  for  the 
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FIGURE  4.6 

CAV  FOR  FARFIELD  PARAMETRIC  RECEIVER  OPERATING  IN  THE  OCEAN 
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values  considered  in  this  example.  It  can  be  seen  from  Fig.  4.6  that  CAV 

r  K 

is  greater  than  0.1  only  at  array  lengths  longer  than  about  200  m.  This 

means  that  an  rms  variation  in  amplitude  of  approximately  10%  is  expected 

for  the  sideband  signals  of  a  200  m  long  parametric  receiver  operating  in 

the  ocean.  This  is  greater  than  the  level  of  amplitude  fluctuations  for  a 

1  kHz  signal  received  by  a  point  hydrophone,  which  would  have,  ignoring  the 

effects  of  multipath  signals,  rms  amplitude  variations  of  approximately 

1.4%  at  Z  =  10  km  and  4.4%  at  Z  =  100  km. 
s  s 

4.6  Summary  and  Discussion 

In  the  preceding  sections  an  analysis  has  been  developed  for  the 
parametric  receiver  operating  in  an  inhomogeneous  medium.  The  principal 
assumptions  used  in  the  analysis  are:  (1)  the  signal  source  is  located  on 
the  main  beam  of  the  parametric  receiver,  (2)  there  is  complete  transverse 
correlation  of  amplitude  and  phase  fluctuations  in  the  interaction  region 
of  the  parametric  receiver,  and  (3)  the  amplitude  and  phase  fluctuations  are 
small  compared  to  unity  (i.e.,  the  medium  is  weakly  scattering).  With 
these  assumptions  we  were  able  to  derive  integral  expressions  for  the 
coefficient  of  amplitude  variation  for  the  nearfield  receiver  [Eq.  (4.22)] 
and  the  farfield  receiver  [Eq.  (4.29)].  By  making  a  number  of  assumptions 
regarding  the  spatial  correlation  functions  contained  in  these  expressions, 
it  was  possible  to  obtain  results  in  terms  of  the  parameters  of  the  medium 
in  which  the  parametric  receiver  operates.  In  this  section  we  will  discuss 
the  results  of  the  analysis  from  a  qualitative  viewpoint. 

For  purposes  of  discussion,  the  coefficient  of  amplitude  variation  for 
both  nearfield  and  farfield  receivers  can  be  written  in  the  following  fctiT 
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OAV 
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=  J  //(<-iVh  +  <Ws)  d2l  dz2* 

o 

where  <B,R,>;  -  ^Rp  (  zf  )  Bp  (z,,)  >  +  <B[)(z^)B+(L  -  z^)> 


(4.32) 


<B+(I.  -  Z1)B  (z2)>  +  <B+(L  -  z1)B+(L  -  z2)>,  and 


<B1B2>s  = 


Here  the  correlation  functions  have  been  divided  into  high  frequency  com¬ 
ponents,  <B.B0>U,  that  are  due  to  scattering  in  the  interaction  region  of 
1  L  H 

the  parametric  receiver,  and  signal  frequency  components,  <B1B2>s’  that  are 
due  to  scattering  of  the  signal  wave. 

The  high  frequency  components  increase  with  fp  and  L,  and  the  signal 
frequency  component  increases  with  fg  and  Zg .  For  sufficiently  long  ranges 
and  short  array  lengths,  the  high  frequency  components,  <B1B2>^,  will  be 
negligible  compared  to  <B1B7>S>  and  the  coefficient  of  variation  becomes 


(z1)Bg(z2)>  dz^dz2 . 


(4.53) 


It  can  be  shown  that  Eq .  (4.53)  is  also  the  coefficient  of  variation  for  a 
continuous  end-fire  array  of  length  L.  Thus  it  can  be  concluded  that,  in  an 
inhomogeneous  medium  the  parametric  receiver  performs  like  an  end-fire 
array  if  the  correlation  terms  due  to  pump  and  second-order  waves  are  neg¬ 
ligible.  It  can  also  be  shown  that,  if  the  array  length  is  very  short, 
then  the  .-m  f f i t: ient  of  variation  becomes  identical  to  that  for  a  point 


hydrophone;  viz.. 
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CAvL  =  <B  2>  (L  -  0)  . 

r  K  S 

In  a  practical  example,  such  as  the  one  discussed  in  Section  4.5.2,  the 

high  frequency  components,  <B  B„>„,  may  not  be  negligible,  and  may  even 

1  L  H 

dominate  the  value  of  CAV  .  In  this  case  the  parametric  receiver  will 

FK 

not  perform  as  well  as  a  continuous  end-fire  array,  inasmuch  as  the  para¬ 
metric  receiver  will  have  a  higher  coefficient  of  amplitude  variation. 

Having  developed  a  theory  for  predicting  the  performance  of  a 
parametric  receiving  array  in  an  inhomogeneous  medium,  we  will  discuss  in 
the  next  chapter  some  experiments  that  were  conducted  to  test  the  theory. 


CHAPTER  5. 


EXPERIMENTAL  STUDY 


_  ,  .  _  ,  _  67,84,91,98-101  ,  .  ,,  _  , 

Several  experimental  studies  of  the  effects  of  medium 

inhomogeneities  on  acoustic,  wave  propagation  have  been  conducted  using 
modelled  acoustic  and  thermal  conditions.  An  advantage  of  model  experiments 
(i.e.,  experiments  at  small  acoustic  wavelengths)  is  that  the  properties  of 
the  acoustic  medium  can  he  carefully  controlled.  It  is  possible  in  a  model 
experiment  to  minimize  the  effects  of  time-varying  surface  reflections  and 
transducer  movement,  which  produce  fluctuations  in  an  acoustic  wave  that  are 
extraneous  to  the  study.  Also,  the  characteristics  of  the  medium  in  a  model 
tank  may  be  assumed  to  be  essentially  constant  from  day  to  day,  so  that  dif¬ 
ferent  experiments  can  be  conducted  under  similar  conditions. 

In  the  present  investigation,  model  experiments  were  conducted  to 
measure  the  effect  of  medium  inhomogeneities  on  parametric  reception.  The 
investigation  will  be  described  in  this  chapter  as  follows.  First,  a  des¬ 
cription  is  given  of  the  experimental  medium  and  of  the  apparatus  used. 
Measurements  made  of  the  coefficient  of  amplitude  variation 
for  the  signal,  pump,  and  upper  sideband  waves  of  a  model  parametric  acous¬ 
tic  receiving  array  are  then  discussed.  These  measurements  are  compared  to 
theoretical  predictions  using  the  analysis  developed  in  Chapter  4.  Finally, 
the  experimental  results  are  summarized  and  discussed. 

5.1  The  Experimental  Medium 

A  water  tank  with  a  width  of  0.9  m,  length  of  1.8  m,  and  a  depth  of 
0.8m  was  used  in  this  experimental  study.  The  thermal  microstructure  was 
generated  by  an  array  of  immersion  heaters  located  along  the  bottom  of  the 
tank  (see  Fig.  5.1),  the  total  heating  power  of  the  array  being  4.5  kW.  The 
flow  of  rising  heated  water  was  broken  up  into  patches  by  a  perforated 
aluminium  sheet  mounted  about  5  cm  above  the  heated  array. 
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The  model  tank  just  described  has  been  used  as  a  research  tool  for  a 

number  of  years,  and  the  characteristics  of  its  thermal  microstructure  are 

well  documented.  During  the  course  of  this  study,  measurements  of  the 

rate  of  temperature  increase  and  the  standard  deviation  of  temperature  were 

made  to  insure  that  the  tank  was  operating  as  in  previous  studies.  The 

rate  of  temperature  increase  was  measured  using  a  mercury  thermometer  sus- 

-3  -1 

pended  in  the  water  at  mid-depth.  The  result,  1.1  x  10  °Cs  ,  agrees  to 

95 

within  10 %  of  the  results  of  Smith  and  U'eston-Bartholomew  and  Chotiros  and 
84 

Smith.  A  measurement  of  temperature  variations  was  made  using  a  thermis¬ 
tor  placed  in  the  center  of  the  tank.  Measurements  were  begun  1  hour  after 
turning  on  the  tank  heaters.  This  permitted  thorough  mixing  of  the  water 

so  that  the  thermal  activity  would  be  uniform  throughout  the  tank.  The 

standard  deviation  of  temperature  computed  from  nine  such  measurements  is 
0.023  °C,  which  is  less  than  the  value,  0.032  °C,  measured  by  Chotiros  and 
Smith, but  is  of  the  same  order  of  magnitude.  The  results  of  these  mea¬ 
surements  of  heating  rate  and  standard  deviation  of  temperature  allow  us  to 
conclude  that  the  tank  is  operating  essentially  as  reported  previously. 

91,84 

Consequently,  we  will  adopt  the  methods  developed  by  Chotiros  and  Smith 

for  describing  the  field  of  inhomogeneities  in  the  tank. 

As  discussed  in  Chapter  3,  an  inhomogeneous  medium  can  be  described  by 

a  turbulent  power  density  spectrum.  The  magnitude  of  the  spectrum,  4^,  is 

2 

a  function  of  the  mean-squared  refractive  index  change,  <p  >.  Each  of  these 
parameters  depends  upon  the  standard  deviation  of  the  temperature,  and  hence 
will  be  calculated  from  the  measurements  of  temperature  variations  made  in 
this  study.  The  boundary  wavenumbers,  k^,  and  k  depend  upon  the  physi¬ 

cal  dimensions  of  the  tank,  the  dif fusitivity,  and  the  rate  of  kinetic 

energy  dissipation  per  unit  mass.  We  will  assume  that  these  parameters  are 

34 

unchanged  since  the  investigation  of  Chotiros  and  Smith,  and  will  use 
their  values  of  k  ,  *  ,  and  k  . 

m  t’  p 
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A  list  of  parameters  used  to  describe  the  inhomogeneous  medium  in  the 
model  tank  is  given  in  Table  5.1.  These  parameters  varied  with  position 
throughout  the  tank;  their  values  at  the  center  of  the  tank  are  shown  in  the 
table.  These  values  will  be  used  as  discussed  below  in  making  theoretical 
predictions  for  the  acoustic  waves  propagating  in  the  tank. 

5.2  Description  of  Apparatus 

A  model  parametric  receiver  was  constructed  by  placing  a  pump  transducer 
and  a  hydrophone  at  mid-depth  in  the  tank,  as  shown  in  Fig.  5.1.  The  array 
length  was  varied  by  changing  the  separation,  L,  between  the  two  transducers. 
A  source  of  acoustic  signals  to  be  detected  with  the  parametric  receiver  was 
located  on  the  main  axis  of  the  parametric  receiver,  and  at  a  distance,  Z^, 
from  the  pump  transducer. 

The  apparatus  used  to  generate  and  receive  the  acoustic  waves  for  this 

study  is  shown  in  the  block  diagram  of  Fig.  5.2.  To  minimize  the  problems 

of  electronic  feedover  and  acoustic  multipaths,  both  the  signal  and  pump 

waves  were  pulsed,  the  pulses  being  typically  100  ysec  in  duration.  There 

was  a  time  delay,  t  ,,  in  the  generation  of  the  pump  frequency  wave;  this 
d 

was  to  synchronize  the  pulses  so  that  they  occur  simultaneously  in  the 
interaction  region  of  the  array.  The  pump  and  signal  frequencies  were  main¬ 
tained  at  10  MHz  and  1  MHz,  respectively,  throughout  the  experimental  study. 
Details  of  the  transducers  used  to  generate  the  signal  and  pump  waves  are 
given  in  Appendix  3. 

The  pressure  detected  by  the  hydrophone  was  separated  into  three 

components  (f  ,  f  ,  and  f  )  by  the  arrangement  of  filters  and  tuned  ampli- 
u  p  s 

fiers  shown  in  Fig.  5.2.  A  signal  processing  unit  sampled  a  10  ysec  segment 
of  each  pulse  and  detected  the  peak  value  of  this  segment.  The  amplitude 
of  each  sampled  pulse  was  then  recorded  by  a  data  logging  unit  for  later 
off-line  computer  processing.  Further  details  of  the  receiving  apparatus 
may  be  found  in  Appendix  4. 
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In  this  section  we  will  discuss  some  experiments  conducted  in  the  model 


tank  using  the  apparatus  just  described.  Since  we  are  particularly  interested 

in  the  behavior  of  the  upper  sideband  component,  p  ,  of  the  second-order 

u 

pressure  field,  we  will  first  measure  the  amplitude  of  p  in  a  homogeneous 
medium  to  verify  that  it  is  being  generated  as  predicted  by  theory.  Then  we 
will  examine  the  amplitude  fluctuations  in  the  pur.ip,  signal,  and  upper  side¬ 
band  waves  for  a  variety  of  experimental  cond: t Lons  and  compare  the  results 
to  theoretical  predictions  based  on  the  analysis  of  Chapter  4. 

5.3.1  Measurements  of  Upper  Sideband  Pressure 

Aside  from  nonlinear  acoustic  interaction,  there  are  four  ways  in  which 
voltages  at  the  upper  sideband  frequency  may  appear  at  the  input  of  the  sig¬ 
nal  processor.  When  two  pressure  waves  are  incident  on  the  face  of  the 
hydrophone,  mechanical  mixing  occurs  and  an  upper  sideband  voltage  is  gen¬ 
erated  by  the  hydrophone.  (This  extraneous  si; nal  is  called  pseudosound.) 
Also,  the  pressure  waves  may  be  received  by  the  hydrophone  and  then  mixed  in 
the  electronic  receiving  apparatus.  Another  possibility  is  that  the  two 
signals  may  be  transmitted  electromagneticallv  and  then  mixed  in  the  elec¬ 
tronic  receiving  system.  Finally,  a  harmonic  of  the  signal  source  may  be 
radiated  which  is  at  the  upper  sideband  frequency. 

Tests  were  conducted  to  ensure  that  the  11  MHz  signal  detected  by  the 
hydrophone  was  being  generated  by  acoustic  nonLinear  interaction.  The  sig¬ 
nal  disappeared  when  either  the  signal  or  pump  transducer  was  blocked  off, 
indicating  that  it  was  not  produced  bv  electromagnetic  pickup  or  by  har¬ 
monic  radiation  of  the  signal  source. 

A  measurement  of  mixing  in  the  electronic  receiver  was  conducted  using 
the  arrangement  shown  in  Fig.  5.3.  Wii h  oscillator  2  set  to  zero  output 
voltage,  oscillator  1  was  adjusted  to  provide  a  voltage  that  was  of  the 

order  of  the  hydrophone  output  voltage  when  receiving  the  pump  wave  directly. 
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Similarly,  oscillator  2  was  adjusted  to  simulate  the  hydrophone  output  when 
receiving  the  1  MHz  signal.  The  68  i!  resistors  served  to  decouple  the  os¬ 
cillators.  With  both  oscillators  on,  the  upper  sideband  voltage,  V  was 

ou  t 

more  than  50  dB  below  tile  level  observed  when  the  input  voltages  Were  gen¬ 
erated  aioustically .  This  test  confirmed  that  the  upper  sideband  signal  de¬ 
tected  by  the  hydrophone  was  not  produced  by  mixing  in  the  electronic 
rece iver . 

The  interaction  frequency  pressure  component  of  pseudosound  produced  by 

two  collinear  plane  waves  incident  at  the  hydrophone  surface  has  a  pressure 
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amplitude,  P',  given  by 


P' 


P].P 


,/p  c  * 
’  o  o 


(5.1) 


where  P^  and  P^  are  the  pressure  amplitudes  of  the  two  waves.  As  the  upper 
sideband  pressure  amplitude  measured  in  the  present  experiments  was  four 
orders  of  magnitude  greater  than  the  level  predic  ted  bv  F.q .  (5.1),  it  was 
demonstrated  that  the  contribution  of  pseudosound  was  negligible  in  these 
experiments . 

The  tests  and  measurements  described  in  the  previous  three  paragraphs 
assure  us  that  the  upper  sideband  pressure  wave  measured  in  this  study  was 
generated  by  nonlinear  interaction  of  the  pump  and  signal  waves,  and  not  by 
extraneous  effects.  Next  we  will  discuss  the  amplitude  of  this  wave.  The 
on-axis  pressure  amplitude  of  the  interaction  frequency  wave  of  a  nearfield 
parametric  receiver  is  given  by  Eq .  (2.17),  and  is 

D  _  w .  p  Pc  3  L 

FN  1  --B-  Ac -  exp  (-a  L)  .  (5.2) 

2o  c  - 

o  o 

4 

This  expression  predicts  a  pressure  amplitude  of  P  -  1.0  *  10  Pa  at  the 

hydrophone  of  the  parametric  receiver  for  the  following  conditions: 

f  =  11  MHz,  P  =  7.01  x  104  Pa,  P  =  3.85  x  104  Pa,  and  L  =  0.5  m.  This 
u  p  s 

'  ii—— 


result,  however,  Joes  not  take  into  account  the  shadowing  effect  of  the 


pump  transducer,  which  may  be  expected  to  reduce  the  effective  array  length 

,  47,51 

by  an  amount 


d  = 


1.25b" 


where  b  is  the  radius  of  the  pump  transducer  housing.  Also,  the  presence 

of  the  pump  transducer  was  observed  to  reduce  the  pressure  of  the  1  MHz  sig- 

4 

nal  wave  at  the  hydrophone  by  1.3  dB,  to  P  =  2.63  x  10  Pa.  Taking  these 

s 

effects  of  shadowing  into  account,  and  using  a  measured  value  of  attenuation 

of  a+  =  3.79  neper/m,  the  upper  sideband  pressure  theoretically  should  be 
—  3 

=6.27  x  10  Pa.  The  upper  sideband  pressure  amplitude  was  measured  to 
3 

be  P  =  6.11  x  10  Pa,  which  is  in  good  agreement  with  the  predicted  value. 
Similar  results  were  obtained  for  a  10  cm  array  length. 

We  have  determined  by  the  measurements  discussed  in  this  section  that 
the  parametric  receiver  is  functioning  as  expected  in  a  homogeneous  medium. 
Now  we  will  proceed  to  a  discussion  of  some  experiments  conducted  to  deter¬ 
mine  the  effects  of  inhomogeneities  on  its  operation. 

5.3.2  Procedure  for  Obtaining  Data 

The  procedure  used  in  measuring  amplitude  fluctuations  in  the  signal, 
pump,  and  upper  sideband  waves  may  be  described  as  follows.  The  trans¬ 
ducers  shown  in  Fig.  5.1  were  aligned  in  the  desired  geometry  at  mid-depth 
in  the  tank.  Using  the  apparatus  described  in  Section  5.2,  and  with  the 
tank  heaters  turned  off,  recordings  were  made  of  the  amplitudes  of  the  sig¬ 
nal,  pump,  and  upper  sideband  pulses.  These  measurements  obtained  with  the 
heaters  off  are  amplitude  fluctuations  produced  by  the  electronic  ap¬ 
paratus  and  by  any  'ambient'  inhomogeneities  in  the  tank  such  as  those  due 
to  microbubbles,  biological  matter,  or  residual  thermal  patches.  Typically, 


three  sets  of  data  were  recorded  with  the  heaters  turned  off,  and 
approximately  150  samples  per  data  set  were  taken  at  a  sample  rate  of  about 
0.25  Hz.  After  recording  these  data,  the  heaters  were  turned  on  and  allowed 
to  warm  up  for  1  hour  so  that  the  turbulent  mixing  could  reach  steady  state. 
Then  an  additional  five  sets  of  data,  similar  to  those  juot  described,  were 
recorded  with  the  tank  heaters  turned  on.  The  water  temperature  was  mea¬ 
sured  at  the  beginning  and  end  of  each  data  set.  When  the  mean  water  tem¬ 
perature  in  the  tank  reached  about  30°  C,  the  heaters  were  turned  off  and 
the  tank  was  allowed  to  cool  until  the  following  day.  The  reason  for  stop¬ 
ping  the  experiments  at  30°  C  is  that  the  level  of  amplitude  fluctuations 
decreases  as  the  temperature  increases;  thus  poor  1 f luctuation-to-noise ' 
ratios  usually  occurred  at  temperatures  above  30°  C.  We  will  discuss  this 
point  in  more  detail  below. 

The  measurements  made  with  heaters  on  and  off  were  analyzed  off-line 
by  a  PDP-11  computer,  using  the  computer  programme  AFLUCT,  which  is  listed 
in  Appendix  5.  The  quantities  calculated  by  this  program  are  the  mean, 
standard  deviation,  and  coefficient  of  variation  for  the  amplitude  of  each 
wave  associated  with  the  parametric  receiver.  A  high-pass  filter  with  cut¬ 
off  frequency  of  0.02  Hz  is  incorporated  in  the  program  to  minimize  the 
effects  of  slow  drifts  in  amplifier  gains  and  changes  in  received  signal 
amplitude  with  temperature. 

5.3.3  Signal  Wave  Amplitude  Fluctuations 

As  an  example  of  how  amplitude  fluctuations  of  the  acoustic  waves  were 
measured  and  analyzed,  we  will  consider  data  obtained  for  the  1  MHz  signal 
wave.  Data  were  taken  over  a  five  day  period  using  the  procedure  described 
in  Section  5.3.2.  In  all,  25  sets  of  data  were  taken  with  the  heaters  turned 
on,  and  15  sets  were  taken  with  the  heaters  off.  For  all  of  these  measure¬ 
ments,  the  separation  between  projector  and  hydrophone  was  +  L  =  0.75  m. 
Results  of  the  measurements  are  shown  in  Table  5.2.  The  coefficients  of 
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Data 

Set 

Temperature 

°C 

CAM 

on 

CAV  * 

or  t 

FNR 

dH 

i 

L 

f'AV 

1 

23.  > 

0.0077 

0.0030 

8.2 

0.7303 

0.0052 

2 

24.2 

0.0133 

0.0030 

12.9 

0.7554 

0.0098 

3 

25.0 

0.0088 

0.0030 

8.5 

0.7852 

0.0058 

4 

25.7 

0.0090 

0.0030 

9.5 

0.8123 

0.0069 

5 

26.5 

0.0069 

0.0030 

7.2 

0.8443 

0.0052 

6 

25.6 

0.0081 

0.0033 

7.8 

0.8084 

0.0060 

7 

26.2 

0.0075 

0.0033 

7.1 

0.8321 

0.0056 

8 

27.2 

0.0251 

0.0033 

17.6 

0.8734 

— 

27.9 

0.0071 

0.0033 

6.7 

0.9034 

0.0057 

28. b 

0.0072 

0.0033 

6.8 

0.9345 

0.0060 

ii 

25.1 

0.0112 

0.0043 

8.3 

0.7890 

0.0081 

12 

25.8 

0.0111 

0.0043 

8.2 

0.8162 

0.0083 

13 

26.5 

0.0088 

0.0043 

6 . 1 

0.8443 

0.0065 

14 

27.2 

0.0077 

0.0043 

5.1 

0.8734 

0.0056 

15 

27.8 

0.0120 

0.0043 

8.9 

0.8991 

0.0101 

16 

28.1 

0.0089 

0.0052 

4.7 

0.9122 

0.0066 

17 

28.8 

0.0110 

0.0052 

6.5 

0.9436 

0.0092 

18 

29.3 

0.0094 

0.0052 

5.1 

0.9667 

0.0075 

19 

30.1 

0.0091 

0.0052 

4.9 

1.0048 

0.0075 

20 

30.6 

0.0103 

0.0052 

5.9 

1.0294 

0.0092 

21 

27.9 

0.0065 

0.0036 

5.1 

0.9034 

0.0059 

22 

28.6 

0.0060 

0.0036 

4.4 

0.9345 

0.0045 

23 

29.3 

0.0057 

0.0036 

4.0 

0.9667 

0.0043 

24 

30.0 

0.0052 

0.0036 

3.2 

1.0000 

0.0038 

25 

30.7 

0.0059 

0.0036 

4.3 

1.0344 

0.0048 

*CAV  ....  is  the  average  of  3  data  sets.  Range  =  0.75  m 
ot  t 


TABLE  5.2 

DATA  OBTAINED  FOR  1  MHz  SIGNAL  WAVE 
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amplitude  variation  for  the  J  Mllz  wave  for  heaters  on  and  off  are  labelled  in 
the  table  as  CAV^  and  CAV respectively.  CAV  is  a  measure  of  amplitude 
fluctuations  in  the  voltage  recorded  bv  the  data  logging  unit  shown  in  Fig.  5.2. 
However,  these  fluctuations  may  be  produced  by  a  number  of  causes.  Both  thermal 
inhomogeneities  in  the  medium  and  other  inhomogeneities  such  as  microbubbles 
and  biological  matter  will  produce  amplitude  fluctuations  in  a  propagating 
pressure  wave.  Further  amplitude  fluctuations  will  be  produced  if  there 
are  gain  variations  in  the  electronic  system  that  generate  and  detect  a 
pressure  wave.  Therefore  we  can  say  that  the  measured  amplitude  fluctuations 
comprise  (1)  an  effect  which  we  wish  to  measure,  a  signal,  which  is  caused 
by  thermal  inhomogeneities;  and  (2)  extraneous  effects,  or  noise,  which  is 
caused  by  nonthermal  inhomogeneities  and  system  gain  variations.  It  is  rea¬ 
sonable  to  assume  that  the  noise  is  present  whether  the  tank  heaters  are  on 

or  off,  so  CAV  is  a  measure  of  the  noise.  If  the  signal  and  noise  are 
or  t 

uncorrelated,  then  they  will  add  on  a  mean-squared  basis,  and  the  desired 
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signal  will  have  a  coefficient  of  amplitude  variation  given  by 


CAV 


CAV 


off 


1/2 


(5.3) 


Values  of  CAV  for  the  present  experiment  are  shown  in  Table  5.2,  and 
are  plotted  as  a  function  of  water  temperature  in  Fig.  5.4.  In  making  this 
plot,  some  of  the  data  have  been  omitted.  Specifically,  data  with  a  fluc- 
tuation-to-noise  ratio  (FNR)  less  than  4.5  dB  have  been  omitted,  where 


FNR  =  20 


CAV 

,  on 

logl0  CAV 

off 


Low  values  of  FNR  generally  occur  at  higher  values  of  temperature;  this  is 
demonstrated  by  comparing  data  sets  21-25  in  Table  5.2  to  data  sets  1-5. 
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A  theoretical  prediction  for  the  variation  of  CAV  with  temperature  may 
be  obtained  from  the  sum  of  Eqs.  (3.30)  and  (3.32),  which  is: 


CAV  =  0.13  C 


7/6 


<Zs  + 


L) 


11/6 


ii 

480 


4>  (Z 
m  s 


i  \  ^  /  * 

D  (<t  - 


)  •  (5.4) 


It  may  be  seen  from  the  equations  given  in  Table  5.1  that  C  and  $  are 
J  n  m 

functions  of  the  mean -squared  refractive  index,  which  in  turn  is  given  by 


2  1/2 

V  - 


1  3_c 
c  JT 


9  1/2 

< (AT) ~> 


(5.5) 


where 

c  =  1403  +  5T  -  0.06T2  +  0.0003T3,  (5.6) 

=  5  -  0.12T  +  0.009T2,  (5.7) 

dT 

T  =  mean  temperature  in  °C,  and 

AT  =  deviation  from  mean  temperature. 

Substitution  of  Eqs.  (5 . 5) —  (5 .7)  into  Eq.  (5.4),  and  use  of  the  equations 

defining  C  and  $  in  Table  5.1,  gives  the  theoretical  curve  plotted  in  Fig. 
n  m 

5.4.  In  order  to  make  a  comparison  between  the  experimental  and  measured 

results,  a  regression  line  for  the  measured  data  is  also  shown  in  the  figure. 
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This  line  is  plotted  using  the  formulae  in  Crow  e_t  al.  and  is  based  upon 
the  method  of  least  squares.  It  can  be  seen  that  there  is  good  agreement 
in  both  magnitude  and  slope  between  the  regression  line  and  the  theoretical 
curve. 

The  agreement  between  theory  and  experiment  indicates  that  the 

2 

temperature  dependence  ot  <p  >,  as  expressed  by  Eqs.  (5.5)  —  (5.7) ,  is  valid. 

By  knowing  this  temperature  dependence,  we  are  able  to  compare  measurements 

of  CAV  taken  over  a  range  of  temperatures  simply  by  normalizing  all  the 
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measurements  to  the  same  temperature.  Chotiros  has  shown  the  factor, 
f  ,  for  normalizing  measurements  to  30°  C  to  be,  approximately, 


£„  “  ,OK10 


|  0.  021 (T  -  30)]. 


As  an  example,  we  calculate  the  coefficient  of  amplitude  variation  for  the 
1  MHz  signal  wave  at  the  range  0.75  m  as  follows.  For  each  data  set  in 
Table  5.1,  a  value  of  CAV  is  calculated  by  compensating  for  ambient  noise 
with  Eq.  (5.3),  and  then  normalizing  the  result  to  30°  C  with  multiplication 
by  f  .  Thus, 

[2  2  1 1/2 

CAV  -  CAv  (5.8) 

on  offl 

is  calculated  for  each  data  set.  The  mean  and  standard  deviation  of  CAV  is 
then  calculated  for  all  data  sets,  with  the  results  plotted  in  Fig,  5.5. 
Similar  results  at  other  ranges  are  also  shown  in  the  figure,  along  with 
some  curves  that  show  the  theoretical  range  dependence  of  CAV. 

The  nearfield  theoretical  curve  shown  in  the  figure  is  calculated  from 
Eq.  (5.4),  where  the  range  is  defined  to  be  Z^  +  L,  and  values  of  the  medium 
parameters  at  the  center  of  the  tank  are  taken  from  Table  5.1.  This  equa¬ 
tion  is  strictly  valid  only  for  ranges  much  less  than  Rq  where 


For  ranges  much  greater  than  Rq)  the  farfield  result  [Eq.  (3.24)]  applies. 
For  the  data  in  Fig.  5.5,  Rq  =  0.8  m,  so  some  of  the  data  are  in  a  region 
where  the  near field  theory  is  only  an  approximation.  However,  as  comparison 
between  theory  and  experiment  indicates,  the  nearfield  approximation  is  in 
reasonable  agreement  with  the  measured  results. 

5.3.4.  Pump  Wave  Amplitude  Fluctuations 


Measurements  of  the  coefficient  of  amplitude  variation  for  the  10  MHz 
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pump  wave  (CAV  )  were  made  as  a  function  of  separation  I.  between  pump 

transducer  and  hydrophone.  The  results,  shown  in  Fig.  5.6,  were  obtained 

using  the  procedure  described  in  Section  5.j .2.  The  data  were  adjusted  for 

ambient  noise  and  normalized  to  30°  C  using  Kq.  (5.8). 

In  making  a  theoretical  prediction  for  CAV^  there  are  two  factors  that 

need  to  be  taken  into  account.  One  is  that  the  mean -squared  temperature 
2 

variations,  <AT  >,  near  the  hydrophone  were  measured  to  he  different  from 
2 

<AT  >  in  the  center  of  the  tank.  This  is  due  to  the  tact  that  the  hydro¬ 
phone  was  located  near  the  edge  of  the  heater  array  and  perforated  screen 

which  produced  the  thermal  microstructure.  Because  CAV  is  a  function  of 

P 

2 

turbulence  parameters  which  depend  upon  <AT  >,  this  change  in  mean -squared 

temperature  will  affect  CAV  .  Another  factor  that  affects  the  value  of  CAV 

P  P 

is  due  to  the  finite  cross  sectional  area  of  the  hydrophone.  Chotiros^1^ 
has  shown  that  for  a  square  transducer  of  dimension  b  on  each  side,  a  re¬ 
duction  occurs  in  the  coefficient  of  amplitude  variation  of  the  amount 


CAV  2 
h 

CAV2 


1  -  0.61  b  (AL)  ’2 


(5.9) 


where  CAV^  and  CAV  are  the  coefficient  of  amplitude  variation  for  an 

acoustic  wave  received  with  a  square  hydrophone  and  a  point 
hydrophone,  respectively, 

A  is  the  acoustic  wavelength,  and 

L  is  the  separation  between  the  source  of  the  wave  and  the 
hydrophone. 

Physically,  this  reduction  in  CAV  is  due  to  the  lack  of  correlation  in 
amplitude  fluctuations  across  the  face  of  the  hydrophone.  Thus,  the  hydro¬ 
phone  measures  the  average  level  of  fluctuations  across  its  aperture.  Equa¬ 
tion  (5.9)  can  be  taken  to  be  the  definition  of  an  'aperture  factor', 
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t  -  1  -  O.t.l  h('I.) 
a 

The  coefficient  of  amplitude  variation  for  tlie  pump  wave  can  he 
predicted  by  multiplying  Eq.  (3.31)  by  Eq .  (5.10),  and  is 


0.31 


£  C2k7/I'  Llllb 

an  p 


(5.11) 


The  value  of  in  Eq.  (5.11)  is  determined  from  the  defining  equations  in 

2 

Table  5.1,  using  the  measured  value  of  <AT  >  in  the  region  between  pump 
transducer  and  hydrophone.  The  dependence  of  CAV  upon  separation  L,  as 
given  by  Eq.  (5.11),  is  plotted  as  the  solid  line  in  Fig.  5.6.  It  can  be 
seen  from  the  figure  that  there  is  reasonable  agreement  between  theory  and 
experiment. 


5.3.5  Amplitude  Fluctuations  in  the  Upper  Sideband  Wave 
In  this  section  we  will  discuss  an  experimental  study  of  the  effects 
of  medium  inhomogeneities  on  the  upper  sideband  wave  of  a  nearfield  para¬ 
metric  receiver.  Three  types  of  results  were  obtained  in  this  study. 

First,  measurements  were  made  of  the  coefficient  of  amplitude  variation  for 

the  upper  sideband  wave  (CAV  )  with  the  array  length  fixed  and  with  vari- 

USB 

able  range,  Z  +  L.  Second,  measurements  were  made  of  CAV  with  the  range 
s  USB 

fixed  and  with  variable  array  length.  Third,  a  comparison  was  made  between 
the  amplitude  fluctuations  in  the  signal  wave  and  in  the  upper  sideband  wave. 
All  data  discussed  in  this  section  were  obtained  by  following  the  ex¬ 
perimental  procedure  described  in  Section  5.3.2.  Results  were  normalized 
to  30°  C  and  adjusted  for  ambient  noise  using  Eq.  (5.8). 

The  first  set  of  results  we  will  consider  is  shown  in  Fig.  5.7.  For 
this  data,  the  array  length  was  0.1  m  and  the  range,  Zg  +  L,  was  varied  from 
0.25  m  to  1.55  m.  As  above,  circles  in  the  figure  indicate  mean  values  of 
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CAV u ^  obtained  from  several  data  sets  and  the  vertical  bars  indicate  the 

standard  deviation  of  CAV  for  the  data  sets.  In  making  theoretical 

Jon 

predictions  for  the  experiment,  two  different  methods  were  used.  One  pre¬ 
diction,  shown  as  the  dashed  line  in  Fig.  5.7,  is  based  on  Eq .  (4.49).  It 

can  be  seen  from  the  figure  that  Eq.  (4.49)  gives  values  of  CAV  that  are 

U  S  B 

significantly  lower  than  the  measured  values  for  most  ranges.  However,  the 
development  leading  to  Eq .  (4.49)  was  based  on  the  assumption  that  the  pump 
wave  can  be  approximated  as  a  plane  wave  of  infinite  extent.  A  more  ac¬ 
curate  theoretical  model,  which  takes  into  account  the  effects  of  the  finite 
apertures  of  the  pump  transducer  and  hydrophone ,  has  been  developed  by 
Cho tiros . It  is  shown  in  Appendix  6  that  the  coefficient  of  amplitude 
variation,  CAV ,  for  the  near  field  parametric  receiver  can  be  written  as 

L  L 

“Vf.R  *  l'’„l"2/7]'>«<K1)Re(K2)<B]lB2> 

O  "O  ’ 

+  Im(K1)Im(K2)<S1S2>  +  ReO^)  Im(K,)<B  S2> 


+  IM(K1)Re(K2)<S1B2> 


d  z ,  d  z 
I  z 


(5.12) 


where 


L 


dz  is  the  normalized  voltage  at  the  terminals  of  the 
hydrophone , 


9 


6  is  the  phase  angle  of  the  second  order  pressure. 
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,  and 


h  and  b  are  characteristic  radii  of  the  pump  transducer  and  hydrophone, 
respectively . 

This  result  requires  numerical  integration  for  its  evaluation.  The  pro¬ 
gram  CAVNUM,  listed  in  Appendix  5,  is  used  to  evaluate  Eq.  (5.12)  for  the 
experiments  discussed  in  this  section.  Results  obtained  by  the  program 
are  plotted  as  the  solid  line  in  Fig.  5.7.  It  can  be  seen  from  the  figure 
that  the  numerical  results  based  on  Eq .  (5.12)  are  generally  in  better 
agreement  with  the  data  than  the  analytical  result,  Eq.  (4.49). 

A  similar  set  of  results,  shown  in  Fig.  5.8,  was  obtained  for  a  fixed 
range  of  Zg  +  L  =  0.75  m,  and  a  variable  array  length.  In  this  case,  there 
is  less  difference  between  the  numerical  and  analytical  results,  although 
the  numerical  results  are  generally  in  better  agreement  with  the  data. 

Another  useful  way  of  studying  the  parametric  receiver  in  an  inhomo¬ 
geneous  medium  is  to  compare  the  coefficients  of  amplitude  variation  for 
the  upper  sideband  (CAVygg)  and  signal  (CAVg)  waves.  The  ratio,  CAVygg/CAV^, 
compares  the  'noise  level'  of  the  parametric  receiver  in  an  inhomogeneous 
medium  to  the  noise  level  of  a  point  hydrophone  in  the  same  medium.  Here 
'noise  level'  refers  to  the  random  amplitude  variations  appearing  at  the 
output  terminals  of  either  the  hydrophone  of  the  parametric  receiver  or  the 
point  hydrophone. 

Measurements  of  CAVygg  and  CAV^  were  made  simultaneously  for  parametric 
receiver  array  lengths  of  0.1  m  and  0.5  m.  The  hydrophone  of  the  para¬ 
metric  receiver  was  used  to  detect  both  upper  sideband  and  signal  waves,  as 
9hown  in  Fig.  5.2.  A  mean  value,  <CAV^gg/CAVs> ,  was  calculated  for  several 
data  sets  taken  at  each  range,  with  the  results  shown  in  Fig.  5.9.  The 
theoretical  curves  shown  in  the  figure  are  based  on  values  of  CAV()l.^ 
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from  Eq.  (4.49)  l the  dashed  lines)  or  by  numerical  integration  of  Eq .  (5. Id), 
and  on  values  of  CAV^  obtained  from  Eq .  (5.4).  Discussion  of  the  physical 
meaning  of  these  results,  as  well  as  those  shown  in  Figs.  5.6  and  5.7,  will 
be  given  below. 


5.3.6  Results  for  Spherically  Spreading  Pump  Waves 

The  experiments  described  in  the  preceding  section  were  conducted  with 

the  hydrophone  placed  in  the  nearfield  of  the  pump  transducer.  Some  data 

were  also  obtained  with  the  hydrophone  in  the  farfield  of  the  pump,  so  that 

the  pump  waves  were  spherically  spreading.  The  pump  transducer  used  in  the 

present  experiment  was  of  a  smaller  cross  sectional  area  than  that  used  in 

the  experiments  described  in  Sections  5.3.4  and  5.3.5  (see  Appendix  3  for 

a  description  of  the  transducers) .  The  experimental  apparatus  and  the 

procedure  used  to  obtain  the  data  were  similar  to  those  described  in 

Sections  5.2  and  5.3.2;  however,  amplitude  fluctuations  in  the  lower 

sideband  wave  (CAV  )  were  measured  rather  than  those  in  the  upper  side- 
LSB 

band  wave.  Details  of  the  electronic  receiving  system  used  for  this 
experiment  are  given  in  Appendix  4. 

Measurements  of  CAV  were  made  for  a  fixed  range,  Z  +  L  =  1.4  m, 

LbB  S 

and  a  variable  array  length.  Data  from  the  experiment  are  plotted  in  Fig. 

5.10.  Each  open  dot  in  the  figure  represents  a  mean  value  of  CAV  ob- 

tained  for  a  set  of  approximately  150  samples.  These  experimental  values 

of  CAV  have  been  compensated  for  temperature  and  ambient  noise  using 
LbB 

Eq.  (5.9).  The  theoretical  curve  shown  in  the  figure  is  obtained  from  Eq. 
(4.44),  which  for  the  present  experiment  can  be  written  as 
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Integration  of  Eq .  (5.13)  gives 


CAV  =  (J.  1415  C  2k  7/bx}1/b 
LSB  n  p 


,  1  4  f,  ,3/2  3/2-,' 


(5.14) 


where 


n  -  n  i  t  r  2.  7/6v  H/(>  3it  .4  4  3 

a  -  0.13  C  k  L  +  — —  4’  (k-  -  k  )Z 

ns  s  480  m  t  m  s 


,  and 


b  =  ^  (0.13C2k  1/b)7.  5/6+^-C  (r4  -  k4)Z2 
6  ns  s  480  m  t  m'  s 


It  can  be  seen  from  Fig.  5.10  that  values  of  CAV  predicted  by  Eq .  (5.14) 

LSB 


are  in  reasonable  agreement  with  the  data. 
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5.4  Discussion 

There  are  a  few  simple,  but  useful,  conclusions  that  can  be  drawn  from 
the  results  presented  in  Figs.  5.4-5.10. 

Theoretical  values  of  the  coefficient  of  amplitude  variation  for  the 
linear  waves  are  in  fairly  good  agreement  with  measured  values  (Figs.  5.4- 
5.6).  This  agreement  between  theory  and  experiment  allows  us  to  conclude 
that  the  theoretical  model  for  the  acoustic  medium  discussed  in  Section  5.1 
is  a  reasonable  one. 

Examination  of  Figs.  5.7,  5.8,  and  5.10  shows  that  CAV  tends  to 

U  bb 

increase  with  either  increasing  array  length  or  increasing  range.  This  re¬ 
sult  seems  reasonable  if  we  recall  the  physical  considerations  discussed  in 
Section  4.6.  It  was  noted  there  that  the  coefficient  of  amplitude  variation 
for  both  nearfield  and  farfield  parametric  receivers  depends  upon  two  con¬ 
tributions:  one  from  the  high  frequency  waves  in  the  interaction  region 
and  a  second  from  the  signal  wave.  As  the  range  between  signal  source  and 
parametric  receiver  is  increased,  the  contribution  to  CAV  from  the  signal 

UoD 

wave  is  increased.  This  is  simply  because  more  scatterers  are  present  in 
the  propagation  path  of  the  signal  wave  as  the  range  is  increased.  Simi¬ 
larly,  as  the  array  length  increases,  there  is  an  increase  in  the  contribu¬ 
tion  to  CAV  due  to  the  pump  and  upper  sideband  waves,  because  they 
USn 

propagate  through  an  increased  volume  of  scatterers. 

Similar  considerations  apply  to  the  results  shown  in  Fig.  5.9.  It  can 
be  seen  in  the  figure  that,  for  an  array  length  of  0.1  m,  the  amplitude 
fluctuations  in  the  upper  sideband  wave  and  in  the  signal  wave  are  approxi¬ 
mately  equal;  i.e.,  CAV  /CAV  =  1  for  long  ranges.  It  can  be  shown  theo- 

U  bb  s 

retically  that  this  is  also  true  for  the  0.5  m  array  length,  although  CAV 

USB 

and  CAVg  become  approximately  equal  at  larger  ranges  for  the  0,5  m 
long  array  than  for  the  0.1  m  long  array.  This  result  Is  reasonable,  be¬ 
cause  as  the  range  of  the  signal  source  is  Increased,  the  fluctuations  in 
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the  signal  wave  will  increase  while  the  high  frequency  fluctuations  in  the 

interaction  region  will  remain  constant.  At  sufficiently  long  ranges,  the 

contribution  to  CAV  due  to  scatterers  in  the  interaction  region  will  be- 
USB 

come  negligible,  and  only  scattering  of  the  signal  wave  will  contribute  sig¬ 
nificantly  to  CAV  .  For  array  lengths  greater  than  the  spatial  correlation 
USB 

distance  of  the  signal  wave  fluctuations,  it  may  be  expected  that  CAV 

lib  D 

would  become  less  than  CAV  at  sufficiently  long  ranges.  This  effect  would 

s 

occur  if  the  parametric  receiver  summed  uncorrelated  'noise'  throughout  the 
interaction  volume,  thus  'averaging  out'  some  of  the  fluctuations.  The  size 
limitation  of  the  model  tank  and  choice  of  signal  frequency  prevented  a 
demonstration  of  this  effect;  it  would  be  an  interesting  point  to  pursue  in 
future  research. 

5 . 5  Summary 

The  experimental  study  described  in  this  chapter  investigated  some 
effects  of  medium  inhomogeneities  on  the  parametric  receiver.  The  experi¬ 
ments  were  conducted  in  a  model  tank  in  which  an  array  of  immersion  heaters 
and  a  perforated  screen  produced  a  thermal  microstructure.  Measurements 
were  made  of  the  coefficient  of  amplitude  variation  for  the  signal,  pump, 
and  upper  sideband  waves  associated  with  a  nearfield  parametric  receiver, 
the  measurements  being  made  as  functions  of  array  length  and  of  range  from 
signal  source  to  hydrophone.  It  was  found  that  the  amplitude  fluctuations 
in  the  signal  and  pump  waves  were  reasonably  well  predicted  by  the  theo¬ 
retical  work  reviewed  in  Chapter  3.  it  was  found  that  predictions  of 
amplitude  fluctuations  in  the  upper  sideband  wave  required  a  theory  which 
takes  into  account  the  phase  variations  of  the  pump  radiation  in  the  near¬ 
field  of  the  pump  transducer.  Some  simple  conclusions  were  drawn  from  the 
experimental  results  and  from  physical  considerations:  CAV^g  increases 
with  range  and  with  array  length;  and  for  sufficiently  large  ranges,  signal 
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wave  fluc-.ruat.ions  i  easured  by  a  point  hydrophone  will  approx  imatel  y  equal 
the  sideband  fluctuations  measured  by  the  parametric  receiver. 


CHAIM  I  I;  ••  :,l'Tt\i:Y  A  c  D I  SCi.SM  1 1;  i 

(y.I  Summary  <i_i  i  lie  Thesis 

The  aim  of  the  investigation  reported  in  this  thesis  was  to  determine 
the  effects  that  medium  inhomogeneities  have  on  the  performance  of  a  para¬ 
metric  receiver.  An  introduction  to  the  study  was  given  in  Chapter  1, 
where  the  topics  of  nonlinear  acoustics  and  wave  propagation  in  an  inhomo¬ 
geneous  medium  were  briefly  discussed. 

In  Chapter  2,  the  basic  principles  of  parametric  receiving  arrays  were 
reviewed.  It  was  shown  that  nonlinear  interaction  between  an  incoming 
'signal'  wave  and  a  locally  generated  'pump'  wave  produces  an  array  of  vir¬ 
tual  sources  in  the  region  in  front  of  the  pump  transducer.  For  a  homo¬ 
geneous  acoustic  medium,  this  array  has  the  same  directivity  characteristics 
as  a  conventional  end-fire  array  realized  from  elements  operating  at  the 
frequency  of  the  wave  to  be  detected.  F.xpressions  were  given  which  describe 
the  interaction  frequency  pressure  field  when  the  hydrophone  of  the  para¬ 
metric  receiver  is  placed  in  the  nearfield  and  in  the  farfield  of  the 
pump  transducer.  Practical  considerations  regarding  processing  of  the 
difference  frequency  signal,  shadowing,  and  finite  amplitude  attenuation 
were  discussed  briefly. 

A  discussion  of  acoustic  wave  propagation  in  an  inhomogeneous  medium 
was  given  in  Chapter  3.  Brief  consideration  was  given  to  the  historical 
development  of  the  subject.  It  was  seen  that  acoustic  scattering  in  the 
ocean  is  primarily  due  to  the  thermal  microstructure,  which  produces  random 
variations  in  the  refractive  index  of  the  acoustic  medium.  Two  functions 

were  discussed  that  can  be  used  to  describe  an  inhomogeneous  medium:  the 
spatial  correlation  function,  and  the  spatial  wavenumber  spectral  density. 

Because  the  spectral  composition  of  an  inhomogeneous  medium  may  be  deter¬ 
mined  by  turbulence,  some  basic  concepts  from  turbulence  theory  were 


116 


117 


presented.  A  summary  of  expressions  was  given  for  amplitude  and  phase 
fluctuations  and  their  correlation. 

In  Chapter  4,  the  analytical  methods  discussed  in  the  two  previous 
chapters  were  used  to  develop  a  theory  predicting  fluctuations  in  the  inter¬ 
action  frequency  pressure  wave  for  both  nearfield  and  farfield  parametric 
receivers.  It  was  assumed  in  developing  the  analysis  that  (1)  the  medium  is 
weakly  inhomogeneous,  so  that  the  amplitude  and  phase  fluctuations  are  small 
compared  to  unity;  (2)  there  is  complete  transverse  correlation  of  amplitude 
and  phase  fluctuations  in  the  interaction  region  of  the  array;  and  (3)  the 
signal  source  is  located  on  the  main  beam  of  the  parametric  receiver.  By 
approximating  the  pump  wave  as  being  planar  and  collimated,  it  was  shown 
in  Eq.  (4.22)  that  mean-squared  amplitude  fluctuations  for  the  nearfield 
parametric  receiver  may  be  written  as 


<BjB,,>  dz^dz2 


(6.1) 


where  <B^B2>  is  a  convenient  notation  for  the  nine  correlation  terms 
specified  by  Eq.  (4.30).  A  similar  expression  was  derived  for  the  farfield 
parametric  receiver  [see  Eq.  (4.29)].  These  integral  expressions,  Eqs . 

(4.22)  and  (4.29),  were  evaluated  by  making  a  number  of  assumptions  regarding 
the  spatial  correlation  of  the  various  amplitude  fluctuations.  The  following 
results  were  obtained.  For  the  nearfield  parametric  receiver,  the  coef¬ 
ficient  of  amplitude  variation  is  given  by  Eq .  (4,49),  which  is 


CAV  2  =  0.3375  C  2  k  7/6  LU/6  + 
PR  n  p 
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where  a  =  0.13  C  2  k  7/6  Z  11/5  +  *  (k  4  -  k  A)  Z  3, 

ns  s  480  m  t  ms 


b  -  U(o.n  c  2  k  "V  5/6  +  ^  .  <«*  -  «  Sz  2 

6  n  s  s  480  m'  t  ms 


For  the  farfield  receiver,  the  result  is  [Eq.  (4.51)] 


2  =  k  3  a3  <y2>  Z  tan-1|-^\  +  0.1415  C  2  k  7/6  .11/6 


PR  2L  s 


n  p 


\s  /  (6.3) 

where  a  Is  the  correlation  distance  of  the  refractive  index  variations. 

The  experiments  discussed  in  Chapter  5  were  conducted  in  a  model  tank 
in  which  an  array  of  immersion  heaters  and  a  perforated  screen  produced  a 
thermal  microstructure.  Measurements  were  made  of  the  coefficient  of  ampli¬ 
tude  variation  for  the  pump,  signal,  and  upper  sideband  waves  associated 
with  a  nearfield  parametric  receiver.  It  was  found  that  the  amplitude  fluc¬ 
tuations  in  the  signal  and  pump  waves  were  reasonably  well  predicted  by  the 
theoretical  results  reviewed  in  Chapter  3.  Amplitude  fluctuations  in  the 
upper  sideband  wave  were  compared  to  theoretical  predictions  based  on  the 
analysis  of  Chapter  4,  and  to  numerical  results  based  on  a  theory  which 
uses  a  more  exact  model  for  the  pump  radiation  than  that  of  Chapter  4.  Both 
theoretical  models  gave  predictions  which  were  in  reasonably  good  agreement 
with  the  experimental  results,  although  the  numerical  results  were  generally 
in  better  agreement. 


6.2  Discussion  of  Results 

The  coefficients  of  amplitude  variation  for  the  nearfield  and  farfield 
parametric  receiver  each  contain  a  high  frequency  component  due  to  scattering 
of  the  pump  and  interaction  frequency  waves,  and  a  low  frequency  component 
due  to  scattering  of  the  signal  wave. 

If  the  array  is  sufficiently  long,  and  if  the  pump  frequency  is  much 
higher  than  the  signal  frequency,  then  the  low  frequency  contribution  will 
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be  negligible,  and  the  coefficients  of  amplitude  variation  for  the  para¬ 
metric  receivers  will  be  approximately  equal  to  their  high  frequency  com¬ 
ponents.  This  point  is  demonstrated  by  the  experimental  results  shown  in 
Fig.  5.7,  where,  at  shorter  array  lengths,  CAVp^  tends  to  'level  out'  at  a 
value  determined  by  the  high  frequency  contribution.  Also,  it  can  be  seen 
from  Fig.  5.8  that  CAV  tends  to  increase  with  increasing  array  length. 

In  the  numerical  example  for  the  farfield  receiver  (see  Fig.  4.6),  it  was 
seen  that,  for  longer  array  lengths,  the  high  frequency  contributions  were 
dominant . 

If  the  array  length  is  sufficiently  short,  and  if  the  range  of  the 

signal  source  is  large,  then  the  high  frequency  contribution  to  CAV  will 

PK 

be  of  the  form 

CAVPRi^Jjf<Bs<2l)Bs(z2>>  dzld22  • 
o 

which  is  identical  to  the  coefficient  of  amplitude  variation  for  a  continuous 
end-fire  array  of  length  L.  Therefore,  it  is  only  when  the  low  frequency 
component  of  the  fluctuations  is  dominant  that  the  parametric  receiver  per¬ 
forms  like  an  end-fire  array  in  an  inhomogeneous  medium.  Otherwise  the  value 
of  CAV  is  greater  for  the  parametric  receiver  than  that  for  the  end-fire 
array  by  an  amount  determined  by  the  scattering  of  the  high  frequency  waves 
in  the  interaction  region.  For  very  short  array  lengths,  the  parametric 
receiver  and  the  end-fire  array  will  behave  like  point  hydrophones,  so  that 

CAVPR  =  CAVEF  *  <Bs2>  (L  =  °>  • 

This  conclusion  is  supported  by  the  experimental  results  shown  in  Fig.  5.9, 
where  it  is  shown  that  the  coefficients  of  amplitude  variation  for  the 
parametric  receiver  and  for  the  'point'  hydrophone  tend  to  become  equal  at 


longer  ranges. 


120 

K”  1  >  1  ■ » '  j  •  *  J  J  >  ’  l  ll.S  gleat.M  I  ll.'lll  111.-  -.)■>;»  L  i  .1  I  I’OI  trial  j.lll  ,!l  ,  i  !  I:,, 

:  1  ’  ■ 1  * vr  I  I  in- 1  till l  inns,  hui  (ill  sui  |  icient  1  v  slm rl  lli.it  i  In-  lii.h  ;  :  - 

I  I  n*- 1  ii  it  ions  .11,-  negligible,  ,|  n,1%  be  i-xpected  that  (  W  v.mM  I.  |. 

r  K 

> 

i  linn  • .  'this  tl  I  iH't  woul,.  m  ,-ur  it  tin-  parametric  rm- 1  v<- r  -nnim-eil 

nneorre  late,!  'noise'  throughout  the  interaction  volume,  thus  ',iv<  rat-.iiir  a  ' 
some  of  the  t 1 uc tun t ions.  In  this  way  the  parametric  receiver  would  exhih i i 
the  array  gain  of  a  conventional  end-fire  array  having  an  array  1  met  1 1  gi,-it< 
than  the  correlation  distance  of  the  signal  wave  fluctuations.  inis  ettv.  • 
was  not  demonstrated  in  the  experimental  study,  due  to  the  si  :e  1  imit  iti.  ■  .  ■ 

the  model  tank  and  the  choice  of  signal  frequency:  it  would  In-  an  interest  it,, 
topic  tor  future  research, 
h  .  Cone  1  us  ions 

I  he  purpose  of  the  invest  ii.at  i « -n  reported  in  this  thesis  has  been  aehi, 
that  purpose  being  to  determine  the  ,  t  t  ei  t  s  that  thermal  inhomogene i  t  it*:,  h  ic» 
upon  the  performance  of  the  pummel  ic  receiver.  It  was  determined  hot  j 
theoretically  and  experimental  I  t  hat  t  he  thermal  microstruot  uri  gem-  tut.  e.  |  i  i, 
model  tank  produced  a  negligible  small  level  of  fluctuations  in  the  sideband 
waves  of  both  nearf it  Id  and  fat  Meld  parametric  receivers.  The  maximum  1 c  cl 
of  sideband  fluctuations  measui ed  was  0.03,  or  five  percent.  1-  was  shown 
theoretically  that  the  effects  of  a  typically  encountered  thermal,  mi  c  rest  -m  t. 
in  an  ocean  application  are  also  likelv  to  be  small.  Only  at  array  length' 
greater  than  about  123  m  (for  a  100  klU  pump  frequency)  are  the  1 luot eat  mi  - 
sideband  pressure  likely  to  exceed  0.1,  or  ten  percent.  The  main  eon, lesion 
of  this  work,  therefore,  is  that  the  thermal  microstructure  will  general  1'  ha 
a  negligible  effect  upon  parametric  reception. 

[’here  are  several  other  conclusions  that  can  be  drawn  from  tin  sttidv. 


reconciliation  of  experimental  and  theoretical  results  provides  evidence  i h  it 
the  spectral  description  of  the  medium,  as  formulated  bv  Choi i rem  and  Smith, 
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is  correct.  Furthermore,  as  this  spectra]  description  is  based  upon  Pao's 
model  of  the  dissipative  subrange,  the  present  results  provide  further  evidence 
to  support  Pao's  model. 

In  comparing  results  for  the  nearfield  and  farfield  parametric  receivers, 

it  may  be  concluded  that  the  thermal  microstructure  will,  with  all  other 

conditions  being  the  same,  have  less  effect  on  the  farfield  receiver  than  the 

nearfield  one.  This  assumes  that  the  condition  b  J  \  L  is  sari;  Hid,  so 

o  v  p 

that  the  results  in  sections  4.5.1  and  4.5.2  apply. 

It  was  determined  that  varying  the  array  length  of  a  parametric  receiver 
has  two  effects  upon  its  performance  in  an  inhomogeneous  medium.  firstly,  if 
the  array  is  long  compared  to  the  correlation  distance  of  the  signal  wave  that 
is  to  be  detected,  then  the  array  will  tend  to  reduce  the  level  of  the  signal 
wave  fluctuations  by  'averaging  them  out'  over  the  array.  In  this  regard,  the 
parametric  receiver  performs  like  a  conventional  end-fire  array  of  the  same 
length.  Secondly,  the  fluctuations  of  the  pump  and  interaction  frequency  waves 
increase  with  array  length,  so  the  performance  of  the  parametric  receiver  will 
be  degraded  in  comparison  to  that  of  the  end-fire  array  at  very  long  array 
lengths.  If  the  randomizing  effect  of  the  medium  on  the  incoming  signal  has 
already  been  large  by  the  time  the  signal  reaches  the  parametric  receiver, 
then  the  signal  wave  fluctuations  will  dominate  the  medium  effects,  and  there 
will  be  little  difference  between  the  performance  of  the  parametric  receiver 


APPENDIX  1.  DERIVATION  OF  EQUATIONS  (4.20)  and  (4.21) 


The  integral  we  wish  to  evaluate  is  the  transverse  component  of  the 
second  order  pressure  in  the  interaction  region  and  is  given  by 

b  b 

I  =  I  I  ~  exp(-jk  r)dxdy,  (Al.l) 


-b  -b 

where  r  =  [ (L  -  z)2  +  x2  +  y2]5 
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Using  a  procedure  developed  by  Berktay,  we  will  find  an  approximate 
expression  for  I,  cast  that  expression  into  the  form  of  Fresnel  integrals, 
and  then  use  asymptotic  values  of  the  Fresnel  integrals  to  obtain  a  final 
result. 

In  Eq.  (Al.l)  make  the  following  approximations:  in  the  1/4  term,  use 
r  =  L  -  z;  and  in  the  argument  of  the  exponential  term,  use 


,  2  .  2 

.  T  .  1  x  +  y 

r  =  L  -  z  +  -z  ■ 

L  L  -  z 


Making  use  of  the  symmetry  in  x  and  y,  Eq.  (Al.l)  can  then  be  written  as 


2e 


-jk+(L  -  z)  b 


L  -  z 


u 

J"  exp[-jk+x2/ (I.  -  z)]dx J" exp [-j k+y2/ (L  -  z)]dv 


(A1.2) 


This  expression  contains  two  integrals  of  the  form 


l 

J 


F  =  I  exp[-jk+u  (I.  -  z)  ]du 

w 
O 


(A1.3) 


Making  the  change  of  variable, 


/  if  (L  -  z) 

u  =  v/  - ; - 

k. 
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12  3 


Eq.  (A1.3)  becomes 


where 


F , 

o 

b'  =  b[w(L  -  z)J  •  The  integral 


exp(-j|  v2]dv  , 


(A1.4) 


in  Eq.  (A1.4)  is  the  Fresnel 


integra’  and  in  the  nearfield  of  the  pump  transducer  the  following  asymp- 

„  ,  .93,47 

totic  value  n iy  be  used 


lim 

b 


/ 


exp[-j~  v2]dv  =  -|(1  -  j) 


(A1.5) 


Substitution  of  Eq.  (A1.5)  into  Eq.  (A1.4)  and  using  the  result  in  Eq. 
(A1.2)  gives 


1 


-ik  (L  -  z) 

2e  _ 

(L  -  z) 


/tKL  - 


* 


or 


I  = 


-jk±(L 

e 


z) 


(A1.6) 


This  appears  in  the  text  as  Eq.  (4.20).  A  similar  result  can  be  obtained 
for  the  normalized  homogeneous  component,  H,  of  the  observed  pressure. 

The  definition  of  H  is  given  in  Eq.  (4.16)  and  is 


H 


+  r') 


This  can  be  written  in  terms  of  I  as 


L 


o 


(A1.7) 


Substitution  of  Eq.  (Al.b)  for  I  gives 


R 


which  appears  in  the  text  as  Eq.  (4.21). 


APPENDIX  2.  TRANSVERSE  CORRELATION  OF  FLUCTUATIONS 

One  of  the  assumptions  made  in  the  theoretical  analysis  of  Chapter  4 
was  that  there  is  complete  transverse  correlation  of  amplitude  and  phase 
fluctuations  in  the  interaction  region  of  the  parametric  receiver.  We  will 
discuss  this  assumption  for  both  the  nearfield  and  farfield  parametric  re¬ 
ceiver  in  this  appendix. 

A2.1  Transverse  Correlation  for  the  Nearfield  Case 

For  a  plane  wave,  the  volume  of  inhomogeneities  that  contribute  to  the 
scattered  pressure  is  a  cone  with  its  vertex  at  the  observation  point  and 
an  aperture  angle  of  1/ka  radians,  where  a  is  the  radius  of  the  inhomogenei- 
ties.  If  the  medium  contains  a  distribution  of  various  sizes  of  inhomo¬ 
geneities,  there  is  a  cone  associated  with  each  size,  the  largest  cone  being 
associated  with  the  smallest  inhomogeneities.  At  each  point  in  the  inter¬ 
action  region  of  a  nearfield  parametric  receiver,  therefore,  only  the 
scatterers  contained  within  their  associated  cone  contribute  to  the  ob¬ 
served  scattered  pressure.  Furthermore,  only  scatterers  insonified  by  the 
collimated  pump  beam  contribute  to  the  observed  pressure.  The  volume  over 
which  the  fluctuations  are  assumed  to  be  transversely  correlated  is  the 
intersection  of  the  cone  and  the  insonified  region,  as  shown  by  the  diag¬ 
onal  lines  in  Fig.  A2.1. 

As  an  example,  for  the  model  tank  used  in  the  present  study,  most 
patches  contributing  to  the  acoustic  fluctuations  will  be  larger  than 

_3 

a  =  2it/ko  =  4.19  x  10  m.  This  means  that  the  largest  cone  will  have,  for 

-3 

pump  waves  of  10  MHz,  an  aperture  angle  of  5.7  x  10  rad.  For  an  array 
length  of  L  *  0.5  m,  it  can  be  shown  that  the  width  of  the  significant 

-3 

scattering  volume  in  the  middle  of  the  array  (at  z  =  0.25  m)  is  1.43  x  10 
m.  It  can  be  found  from  Tatarski's  results  (Fig.  13  of  Ref.  19)  that  the 
correlation  coefficient  for  the  maximum  transverse  separation,  1.43  mm,  is 
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VOLUME  ASSOCIATED 


FIGURE  A2.1 

SIGNIFICANT  SCATTERING  VOLUMES 
IN  NEARFIELD  OF  PUMP 


ARL:UT 
AS-80-1 120 
CRC  •  GA 
4-28-80 
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R  =  0.65.  The  transverse  correlation  coefficient  for  a  9  MHz  plane  wave 
B 

84  T 

at  a  range  0.255  was  measured  in  the  model  tank  to  be  R  =  0,?  at  a  sepa- 

D 

ration  of  1.4  mm.  These  results  indicate  that  the  assumption  of  complete 
transverse  correlation  of  the  pump  amplitude  fluctuations  is  reasonable  for 


array  lengths  l  <  0.5  a.  A  comparison  of  the  correlation  coefficients  Rfi 

and  Rg  [see  Eqs.  (3.38)  and  (3.39)  or  their  plot  in  Fig.  3.3.1  of  Ref.  91] 

T  T 

indicates  that,  for  0  <  (Ap/a)  <  z,  R  _>  R  .  This  means  the  phase  fluctua- 

S  B 

tions  will  be  correlated  over  a  larger  volume  than  the  amplitude  fluctuations. 


A. 2. 2  Transverse  Correlation  for  the  Farfield  Case 

For  a  spherical  wave,  the  significant  scattering  volume  is  an 

ellipsoid  with  foci  at  the  points  where  the  source  and  receiver  are  lo- 
68 

cated.  The  maximum  transverse  dimension  of  the  significant  scattering 
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volume  for  a  scatterer  of  radius  a  is  given  by  Aiken  to  be  equal  to 
L  tan[ (4tca)-1] .  For  a  point  source,  the  condition  of  having  complete 
transverse  correlation  of  fluctuations  is  therefore  approximately  given  by 


L  tan 


< 


where  lfc  is  the  correlation  distance  of  the  fluctuations;  i.e.,  the 

T 

separation  at  which  R  equals  1/e. 

For  a  source  having  a  finite  aperture  and,  hence,  directivity,  the  volume  of 
scatterers  which  contribute  to  acoustic  fluctuations  is  the  intersection  of 
the  significant  scattering  volume  and  the  volume  insonified  by  the  source 


beam  • 


APPENDIX  3.  TRANSDUCERS 

The  transducers  used  in  the  experimental  study  are  described  in  this 
appendix. 


A.  3.1  Signal  Source 

The  signal  source  transducer  consisted  of  a  2.0  cm  square  ceramic  plate 
mounted  on  a  backing  made  from  a  mixture  of  flyash  and  epoxy  resin.  The 
transducer  was  operated  at  its  resonant  frequency  of  1  MHz.  Its  efficiency 
was  measured  with  a  radiation  balance,  and  found  to  be  37%.  The  3  dB  beam- 
width  was  measured  as  3.2°. 

A. 3. 2  Pump  Transducer  for  the  Nearfield  Receiver 

The  pump  transducer  used  with  the  nearfield  parametric  receiver 
consisted  of  an  air-backed  1  cm  circular  ceramic  plate,  resonant  at  10  MHz, 
and  mounted  in  a  brass  housing.  The  diameter  of  the  housing  was  2.0  cm. 

A  radiation  balance  was  used  to  measure  the  transducer  efficiency,  which 
was  found  to  be  3.4%.  Excess  attenuation  of  the  pump  wave  at  ranges  of 
10  cm  and  50  cm  was  shown  to  be  negligible  by  demonstrating  that  the  pump 
pressure  increased  linearly  as  a  function  of  driving  voltage  applied  to  the 
pump  transducer.  The  3  dB  width  of  the  pump  column  was  measured  to  be 
0.38  cm  at  10  cm  range,  and  0.75  cm  at  50  cm  range, 

A. 3. 3  Pump  Transducer  for  the  Farfield  Receiver 

For  the  farfield  parametric  receiver,  the  pump  transducer  was  a  0.5  cm 
square  ceramic  plate  backed  by  epoxy  resin.  The  efficiency  was  determined 
from  admittance  measurements  made  in  air  and  water,  and  was  measured  to  be 
3%. 


A.  3.4  Hydrophone 

The  hydrophone  used  in  the  experiments  was  a  2  mm  square  ceramic  plate 
backed  with  epoxy  resin.  Its  sensitivity  at  the  frequencies  of  interest  are 
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given  in  Table  Al.l.  The  sensitivity  was  determined  by  placing  the  hydro 
phone  in  a  pressure  field  of  known  amplitude  and  measuring  the  voltage  ap 
pearing  at  the  hydrophone  terminals. 


f 

S 

MHz 

mV/ Pa 

1 

8.37  x  10-5 

9 

3.7  x  10“3 

10 

1.31  x  10-3 

11 

1.15  x  10-3 

f  =  frequency 
S  =  receiving  sensitivity 


TABLE  Al.l 


HYDROPHONE  SENSITIVITY 


APPENDIX  4.  ELECTRONIC  RECEIVING  SYSTEM 

A  block  diagram  of  the  electronic  receiving  system  used  with  the 
nearfield  parametric  receiver  is  shown  in  Fig.  A4.1.  The  hydrophone  is 
transformer  coupled  to  the  filters  and  tuned  amplifiers  that  are  used  to 
separate  the  various  signals. 

A  single  channel  of  the  signal  processing  unit  is  sketched  in  block 

form  in  Fig.  A4.2.  For  the  'upper  sideband'  (USB)  and  1  MHz  'signal'  chan- 
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nels,  a  peak  detector  based  on  Knott's  design  was  used.  For  these  two 
channels,  a  sample  and  hold  circuit  was  constructed  that  uses  a  National 
Semiconductor  LF398  integrated  circuit.  The  peak  detector  and  sample  and 
hold  circuits  used  for  the  'pump'  channel  had  been  used  in  previous 
studies  and  are  described  in  Ref.  110.  The  data  logger  was  a  commercial 
unit  made  by  Solartron. 

The  gain  of  each  channel  at  the  frequencies  of  interest  is  shown  in 
Table  A4.1.  The  voltage  transfer  characteristic  of  the  receiving  system 
was  determined  by  replacing  the  hydrophone  with  a  voltage  source  and  mea¬ 
suring  the  output  voltage  for  each  channel  of  the  signal  processing  unit 
as  a  function  of  input  voltage.  Each  channel  was  determined  to  be  linear 
over  the  voltage  ranges  used  in  the  experiments. 

A  block  diagram  of  the  receiving  system  used  with  the  farfield 
parametric  receiver  is  shown  in  Fig.  A4.3.  An  inductor  Lq  was  placed  in 
parallel  with  the  hydrophone  to  tune  out  the  static  capacitance  of  the 
ceramic  plate  at  the  lower  sideband  frequency,  9  MHz.  The  low-pass  filter 
shown  in  the  figure  attenuated  the  pump  frequency  voltage  by  approximately 
40  dB,  and  attenuated  the  9  MHz  voltage  by  about  11  dB.  The  variable-gain 
tuned  amplifier  was  operated  typically  at  a  gain  of  50  dB,  with  a  bandwidth 
of  approximately  70  kHz,  centered  about  9  MHz.  The  signal  processing  unit 
was  that  described  in  Ref.  110.  It  selected  a  10  to  35  ysec  section  of  the 
received  pulse,  the  section  being  sampled  about  40  psec  after  the  beginning 
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FIGURE  A4.2 

ONE  CHANNEL  OF  SIGNAL  PROCESSING  UNIT 


ARL:UT 
AS-80-1121 
CRC  -  GA 
4-28-80 
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of  the  pulse.  The  sampled  pulse  was  peak-detected  and  converted  to  digital 
form  by  the  data  logger. 


Channel 

Frequency 

Gain 

MHz 

dB 

1 

+  48 

1  MHz 

10 

-  48 

Signal 

11 

-  45 

1 

-  32 

10  MHz 

10 

+  11 

Pump 

11 

-  12 

1 

-  60 

11  MHz 

10 

-  49 

USB 

11 

+  24 

TABLE  A4.1 


FREQUENCY  RESPONSE  OF  SIGNAL  PROCESSING  UNIT 


APPENDIX  5.  COMPUTER  PROGRAMS 


This  appendix  contains  listings  of  the  computer  programs  used  in  this 
study.  All  programs  are  written  in  the  FORTRAN  language.  The  program 
AFLUCT  was  written  for  a  PDP  11  computer.  Program  CAVANA  and  CANVUM  were 
written  for  a  CDC  Cyber  171  computer. 
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*■ 
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C***<**(*<H»<**«O**«<*»<'0;r'>  ‘>9»*<»*>'«4l«‘»'>«l*<»<H*i>***<*****®®**»1H****4*** 

c  p  1, 1 1  ■  ’  f<  r>  M  A  M  - 1 )  I  I 

C  IhIS  PROGp AM  CAl-tMi  ATFS  ThF  MFAN,  STANDARD  OFVtATlnNt 

C  COF.RFlClt.M  Of-  AMr’LlTUOt  VARIATION*  AND 

C  au r olohrel a  f  i on  h  iiw  data  me.  ad  f Mr m  paper  taPf. 

C  provision  is  mam  f  or  t  -  b  chaunEi  s  of  data. 

C  a  HIGH  pass  E  IL  1EP  WITH  VaRIahiE  CllTOFF  FRFQijEmCY 

C  Is  PROVIDED. 

C 

c  Set  up  arrays. 

DIMENSION  RLABEL  (s>  *  I  LABEL  (a  1  *AU1  <3>  »AO?  <  "U 
D I  mens  I un  OATA<b»S)*CROS(b). CPRO (5*5) » CO V ( s *  5 ) » i *  < S I 
DIMENSION  RMf  AN  (  5  t  5  >  ♦  SL»S  (  s  )  .CvS(5)  *RMX(b) 
dimension  Sox (s> *c^x (S) »std (S) 
c  Enter  data  into  arrays, 

DaT  a  PLABFL/XmE AN* ♦ *SOVN*  « *r OVA* ♦ X A uc OF  t*rRcO*/ 

DATA  ILAHFE/F  *,* F  *«*S  *«*RF*/ 

DATA  ADI / *  1  * l  * .* \ *2  *  *  *  )  *  3  * / 

DATA  AD2/*l*2  ***?*3  *»*3*[  */ 

C  start  here  a t  beginning  of  each  *apl. 

1  CONTINUE 

REWIND  t 

C  HEAU  IN  Fit.  TER  CONSTANT 

WRITE  <b,  MO' 

300  format  </./*  number  points  In  filter*/) 

RE  AO  <6*3  1)  IFLTR 

301  FORMAT  ( I i ) 

IFLl=IFLTP-t 

c  Read  in  number  of  channels  with  data 
write  (Of  302) 

302  FORMAT  I //fNudhEP  OF  CHANNEL G  wITh  DATA*/' 

READ  (6*3  .3)  NCH 

303  FORMAT  (HI 

c  set  array  i lements  =  o 

DO  3  I  =  1 T  nCh 
Sl)S(l>=  . 

C VS  <  I  )  =  . 

CROS (  I) =0. 

DO  3  J= I » IF  L  TP 
DAI  A ( I » J) =0.  i 
RMEAN  (  I  *  J)  =  ). 

3  COV<I*J>*£. 

EN*i)  . 

N=0 

ENF=IFlTR 


C  read  in  tape  m  and  number  of  points  to  be  anai  yzEd 

c 

WR I Tf  (6*7  ' ) 

70  FORMAT!//*  ID  */) 

READ (6*Ho ) 

80  FORMAT(X) 

4  WRITE (6*7) ) 

71  FORMAT!//*  NUMBER  OF  POINTS  (USE  I3l*/I 
R£AU ( 6  »  B 1 )  INSP 

81  F nRMA  T ( 1 3  t 
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C 

C  shift  data  through  array*  read  i*  nlw  data 

c 

5  CON T  1 1 JUt 

DO  7  1*1, u C h 
DO  7  J=1,IFlT) 

KsIFLTR-J 

7  DATA  ( I  .«.♦!  )  =paTA  U  »K) 

READ(5»lfll  *END=19*ERR=l9)  D/>  T  A  ( 1  »  1)  ,  I X  <  1  )  . DaT  A  <2 »  1 )  , 
2  X  (2)  *OATA  (3,n  ,  IX(3>  .DATA  (**,1)  .IX(4)  ,DATA  <5,1)  »IX  (5) 
101  FORMAT(f-5. » I3.F7.3.  I  T»F7.3,T3,F7.3,I3) 

C 

c  test  for  end  of  tape. 

6  DO  t)  I  =  1,nCh 

IF(DATA(1  ,1)  .L'O.O.)  00  TO 

c  test  for  had  data 

IF (  lx (  1  )  .nE.V  )  GC  To  So 
IF (DATA < I f  1  )  .GT. 1 ,9«)  GO  TO  5' 

IF (DATA  (  I  ,  1  )  .LT.0.C65)  GO  Tr  50 

8  CONTlNUt 

IF(N.GT.O)  GO  TO  16 
C  print  first  data 

WRITE<6*100)  (DATA (  I  » >  »  , I  =  | ,NCH) 

100  FORMAT ( 1 1 h  FIRST  DATA, 1 X.SF 7.3) 

16  N=N* 1 

IF(N.LE.IFLTl )  GO  TO  5 
C 

c  collect  uaTa  imIo  sums  and  sums  cf  squares  and  products 

c 

En=En* 1 
DO  15  I=J,NCH 
RMX  « I ) =0. 

SDX ( I ) *L , 

no  14  j=i , ifltr 

RMtAN ( I . J ) =RM(  AN ( I . j) *0414 ( I » J> 

SDX  ( I  )  =SDX  (  I  )  ♦  (DATA  <  I  »  J)  **<•>) 

RMX ( I ) =HMX l I) .DATA  « I . J) 

COV ( I » J) =COV ( I »  J ) ♦ ( DA  T  A ( I , ) ) *DATA< I  * J)  ) 

CPRD ( 1 » J) =CPRO ( 1 , J) ♦DATA ( 1 . j) *DATA (2, J) ) 

CPRO ( ?. J) sCPRO ( 2  »  J ) ♦ (DATA (2, J) *DATa (3. J) ) 

CRRO ( 3, ) =cRRO ( 3. J) ♦ (DAT  A ( 3. J  *  *DATA ( 1  ,  J)  I 

14  CONTINUE 

RMX ( I ) »HMX ( I ) /ENF 
RMSsRMX ( I ) **2 

sdx(i)=sdx(d/enf 

sdx(i)=sox(d-rms 

SdS(I)*S[)S(I)*SUX(1) 

I F ( RMS , t Q . C . )  GO  TO  15 
CVS< I ) «LVS< I ) ♦ (SDX  < I ) /RMS) 

15  CONTINUE 

IF CN.LT.INSP)  GO  TO  S 
GO  To  2 

1'?  INSPso 


i  (s 


c  process  data 

?0  CONT  J  NUF 

WR I T  F  (  6  * )  i  )  •  i 

102  F'nRMA  T  (  1  X  »  Is  SaMR|ES//> 

no  rj.  ut,tM 

C  C Al-CUl  a  rt.  Tot  mE,ui  ot-  At  l  hUT  TH:  FIRST  TWO  data  POINTS 
’  >  30  J = 1  »  1  F  I  TP 
RME  AN  (  1  »  J  )  -WM|  am  I  ,  J  )  /LN 
30  CONT  I  Nllt 

C  Calculate  standard  deviation  for  All  but 

C  Tnt  FIRST  FOUR  UaTA  points 

STD ( I)  =  (COV  t I . 1  ) /EN) - (RMFAU  !  I. n #*2> 

STD ( I >  =SOR 1  (STD ( I )  ) 

C  Calculate  coefficient  of  van.  u/-  filter  for 

c  all  out  first  four  points 

RMX ( I ) =STD  t I ) /RMEaN (  I » I  ) 

c  filtered  standard  deviation 

SDX  (  I  i aSQRl  (SOS (  i ) /e ) 

c  Filtered  coefficient  of  variation 

23  CVX  (  I  )  3SQRT  (CVS  ( I ) /E-1) 

WRITE  <6*  It'*! 

104  FORMAMOH  CHANNf.L  ♦  6* » 2M  l,8*t?H  2*8X*?H  3.8X* 

2  2h  4  ,  H X  *  2h  S/> 

WHITE (6*1 j3)  RLABEL ( !  >  *  I  LARI  L ( 1 ) ♦ <RMEftN< I • I  ) • 1*1 »NCH> 
WRITE  ( fe  *  1 0  3 )  RLABEL (?)  *ILAMlL(l» <STD<I> • I S 1  ,NCH) 

WRI Tfc  (6*1  )3»  RLABEL  <3) » I LA'U  L ( 3 ) »  < RMX ( I ) * T a \ »NCH» 

WRl Tfc (6* 1  3)  RLABEL <  2  >  »  1 L  A  Hi  L  ( 2 )  »  ( SDX  ( I )  ,  T  s  | ,NCH» 

WR l Tfc ( 6  » 1  3)  RLABEL < 3) t ILAHI L (4) * (CVX ( I ) ♦ T= 1 .NCH) 

103  FORMAT < tX, Ab*A?*5Fl  >.*» 

c  Calculate  autocorrelation  function  a < i *  1 > *x <  i *  n 

DO  24  J=1  ,  1  F l  TH 
DO  2b  1=1 *NCH 

RMX  (  1  )  a  (  ((  ov  (  1  ,  ,J)  .'£N1  -  IRMfAt  (1*1)  *RMfc  A  N  (  I  .  J  )  )  ) 

2  /<STD(I»»*2) 

25  CONT  INlJt 

WR I Tfc (6*1  3)  RLARfcL <4) » ILAH'L ( J) * (RMX ( IrNr-) 

CONT I NUt 

IFUNSR.fco.il)  GO  TO  1 
GO  TO  4 

C  IF  fcHRORS  ARE  PRESENT  IN  D A ( A  *  CONTROL  COMES  HFRf 
50  WRITE (6»2  l )N 

WRl Tt (6*99! DATA ( 1 , 1 )  *  I  A  (  j  >  *!  ATA (2* 1  )  ♦ IX (2»  ,DATA (3* 1 ) 

?  *  I  X ( 3 ) »  D  A  T  A ( 4  *  I )  *  1  A  I  4 ) ,  U  A  T  A (5*1) *  I  X  (  5  ) 

READ (6.20? I  DATA  <  1  ,  1  !  . I  A  < 1 ) ,r ATA (2* 1 > , IX (2> .DATA (3» 1 ) 

2  *  I  X ( J )  *UATA(4*j)  *  I  X  (  4 )  ,  0  A  T  A ( S  *  1  )  *  I  X  (  5  • 

GO  TO  6 

99  FORMAT  UX,Fb.l»I3,F7.3*Il,F7.3*I3*F7.3*Ii.F7.3»I32) 

20 2  FORMAT  (F5.J*I3»F7.3*J3*F7.3.I3*F7.3*I3»F7.3*H> 

201  FORMAT (6H  ERROR, 15) 

END 


o  oooooonooo  n  o 


C 


PROGRAM  CAVANA  (INPUT. OUTpUT) 

******#•*<♦*#**#*##*•*##***#0# 


REAL  KP»KS*PHI »KT,KM,LtLF.MUP*MUS 

COMMON/ A/  Kp,KS*PHl  »f<T»KM,L«LF»MUP»MUS*ZS.RNG»PI 

COMMON/a/  CAV*CSI0»CPUMP 

KP*4.1691E4 

KS*4.1691E3 

DATA  PHI.kT.kM  /2.8E- 1 2. l ?9 . 0 .22 .9/ 

PI  *  3.14159265 

THIS  PROGRAM  CALCULATES  The  COEFFICIENT  OF  AMPLITUDE 
VARIATION  FOR  SIGNAL. PUMP*  and  USB  WAVES  USING  ANALYTICAL 
EXPRESSIONS  DISCUSSED  IN  CHAPTERS  3  AND  4. 

the  turbulence  parameters  are  varied  with 
POSITION  in  the  tank. 

AN  APERTURE  CORRECTION  FACTOR  IS  USED  To  COM¬ 
PENSATE  FOR  THE  FINITE  SIZE  OF  THE  HYDROPHONE. 

SET  RANGE  AT  RNg=0.75M  AND  VARy  THr  ARRAY  LENGTH. 


RNG=0 . 75 

c  print  heading  for  results  obtained  with  fixed  range* 

PRINT  1  > 

10  FORMAT  (/.16X.WTHE  FOLLOWING  A«E  FOR  RAnGFs-,75  MW) 

PRINT  2” 

20  FORMAT  (/.22X.WLW*10X.WCPljMpW,I0X.WCAVUSBW) 

C  VARY  ARRAY  LENGTH  IN  INCREMENTS  OF  0.05M. 

DO  30  1=1.14 
RI  =  I 
L*  RI*  .05 
ZS  =  RNG-L 

C  CALCULATE  CaV  For  PUMP.UsB.AND  signal  waves, 
call  CAVSUB 
PRINT  21 .L*CPUMP,cAV 

21  FORMAT  </»16X.F8.2»2f15.6) 

30  CONTINUE 
C 

C  SET  array  LENGTH  AT  L«0.iM  AND  VARY  RANGE, 

c 

PRINT  11 

U  format  <////,*  the  following  are  for  l  *  n.i  m#i 

c  print  headings  for  these  results. 

PRINT  25 

25  FORMAT  </.20X,*RanGE*.8x,*Cs1gP»12X.*CAv*) 

L  *  0.1 

c  vary  range  in  increments  of  o.im. 

DO  31  1*1.24 
RI  *  I 

RNG*  RI*  0 . 1 «q  .  1 
ZS  *  RNG-L 

c  calculate  cav  For  pump.usb.and  signal  waves, 
call  cavsub 

PRINT  21 .RNG.CSIG.CAy 

31  CONTINUE 

C  SET  ARRAY  LENGTH  AT  L*0.5M  AND  VARY  RANGE. 


o  o  o  o  o  r> 


HO 


c 

PRINT  111 

111  FORMAT  <////,*  ThE  FOLLOWING  ARE  FOR  L  *  0.5  M# » 

c  print  headings  for  these  results. 

PRINT  125 

125  FORMAT  < / > 20X « FRANCE* * RX . *CS I GF » 1 2X » FC A V* ) 

L  *  0.5 

C  VARy  RANGE  in  INCREMENTS  OF  P.1M, 

DO  131  1=1 *20 
R I  ■  I 

RNG*  RI*  3 •  1 *0 •  5 

ZS  =  rng-l 

C  CALCULATE  cav  for  pump*usb*anl>  signal  waves. 

CALL  CAVSUB 
PRINT  2I*RNG*CSIG*CAv 
131  CONTINUE 

END 


C 


SUBROUTINE  CAVSUB 


THIS  SUBROUTINE  CALCULATES  THE  COEFFICIENT  OF  AMPLITUDE 

variation  for  pump,signal»and  usr  waves*  and  Returns 

THEIR  VALUES  TO  PROGRAM  CAVANA. 

REAL  KP.KS»PHI,KT,KM,L»LF,MUP.MUS 
COMMON/A/  KPtKS»PHl»KT,KM#L.LF»HUP,MUS*ZS.RNG.Pl 
COMMON/B/  CAV,CSIG*CPUMP 
C  VARY  TURBULENCE  PARAMETERS  WITH  POSITION  IN  TANK. 

C  AUXpi ,AUXP2* AUXP3  ARE  BOOK-KEEPING  VARIABLES.  km 
C  AND  KT  ARE  SPATIAL  WAVENUMBERS.  PhI  IS  A  TURBijLFNcE 
C  PARAMETER*  AS  ARE  CNS  AND  CNP*  WHICH  ARE  FOR  SIGWAL 
C  AND  PUMP  WAVES*  RESPECTIVELY. 

C 

C  VARy  mu  with  position  in  tank 

MUP=(1. 5249*( L/2. *0.1) *0 . 83?3 ) *2. IE-5 
MUS* (1,5249* ( (ZS*L>/2.*0.i ) *0 .8323 ) *2. 1F-K 
AUXP1*0.105* (RNG/2.) 

AUXP2*1»695-  ! RNG/2 • ) 

AUXP3»(  12.44*  (l. /AUXPI  )♦  (i./AljXP2>  ) 
KM*(PI/2.)*AUXP3 

PHllwPHl* (4,7854) *KM** (-0.5) 

KT1*KT*SQRT (KM/22.9) 
CNS*MUS/SQRT(6.33«PI*KT1**(-2./3.) ) 

AUXPl«o*lo5^  <L/2.) 

AUXP2«1.*95-(l/2.) 

AUXP3«(12.44*(1,/aUXP1 )*(i./AuXP2) ) 
KM*(PI/2.)*AUXP3 
KT1»KT*SQRT (KM/22,9) 
CNPaMUP/SQRT(‘.33«Pl*KTl**(-2./3.) ) 

c  calculate  coefficents  for  fluctuation  terms. 

A1  *  0,3l*CNPW*2,*KP**(7./6.) 

A2  *  0.13*CNS**2.»KS**(7./6 J 

A3  *  (PI/480, )*PHI1*(KT1**4.-kM*»4.) 

C  CALCULATE  THE  LOW  FREQ  COMPONENT*  lf*  of  the 


c  amplitude  fluctuations  in  the  usb  wave# 

Tl«A2*ZS**(li ,/6.) 

T2*3.*a3*ZS**3, 

A»T 1 *T2 

T3f (11./6, )*A2*ZS«*<5./6.) 

T4«3.*A3*ZS**?. 

B*T3*T4 

T5*<A^B*L)**(3./2.)-A**(3./r.) 
LF*<1./L**2.)#(4./(9.*B##?.j )*T5**2. 

C  CALCULATE  CAV  FOR  THE  PUMP  WAVE.  An  APERTURE  FACTOR 

c  is  used  to  account  for  fjnite  size  of  hydrophone. 

CPUMP  =  0.31*CNP**2.*KP*#(7,/A.>*L**<U./a.  J 
HF*<0.A922/0.31)*CPuMP 

cavsq  *  lf  ♦  hf 

APFAC1  *  1.  -  (9.9613E*2/S0RT(U » 

APFAC2  *  !•  -  I1,oAA7E-1/SQRT(L) ) 

CAV  «  SQRT<CAVSQ*APFAC2) 

CPUMP  *  SQRT(CPUMP*APFAC1) 

CSIG  =  0.l3*CNS**2.*KS**(7./6,)*RNG**(l|./6.) 

CSIG  *  CSI6* ( (PI/48O.) *PHI1* (kT1**4.-kM##a. ) 

*RNS*«3.) 

CSIG  =  SQRT(CSIG) 

return 


no 


C 
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PROGRAM  CAVNUM  ( INPUT .OUTPUT ) 

£•******#**»♦**••<»***«*<»*«*#***  ***************************** 
c 

c  This  program  calculates  the  mean  square  amplitude 

C  FLUCTUATIONS  in  the  SIOEbANd  PRESSURE  FOR  a  PARAMETRIC 
C  RECEIVER,  WHERE  THE  HYDRcPHqNE  15  iN  THE  NEARFlEt D  OF 
c  THE  PUMP  TRANSDUCER.  A  NUMERICAL  INTEGRATION  IS  USED. 

c 

C  This  VERSION  OF  the  PROGRAM  IS  used  TO  COMPUTF  CaVPR 
c  FOR  AN  ARRAY  LENGTH  Le0.lM  AND  VARIABLE  RANGE.  SIMILAR 
C  VERSIONS  OF  CAVNUM  ARE  USED  FOR  L*i *SM  AND  VARIABLE 
C  RANGE,  AND  FOR  FIXED  RANGE  AND  VARIABLE  ARRAY  LENGTH. 

C 

C  SET  up  INITIAL  VALUES  and  CONSTANTS.  KS,KP,AND  KIISB  aRE 
C  SlGNALfPUMP,  and  USB  wavenumbers.  AL  IS  ARRAY  lfngTh. 

C  CNR  AND  MUP  ARE  TURBULENCE  PARAMETrRS  IN  THE  InTfRaCTION 

c  region.  cns»mu,km,kt,  and  phi  are  turbulence  parameters 
c  in  the  centre  of  The  tank. 
c 

REAL  L,KS,KP,KUSB,KMfKT,MU,vUp 
REAL  IMKIi,JMKI2 

COMPLEX  Cl »C?*VUSB»KZ»AUXT, aUX2»AUX3,RInT 
COMMON/ A/ aL » ,2i Ilf) tZ?<l * > 
common/b/cnp,kp*cnS»kS,phi »kT,KM,PI »MU* ZS. A* AP 
COMMON/C /kT 1 ,PHI1 ,A) , A2 » A3 » A4 * A5 » L » I t J* K »m, n*MUP*RBL 

COMMON/D/kUSB,A1I »AUSB, AS,A4A»B1B2jS1S2,HIS2»S1B2 
01  MENS  I  ON  Y1  (10)  »  $UMZ  1  (In)  iSUmRI  U<i)  ,SRl  <M>)»$IM<IO) 
DIMENSION  Z(l')  ,RnG<!M 
DATA  CNP.KP  /I ,16?E-4»^.16RiE4/ 

DATA  CNS.KS  /1 .9594E-4,4.1601e3/ 

DATA  PhI,kT,KM  /2.8E-12, 1 7P. 0,22.9/ 

PI  =  3.U159265 
KUSB«4.586E4 

C  FILL  ARRAYS  WITH  /EROS  AS  INITIAL  vAlUES. 
data  SUMRI.AI  ,z  730*0.0/ 

DATA  Z1 ,Z?* Yl,SUMZl/4  *0./ 

C  AR  AND  HR  ARE  THE  CHARACTERISTIC  RaDiI  OF  THE  PUmP 
c  transducer  and  hyorophone*respectivEly. 

AR*0 • 005 
HR* I » 12&E-3 

c  l-RlIm  is  the  upper  limit  of  Integration. 

RlIM*0.001 

c  ap  is  the  effective  patch  sjze. 

AP*1.4IE-P 

C  Z ( A ) »/l ( 1 >  »Z? ( 1 )  ARE  LOWER  LIMITS  OF  INTEGRATION. 

Z(D»i.E-e 
Zl  m*l.E-8 
ZZU)*l,E-8 

c  print  heaoing  fop  results. 

PRINT  30 

30  FORMAT ( /,5X, WCPRF , 12X» ARANGfFf I5X » *ZS*»  15X  , 

vary  range  WITH  ARRAY  LENGTH  CONSTANT  at  o.im; 


on  o  o  o  o  n  o r> 


ORNG.0.2 
AlU)=0.1 
RNG <  1  >  *  .? 

00  97  N* 1 » 1 0 

calculate  homogeneous  USB  SIGNAL « VliSB »at  THE 
hydrophone  output. 

H»<AL<li-RLlM)/0, 

L*AL< 1» 

ZS*RnG<N)-L 

vary  TURBULENCE  PARAMETERS  with  position  in  tank 
MU*  ( 1 .5249*  (  (ZS*L)/2.40.1)  ♦ i. •  0323)  *2.  IE-5 

MUP* < 1.5249* (l/2. 40.1 >*0.8323 ) *2. IE-5 
AUXP1*0» lfl5*  <  ZS*L) /2. 

AUXP2*1.695-(ZS*L>/2. 
AUXP3*(12.44*1./AUXP1*\./AUXP?) 

KM?  (P 1/2.  )  *AIJXP3 
PHll*PHl* (4.7854) *KM**  <-n.5> 

KTl*KT*SQRT (KM/22.9) 

CNS*MU/SQRT (o.33*PI*KT1*#(-?./3.) ) 
AlJXPl*0»l05*L/2. 

AUXP2*1 .695-L/2. 

AUXP3=(12.A4.1,/AUXP1 ♦l./AUXP?) 
KM*(PI/2.)*AUXP3 
KTI*KT*SORT (KM/22,9) 

CNP=MUP/SQRT(r‘.33«PI*KTl*»(-2./'3.)  ) 

calculate  coefficents  for  fluctuation  terms. 

A 1  =  0,31«CnP**2,*KP** (7,/6. ) 

A2  *  0.13*CNS**2.*KS**(7./6, ) 

A3  «  (PI/48o.)*PHH*(KT1*»4.-KH**4.) 

A4*2.5E-2*KP»*2*MuP**2 
A5*l .25E-2*KS**2*MU**2 

INTEGRATE  HOMOGENEOUS  PRESSURE  Using  SiMPSONtS  Rule. 
DO  50  1*1.9 

CALCULATE  Cl  AND  C2.  AX.AUX  ARE  AUXILARY  VARIABLFS. 
AX1.=  (KP*(AR**2.)  )/<2,*2<I)  ) 
AUXl*CMPLX<l.?.AXi) 

AX2*(KUSB*(HR**2,) )/(2.*(l-Z<T> ) ) 

AUX2*CMPLX ( 1 • Ot AX2 ) 

AX3*(*.*Z(I)*(L-Z(I>  )  >/<  (HR**?.)*UR**2.n 
C2»AX3*AUXl*AUX2 

AXl* ( (KP**2.)*(AR**2, >)/(♦.« IZ(I)**2.) > 

AUX 1«AX 1 /AUX 1 

AX2*( fKUSB##2.)#lHR**2.) >/(4.*‘L-2(I) )**2.) 
AUX2*AX2/AUX? 

AX3*(KP/(?.*Z<in )*<KUSB/(2.#(L-Z(I) ) ) > 

AUX3*CmPLX<0.»AX3) 

Cl*AUX3*AUXl*AUX2 

c  calculate  Real  and  imaginary  parts  of  integraNd.rinT 
RINT*1./ (Cl*C2) 

SRL(D*REAL(RINT) 

SIM<I)«AImAG(RINT) 

C  INCREMENT  Z(I)  by  h, 

ZIIM)*Z(II*H 
90  CONTINUE 


nriooo  onnnn  ooo 


C  PERFORM  INTEGRATION.  ffl  CUlaTE  rEaL  AND  IMAGINARY  PaRTS 
C  SFPaR a  T  Fi  >  . 

1  ♦  ‘■M  ,  ,  4  ,H  (*»♦  »*SRl  (  7  I  ♦  4  #  *  SRI  (8J4SRL(9>) 

1  *‘-m  <  S  » •>.  .  *4  .  *“  IfMft  >  •'  .•bIM<7»  ♦4.«STM(8>  ♦SIM<9)  ) 

VUSrtsf  ME'l  >  i  W|  ,pi  t) 

DETERMINE  MAGMTtiUfc  ANO  PmaST  OF  VUsB. 

PMAGsC  AbS  I  VI  ic  M  I 

PhASF*AT  AN.  (A  IMA Gf  VUSB)  .REAL ( VUSB)  ) 


INTEGRATE  OVER  ThC  APHAYi  Z).Z2sn  TO  Z1 »Z2* (L*»L TM) 
IN  STEPS  OF  M. 

DO  100  Jsl 
00  101  I«l*« 


CALCULATE  ARGUMENT,**  <Z>  *  OF  HOMOGENEOUS  PRESSUPF 

calculate  C 1 1 C2  FOR  Z 1  • 

AXl=  ( KP#  (  aR«**?  ,))/(2.*Zl(T)  > 

AUXl*CMPl.X«l.  »A<1» 

AX2*  (KUSB«  <HP*«2,  J  >/'<2**(l-?1(1>>> 

AUX2aCMPLX < l.  ’  *AX2» 

Ax3*(A,»Z*  ln*IL*Zl  I  !>)>/(  (hR**2*>*<AR**P.)  > 
C2«AX3*AUX1*AUX2 

AXl  =  (  (KP»*2.)««AR**2.)  )/(4.«<Zl  ) 

AUXlsAXl/AUxi 

AX2* ( <KUSB*#?. ) * <hR*«2. ) ) / U •*  lL*Zl  (!)>•*»•> 
AUX2=AX2/AUX? 

AX3*(KP/(P#«Z1  <  U  ) )  4 (FUSB/ (?•*  TL-Zl <  I ) > I > 
AUX3*CMRLX  <0. » AX3) 

C1«AUX34AU*1*AUX2 

c  calculate  k<z»  por  Z 1  *  „  , 

KZ*CMPLX (COS (PHASE* ♦ -SIN (PHASE) > 

KZ*(KZ>/ <Cl*C2> 

REKI1*REAL<KZ) 

IMKI 1*AIMaG (KZ) 

C  CALCULATE  C 1 ♦ C 2  FOP  22. 

AXl*TFP*(AR**?») )/(2.*Z2(j) ) 

AUX1*CMPLX‘1.;  *AX1> 

AX2*(KUSfl*<HR*#2,»  > / (2.* (l-Z2 1 > 

AUX2*CMP<  X< l.'»AX?> 

AX3* ( 4,*Z?  T J)  * (L-Z2  < Jll ) / ( (rR**2»  > • ( AR**P. > ) 
C2*AX3»AUX1«AUX2 

AX1* ( IKP**2, ) * { AR»*2. ) ) / <4 • #  * Z*  * J) **2#>  > 
AUXI«AX1/AUX1  „  ,  t  , 

AX2*{ (KUS8**2.)#(Hr##2»> ) / ( 4 •* TL-Z2 < J) >*•>.) 
AUX2*AX2/AUX2  ,  ,  _ 

AX3*(KP/(?,«Z2(J)»>4 (AUS0/ ( ?•*  <L»Z2 l  J) J  )  » 
AUX3«CMPLX  <0.»AX3) 

Cl  ■AuX34AlJXl  *  AUX2 


oo  ooo  non  n  n  n  n  n  non  non  non  non 
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C  CALCULATE  K<Z>  FOR  Z2. 

KZ*CMPLX(c0S<PHASE>  t-SlNtPHASE) ) 

K7«(KZ»/<C1*C21 

RFKI2sREAL  <KZ) 

IMKI2*AIMA6!KZ) 

c 


CALCULATE  CONTRIBUTION  OF  AmPLITuDF  flucts 

CALL  B1B2SR 
BB*BlB2*REKli*REKl2 

calculate  contribution  of  Phase  fllCts 

CALL  S1S2SR 
SS*SlS2*IHKli#IMKl2 

CALCULATE  CONTRIBUTION  OF  AmP-PhaSf  CORR 

CALL  81S2SR 
BS»REKI1*IMKI2*B1S2 

calculate  contribution  of  phase-amp  cross  corr 

CALL  S1B2$R 
SB=REKI2*IMKI1*SIB2 
COMPUTE  INTEGRAND. 

YHI)=B6.BS.$B*SS 

INCREMENT  zi  by  H. 

103  CONTINUE 

Z 1 < I ♦  1 >  =  Zl  (I)  ♦  H 
101  CONTINUE 

INCREMENT  Z2. 

1  ♦Y](5)  #2  ,.4,*Y|  (6>  .?.*VT(7)  .4,*Yl  («)  1  ?9)  > 

Z2<J*1> SZ?<J>  ♦  H 
100  CONTINUE 

INTEGRATE  OVER  VARIABLE  *2. STORE  RpSULT  IN  SUM??. 

SUMZ2*(H/3)*(SUMZ1  m.4*SUMZl  (2)*?#SUMZ1  <3>.4*SUMZl  (♦>♦ 
1  .4,*SUMZl  <8>  ♦  StIMZ  1  (9)  ) 

CALCULATE  RESULT  OF  NUMERICAL  INTEGRATION 

SUMZ2*SUMZ2/RMAG*«2. 

calculate  the  contribution  of  source  wafers  from 

C  Z»!l-RLIM>  TO  Z«L.  THE  HYDROPHONE  IS  IN  THE  ExTrEmE 
C  NEArFIELD  OF  THESE  WAFERS*  so  THFIR  RADIATION  IS  ApPROX- 
c  imated  by  plane  waves  of  infinite  extent. 

Tl«A2*<ZS*L)**(ll,/6.) 

T2*A3*(ZS.L)**3, 

B1 *  <RL IM*«2,/l*#2 J  # ( T1 ♦ T2) 
T3«L*«<23./12.)-(l-RLIH)**(23./12.) 


14(1 


T4*0.2722*T3**2.  ».j«.2722*R1  I***  123./6) 
T4*T**2*1,27?2*PL  TM«*  <23. /  I  2.  >*T3**?. 

02*  <  A 1 /L **2 . 1 *T4 
0SQ*B1 »B2 

C  SUM  TmE  NUMERICAL  AND  ANALYTICAL  RrSuLTS, 

CPRSQ=SUMZ2*BS0 

Q 

C  CORRECT  FOR  THE  EFFECT  of  FINITE  APERTURE  AT  The  HyROPHONE. 

c 

APFAC  *1 <0.10423/SqRT <L/2. ) ) 

CPRSQ*CPRSU*APFAC 

c 

c  calculate  the  mean-sq  AMPLITUDE  flucts.cpr.in  the 

C  SIOEBaNO  PRESSURE  aT  7*L.  PRINT  ThE  RESULT. 

c 

CPR*SORT (cPRSO) 

PRINT  31*  CPR,RNG(N> »ZS,L 
31  FORMAT (/«E12»4»F14»4»7*«Fi ';.*»F1<>»*> 

c  increment  the  Range. 

79  CONTINUE 

RNG  (NM  )  =PNG  (N)  ♦DRNG 
97  CONTINUE 

END 
C 


C  ENO  OF  MAIN  PROGRAM 

C*********************************************************** 

c 
c 

C  CONTRIBUTIONS  FR OM  THE  CORRELATION  TERMS  B1B2»R1«;?.S1B2»S1S2. 

c 

c  .  . 

£•******#*#*••***********•*•********** ********************** 

SUBROUTINE  BfBZSR  m 

C****************** ************* ************************** 

c 

C  CALCULATE  BlBP»  THE  AMPLITUDE  FLUCT  CqRR  TERM, 

C 

REAL  L»KS»KP,KUSB,KM,KT ,MU»mUP 
real  LCORJ »LCOR2»LCOR3»LCOR4 
COMMON/A/AU2n)tZi<lC)  *  Z?  f  1  -  > 

COMMON/B/CNP.KPfCNS*KS*PHI»KTtKM*Pl»MU*ZS,A*AP 
COMMON/C/KTl ,PHI1 ,Al»A2,A3»AA*A5»L*It J*K»m,n*MUP*RBL 
COMMON/D/KUSB»All,AU5B»AS» A4A*B1B2»S1S2*BTS2*S1B2 
A l  *  0,3l*CNP**2,«KP**(7./6. » 

AUSB*0 • 31 *CNP**2 ,*KUSB** (7./6. ) 

All*SQRT (Al*AUSB>  i 

Tl  »  (Z1  (I)*(AHi>-Z2(JM  )**«]  1./12.)  L| 

T2  *  <Z2<J)*(AL<M-Zl<lm*#(ijl*/1?.>  1 

T3  *  (Zl(I>*Z2(JM**ni./i2,>  1 

T4  *  UAL(I)-Z1<I)>#(AL(i)-Z2(J)>)**(iI,/T2.»  1 


on  o  onoooo 
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ACCOUNT  FOR  LONGITUDINAL  CORRELATION  OF  HIGH  FrEQ  FLUCTS 
TmLRE  are  FOUR  SFPARaTE  terms  for  THE  DIFFERENT 
geometries,  also  there  are  Two  VERSIONS  of  Each 
Term  FOR  zi  GT  Zz  or  Z2  RT  z I. 

IF  1 22  c  J »  .GT.zi  (I)  )60  TO 
DELZ1*Z1<I)-Z2<J> 

LC0Rl*<l.*3.#0ELZl>''<2.«Z?(j>  ) 

LCORj  *LCOP 1 / ( (I.*DELZ1/Z2(J) >*«<3./2.) > 

0ELZ2*  AL ( I ) »Z2 IJ) 

LC0R2*(1 .♦3.*oELZ2»/(2.«Z? ( j) » 

LC0R2=LC0R2/ ( { 1 •♦DELZ2/Z2 ( J) )**(3./2,) ) 

DELZ3»Z1 ( I ) -Z2 <J> 

LC0R3*(1 .♦3.«dELZ3>/(2.*(aL (1 »-Zl <I) ) ) 

LC0R3=LC0R3/ ( l J . ♦DECZ3/ (  AL ( j ) -Zl U ) ) ) ** (3./2, )  I 
OELZ*sALn  )-Zl  <I> 

LC0R4*(I..3,#dELZA»/(2.#zT  C I > ) 

LC0R4*LC0RA/( (1 .♦0ELZ4/Z1  ( I ) > ** ( 3 . /2 . )) 

GO  TO  2  2 
201  CONTINUE 

DELZl*Z2< J)-ZJ (I) 

LCORl*U  .♦3.*oELZl )  /  (2**ZM  I  M 
LCORl*LCORl/( (1.*0ELZ1/Z1 ( I ) > ** < 3. /2. > ) 

DELZ2*  AL  < 1 > -Zl < I ) 

LC0R2-U.*3.#0ELZ2>/(2**Zi  (I) ) 

LC0R2«LC0R2/(  n.*0EEZ2/Zi(I)  )**(3./2.U 
DElZ3*Z2(J>-Zl  <I) 

LCOR3*n.*3.*t>FLZ3)/<2.*tAL(l>’z2tJ>  )  » 

LCOR3=LC0P3/  (  ( 1  .♦0ELZ3/  ( AL  (  t  >  -Z2  ( J)  )  )  **  ( 3./,»«  )  ) 
0ELZ4*AL(] )-Z?( J) 

LCORA-U.^S.^OELZAJ/fH.^Zpt  j)  ) 

LC0RA*LC0R4/ ( ( 1 .*0ELZ4/Z2 ( J) >  ** (3»/2. ) > 

>02  CONTINUE 

OELZ*A0S(Z1  (  U-Z2  (J)  > 

RBL«l./<i.*< (2.«DELZ>/<KS*AP**2.>  >*#2.> 
TBl*All*(Tl»LC0R2*T2*LC0R4) *Ai *T3*LCORl ♦ AtlSB*T4*LC0R3 
CALCULATE  LOW  FREQ  AMP  FLUCTS 
A2  *  0.13*CNS**2.«KS#*<7./6,> 

A3  *  (PI/4B0,)*PHI1*(KT1**4.-kM**4.> 

T1  «  A2*(ZS*Zl(I)>**ni./6.) 

T2  *  A3*  <  Z  S ♦ Z 1 (I > )##3« 

T3  *  A2*(ZS»Z2<Jn#*<n./f,.) 

T4  *  A3*(ZS*Z?( J) >**3. 

TB2  *  SQPT<  <T1*T2>*<T3*T4> ) 

TB2*TB2*RBL 
B1 B2»T81 ♦TB2 
return 
END 


noon 
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C 


SUBROUTINE  S1S2SR 


calculate  siS2»  the  phase  Fluctuation  corr.  term. 


c 


203 


204 


C 


REAL  L»KS*KP#kUSB,KM,KT.Mu»mUP 
REAL  LCORl  •LCOR2.|ECOR3tLCOR4 
common/A/al ( 25 >,zi  do)  »Z2(r  ) 

COMMON/B/cNP,KP»CNStKS»PHI»KTtKM#Pl ♦MUt ZS.A.AP 
COMMON/C/KTl,PHIl,AI *A2*A3»a4,A5»L»I»J»K*m*n»MUP*RBL 
COMMON/0/kUSB»AU*A  SB»AS.A44,BlB2»SlS2.BlS?iSlB2 
A4*2.5E-2*KP##2*MuP#*2 
AS*2.5E-2*KusB*#2.#MuP#*?. 

A44*SQRT (A4#AS) 

T1«SqRT<Z1  < I ) * (L-Z2  <  J)  I  ) 

T2«SQRT(Z?(J)*(L-Z1  U)  > ) 

T3*S0RT(Z1  (n*Z2(J)  ) 

T4*SQRT< (L-zT(I> >«(L-Z2<J> ) ) 


ACCOUNT  FOR  LONGITUDINAL  CORRELATION  OF  HIGH  FREQ  FLUCTS 
IFIZ2(J) .GT.Z1 (I))60  TO  2n3 


DELZl=ZI (I 
LCORl*l • ♦ ( 
LCORl*  I./ 
DELZ2*L-Z? 
LCOR2*l.M 
LC0R2*  1./ 
0ELZ3*Z1 (I 
LCOR3*l • ♦ ( 
LCOR3*  1./ 
DELZ4*L-Z1 
LC0R4*1 • ♦ ( 
LC0R4*  1./ 
GO  TO  204 
CONTINUE 
OELZl=Z2< J 
LCORl*l ( 
LCORl*  1./ 
DELZ2*L-Zl 
LC0R2*1 .♦ ( 
LCOR2*  1./ 
DELZ3=Z2( J 
LC0R3*1 C 
LCOR3*  1./ 
0ELZ4*L-Z? 
LC0R4*1 • ♦ ( 
LC0R4*  1./ 


-Z2CJ) 

2.*DELZl ) / ( KP*AP«*P ) ) ##2 
LCORl*SQ«T(l .♦oELZtVZ2(J) > ) 

J) 

2.*DELZ2W  (KP*AP**?>  >**2 
LCOR2*SQMT(1.*DElZ?/Z2(J) ) ) 
-Z?CJ) 

2.*0ELZ3>/<KP*AP«**?)  >**2 
LC0R3*SQ«T(1 .♦DELZV<L-Z1 <I> J ) l 
I) 

2.*0ELZ4)/(KP*AP«*P) )**2 
LCOR4*SQRT< 1,*DELZ4/Z1<I)  )  ) 


-Zl <  1) 

2.*0ELZ1>/<KP*AP*#?> ) *#2 
LCOR1*SQ«T(1.*dELZT/ZHI)  >  > 

I) 

2.«DELZ2>/(KP«AP*#P) >  *#2 
LCOR2*SQPT<t .♦oelZp/zi <I) ) ) 

-Zl Cl) 

2,*0ELZ3>/<KP#AP**?> )**2 
LC093*SQPT(1  .*nELZ3/<l.-Z2<j)  )  >  ) 

J) 

2«*DELZ4>/(KP*aP*#?) )##? 
LC0R4»SQRT(1 .♦0ELZ4/Z2C j) ) ) 


CONTINUE 

SSHI*A44*(T]*lCOR2*T2«LCOR4) ♦A4*T3*LC0Rl 


♦A4»T3*LCCR1 

CALCULATE  LOW  FReQ  PhASe  FLUcTS 
A5«1.?5E-2#KS**2*PU**2 

T5*A5*SQRT< (ZS*Z1 (I) )*IZS*Z?CJ>  >  > 
T5*T5*RBL 

D1«(4.*ZS«Z] (l))/(KS*AP**p,j 
D2*<4.«ZS*Z2(J) >/<KS#AP**?.) 


k 


o  o  o  o  <2  no  o  or>o 
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T5*SQRT ( ( l««01*ATiN<Dl I ) * (1 . ♦0£*ATAN ( D2> )i*rS 

SlS2  =  S<»r’I*T5 

RETURN 

end 

c 

c 


SUBROUTINE  B1S2SR 


,  v  <•  «  * « o «  )9<,n  »»#***#* 


CALCULATE  B  l  S2 »  Th£  AMPl  -HhaSF  F  u  ‘  • ...  '•  0^^'*'  «»;mPR 

MEAL  L  #KS*NP»KUS8»KM,KT ,mh,hJP 
COMMON/ A  /  AL  (?f,) 

COMMON/B/CNPfKPtCN^tKSfPHj  ,  K  T  ,  KM  t  P  I  «  «Hi*  ;  S  .  A  «  AP 

COMMON/C/K  ri  ,PHI1  5  A1  t  A2  ,  AT»  a  A  ,  A£>»L*  I  *  I  ♦  K  »  M  f  n  *  MUP  »  PBL 

COMMON/u/kOSB  ,  A 1 1  ,  AU$B»  AS,  A4*  .  Pi 82,  SI  $2,  Ml  lj}2 

TBlsSORT ( A2* ( ZS* / i t I > >  ** ( li ,/B* > *A3»  (ZS<  ?i  ,  I ) ) *«3. I 

TS  Is  SORT  ( AS*  <  7  $«/;>{  J)  )  * 

n?s  (4»*ZS*Z2  ( J>>  /  {hS*AP**?#i 

TSI*TS1*SQWT  ( 1  ,  ♦D.^AfANilD*)  ) 

T 1  *  TB1*  T  SI 

tb2=sort  i  Aivi  <  i )  «•»  <  n .  /*.  t  > 

T52-S0PT  U4«;?(j|| 

T2*Tri *Th? 

TB3»SQRT  <  A  l*7l  U  ••  **  <H  .  'h.  )  ) 

TS3-SQRT  (  aS»  0.-/2  (,.>  |  ) 

TJ=T83*TS3 

TB4*SQRT  (  AUSb*  U.-2)  (  I )  >  **  <1  1  *  /t>. >  ) 

TS**SQRT  (  A**Z2  < Ji  \ 

ta*tb4»tsa 

TB5*SQRT  ( AUSR*  <L“/1  <  1 1  >**(11  •  /6«  >  > 

TS5  =  S0RT (AS*(L-Z?( J>  >  > 

T*»*TB5*TS5 

B1S2bTi*T?«T3*T4»T5 

ACCOUNT  for  LONGITUDINAL  CORRELATirN  OF  FiUffS 

BlS2«8IS2*0,n 

return 

cno 


SUBROUTINE  S1R2SP 

!#♦# *•#«*« •#•**#*##*#*#*•• 


CALCULATE  S 1 B2  *  THE  AMPL ITUOE-PHASE  Fl.UCT  CROSS-Corr 

REAL  L»^S,KP,KUSB,KM^T,MU,f‘Up 
common/a/al ( )  ,zinr>,z?u'> 


no 


COMMONi/H/rNP.KPlCNS*Kci#r-t,)  '  ,KM,P!  ,Mi  #ZS.A,AP 

fOMMON/C/r-.  .  A  i  i  :  . . ’  .  *  .■ 

t  OMMON/U/I*  !Sfi»A]J,Ai>f,hiA  ,  ri  1  R?  *  c  1  be 

TBl^SJRT  (  A2#  <  2  W?  ‘  Ok  ;  „  i  ,  .  *  ♦  -  *  ■  -  *•  .■  /?  ( ,i>  \  #*3.  ) 

TS1 “SORT t Ab« (2S*2 j «I )  >  ) 

Dl» (*.*^S*Zl U ) ) / (KS*AP*»>.  . 

TSl*TSl*S««T  n  .♦DI*ATA'’HP!  •  1 
Tl*TBi«TSl 

TB?=SQPT  (  Al*z?  1  !./,-»  ) 

T  S2-SOR  T ( A4*Zi  (  I)  ) 

*  re  i « if?? 

TR3=S3RT  <AUSH»  <L*z2(jM»«'.ji..i  .  »  ‘ 

Ts3*SOPT  (A4«7t  mi 
Ti*TB3*TS3 

T B ASSORT  ;  r.<)  **  t  s  2  »  'f  , 

T5-.A-SOR  T  (  ar,“ 

T4*  f  H4 *  T S4 

rpS  =  SOWT  I  AUSH<»  (U»z2  o  -  I  I 

2  <;5s:'-  rjRT  •  AS«*  U.-/1  ( IM  I 
I  C*T  B-*T  55 
bl B2*Ti ♦T?*T3+Tft*T5 

ACCOOf.r  for  I  ONG  I  Tt.>D  |NAL  COpRti  A  2  ■  ■  r  •  .  T- 

SlSl’af  lB*>»in6 
Rf  TtiRN 


APi.  ENDIX  6.  ANALYSIS  OF  THE  NEARFIELD  PARAMETRIC  RECEIVER 

In  developing  a  theoretical  analysis  for  the  nearfield  parametric 
receiver,  it  was  assumed  in  Section  4.2  that  the  pump  wave  could  be  approxi¬ 
mated  as  a  collimated  plane  wave.  In  evaluating  the  integrals  arising  from 
this  assumption  (see  Appendix  2),  no  account  was  made  of  the  effect  of  the 
finite  radius  of  the  pump  transducer.  Consequently,  the  phase  of  the  inte¬ 
grand  of  Eq.  (4.19)  was  such  that  only  the  amplitude  correlation  term 

2 

contributed  to  the  amplitude  fluctuations  <B  >  of  the  second-order  pressure 

PK 

wave,  as  shown  in  Eq.  (4.20).  A  more  detailed  analysis,  based  on  the  work 

of  Chotiros,^*^  is  presented  in  this  appendix.  It  will  be  shown  that  there 

2 

are  additional  contributions  to  <Bpp  when  the  finite  radius  of  the  pump 
transducer  is  taken  into  <  c count. 

A6.1  Second-Order  S > lution  for  a  Homogeneous  Medium 
The  geometry  for  thi.-  analysis  is  shown  in  Fig.  A6.1.  As  in  the 
analysis  of  Chapter  4,  tl  source  of  the  signaj.  to  be  detected  by  the  para¬ 
metric  receiver  is  located  at  z  =  -Z  .  The  signal  source  is  situated  suf¬ 
ficiently  far  from  the  pump  that  the  signal  pressure  wave,  pg,  is  assumed 
to  be  planar  in  the  vicinity  of  the  parametric  receiver.  Omitting  the  time 
dependence,  the  signal  wave,  may  be  represented  by 

-jk  z 

p  (z)  -  P  e  S  ,  (A6.1) 

*s  s 

where  the  symbols  are  defined  in  connection  with  Eq.  (4.1). 

The  pump  transducer  is  approximated  by  a  plane  wave  having  a  Gaussian 
shaded  transverse  amplitude  function;  i.e.,  the  pressure  amplitude  at  a 
point  (p,9p)  on  the  face  of  the  pump  transducer  is  assumed  to  be  given  by 

- (p2/b2) 

Pp(p,Op)  =  Pq  e  ,  (A6.2) 
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where 


P  and  0^  are  coordinates  on  the  face  ot  the  transducer,  as  shown  in 
Fig.  A6.1, 

Pq  is  the  pressure  amplitude  at  the  center  of  the  transducer,  and 
b  is  the  characteristic  radius  of  the  transducer. 

The  pump  signal  at  some  point  (r,0+,z)  in  the  interaction  region  is  the  sum 
of  contributions  from  all  elemental  surfaces  on  the  face  of  the  piston,  and 
can  be  shown  to  be 


Pp(r.0+.z)  =  p 


2it  " 

//' [( )  (  '  •')] 

n  n  '  ' J 


du  p  dH 
P 


(A6.3) 


2  2 

where  R  =  r  +  p  -  rpcos(0,  -  0  ),  3nd  it  is  assumed  that  R  <<  z. 

±  P 

The  second-order  pressure  at  a  point  (e,8^,L)  on  the  face  of  the 
hydrophone  can  be  found  from  Eq .  (2.7)  to  be 

2  ft  1. 

p  (0,bl,L)  =  I  I  I  - r  q  .  (r , u  .  , z) 


3±(a,0h,L)  ~  J f J  4ML  -  z)  qr(r,V 

o  o  o 

‘j-JkjjL  -  z 


2arr“z)*j/rdr  dVz 


(A6.4) 


where 


2tt 


q+  = 


jBP  P  to 

sot 

2  4 

P  c 
o  o 


2  ,  2 


-iV  f  f  - „2/b 2  f-£Vl  , 

J  J  e  exp |  2z  f  d"+  pdp  ’ 


and 


o  o 


R'  =  r  +o  -  2rocos(0^  -  0+) 

It  is  also  assumed  that  the  voltage  response  of  the  hydrophone  has 
Gaussian  shading  in  the  transverse  dimension,  so  that 


M(o  ,  0^)  =  exp(-oVh^) 


where  M  is  the  voltage  response  in  V/Pa,  and 
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Mq  is  the  voltage  response  at  the  center  of  the  hydrophone. 

The  voltage  produced  at  the  hydrophone  by  the  second-order  pressure  at 
z  =  L  will  therefore  be 


2  2 

exp(-o  /h  )  odo  d6 


(A6.5) 


Substitution  of  Eq.  (A6.4)  for  p+  allows  the  hydrophone  voltage  to  be 
written  as 


®  2  it  L 


rkp  ,  kt 

2z  2 (L  -  z) 


V“  =  A/// ^  eXPfj 

o  o  o 

M  ftp  P  -ik  L 
1  .  ±  o  s  o  J  + 

where  A  =  - - —  e  "  » 

4wp  c 
o  o 

h-fr  ’«pRs)2 *lHsesatt 


^rdr  d0  +  dz 


(A6.6) 


■ff 


-  2rocos$ 


2(L  -  z) 


pdp  diji 


odo  d$ 


^  =  0  -  9  ,  and 

±  P 


*  =  0*  -  V 

Equation  (A6.6)  can  be  greatly  simplified  with  the  aid  of  tabulated 
Integrals,  and  by  using  the  definition  of  the  Bessel  function.  The  angular 
integration  of  I  and  I2  may  be  written  in  the  form 


2tt 


Scosu 


du 


(A6.7) 


where 


I  kprp/z 


I  k+ro/(L  -  z) 


for  I 
for  I2- 
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The  definition  of  the  Bessel  function  of  the  first  kind  is 


111 


Vs) 


j ( 3s i nx  -  nx) 


dx 


Using  the  change  of  variables,  u  =  x  -  -j  ,  and  taking  the  case  for 
n  =  0,  Eq.  (A6.8)  becomes 


J0(B) 


1  /  IdcOKU  , 

*  *  J  “ 


3n 

2 


From  Eqs.  (A6.7)  and  (A6.9)  we  can  write 


I  =  2irJ  (6) 
o 


so  that  1^  and  become 


1  =  2*  /  J0(^-)expf-(§)2-jkp,2/2x] 
o 

f  Jo(L^l)eXpt'(h)2_jk±o2/2(L  “  2)1  0dj 


z]  pdp  , 


and 


I2  =  2tt 


Now 


n: 


and  I2  are  in  the  form  of  a  tabulated  “  integral,  viz. 


C 

/ 


xv  +  1  e-aX2  Jv(8x)dx  =  6V(2a)"(v  +  1}  exp(-62/4a) 


Evaluation  of  1^  and  I2  using  Eq.  (A6.10)  gives 
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and 


T  =  - 

TTh2 

-  exp' 

L  ( v  Y 

h2 

jk±h2 

\L  ■  V  4 

■l  +  jk*  1 

1  +  2(L-z) 

h2  +  2 (L  -  z) 

(A6.8) 

which 

(A6.9) 


(A6.10) 

(A6.ll) 

(A6.12) 
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By  substituting  Eqs.  (Ab.il)  and  (A6.12)  into  Eq.  (A6.6),  the  hydrophone 
voltage  can  be  expressed  as 


VH  =  A± 


L  “  _r 

fjt  •  1 


2rdr  dz 


o  o 


where  A,  = 


™±*PsPoMo  -Jk>L 
e 


±  >  A 

2p  c 
o  o 


c  =  j  /— E  -H  _ _ \ 

C1  J\2z  +  2 (L  -  z)J 


i  2.  2 

k  b 


2  2 
k  h 


(A6.13) 


/  jkpb  \  /  kX 

4|l  +  — - )  All  +  j  — - — |(L 


2z 


2(L  -  z), 


z)' 


k  b 


k ,  h 


and  C2  =  4z(L  -  z)  fl  +  j  -^-\h  j  ^7 


As  and  C 2  are  independent  of  r,  the  integration  with  respect  to 

2 

r  can  be  done  by  the  change  of  variable,  u  =  r  .  The  result  is 


■U 

VH  =  A±  J" C^C 


dz 


(A6.14) 


Equation  (A6.14)  represents  the  hydrophone  voltage  produced  by  the  second- 
order  pressure  field  in  a  homogeneous  medium.  In  the  next  section  this 
result  will  be  used  to  determine  the  amplitude  fluctuations  in  the  second- 
order  pressure  when  inhomogeneities  are  present  in  the  medium. 

A6.2  Amplitude  Fluctuations  in  an  Inhomogeneous  Medium 
In  an  inhomogeneous  medium  the  hydrophone  voltage  will  fluctuate  in 
amplitude  and  phase,  and  may  be  written  in  a  form  similar  to  Eq.  (4.14)  ,  viz. 


V±  =  VH^1  +  bPR  ~  ^ SPR^ 


(A6.15) 


where  BnD  and  SDD  are  the  amplitude  and  phase  fluctuations,  respectively, 
r  K  rK 


of  the  interaction  frequency  pressure  at  the  hydrophone.  A  development 
similar  to  that  in  Section  4.2  gives  the  following  expression  for  v+  in 
terms  of  the  total  amplitude  and  phase  fluctuations,  B  and  S: 


‘  '  f 


(1  +  B  -  jS)  -—r-  dz 
L1  z 


(A6.16) 


Following  the  methods  used  in  Section  4.2,  we  substitute  Eq.  (A6.16)  for  v+ 

in  Eq .  (A6.15),  and  at  the  same  time  write  v  in  terms  of  its  magnitude  and 

H 

phase: 


I VH I eJ  9  (1 


+  BPR  -  JSPR>  =  A±  j 


(1  +  B  -  jS) 


(A6.17) 


Subtracting  the  homogeneous  solution,  Eq.  (A6.14),  from  both  sides  of  Eq. 

(A6.17),  and  dividing  by  the  hydrophone  voltage  v  ,  we  obtain 

H 


jSpR  =  T^T  /(B " jS)  e"j9  c^f: 


0  L 

'«!  lJ <B  -  J 


jS)  K(z)  dz 


(A6.18) 


where  K(z)  = 


The  amplitude  fluctuation  term  can  be  found  by  taking  the  real  part  of  Eq. 
(A6.18),  and  is 


Kr1/ 


[B  Re (K)  +  S  Im(K) ]  dz 


The  mean-squared  amplitude  fluctuations  of  the  second-order  pressure  at  the 
hydrophone  are  therefore 
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<BPr" 


I.  1. 

VH'  if  <(Bl  +  Sl  Im(Kl 


)] 


o  o 


x  [B,  Re(K0)  +  S9  Im(K2)]>  dz  dz2 


Writing  this  result  in  terms  of  the  coefficient  of  amplitude  variation,  and 
expanding  the  terms  in  angular  brackets,  we  obtain 

L  L 


CAV 


PR 


=  |vH,-2//  /ReO^)  Re(K2)  <BXB2> 


+  ImCKj)  Im(K2>  +  ReO^)  Im(K2)  <B]lS2> 


+  ImCKj)  Re(K2)  <S 


lV}  dZjdz, 


which  appears  in  the  text  as  Eq .  (5.12). 
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